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1.0 SUMMARY 


An analytical model has been developed to predict the trans- 
mission of propeller noise into the cabin of a high-speed 
propeller-driven airplane. The model is then used to determine 
the noise control treatments required to achieve a goal of an 
interior A-weighted sound level of 80 dB for three study air- 
planes with different fuselage diameters but with a common 
cruise Mach number of 0.8 at an altitude of 30,000 ft. 

The analytical model describes the excitation in terms of a 
propagating pressure field with a pres.sure amplitude which 
decays rapidly with distance. The fuselage structure is re- 
presented as a series of curved, orthotropic plates with 
in-plane loads, the orthotropic characteristics resulting from 
the effects of smearing out the stiffness and mass of the 
longitudinal and circumferential stiffeners. Finally, the 
receiving cavity or cabin is basically a rectangular parallel- 
piped with ’deformed' curved surfaces. Sound transmission 
into the cabin is then calculated using the concept of acous- 
tic power flow with allowance being made for resonant and non- 
resonant response. The influence of the cabin sidewall treat- 
ment is introduced as an additional transmission loss experi- 
enced by the transmitted power. Sound absorption in the 
cabin results in power dissipation losses which balance the 
net inflow of acoustic power. 

Three baseline aircraft with -wide-body , narrow-body and small- 
diameter fuselages, respectively, are selected based on exist- 
ing airplane designs, and the analytical model is used to 
calculate the propeller noise levels in the passenger cabin. 

It is found that noise reductions of about 25 dB, in the fre- 
quency range of ^.30-230 Hz, are required to achieve the goal of 


80 dBA. Noise reduction methods are investigated in the form 
of either add-on treatments, which can be applied to existing 
fuselage structures, or advanced techniques, which require 
changes to the primary structures. For the add-on methods, 
only a double-wall sidewall with a heavy, limp trim panel 
isolated from the fuselage structure is found to give the 
required noise reduction. No advanced technique is found 
which, on its own, can provide the required additional 25 dB 
transmission loss over the frequency range of interest. 

The elements of the analytical model are presented without 
direct validation. Such validation is, of course, highly 
desirable and the design and performance of validation experi- 
ments is discussed. All the test methods considered have 
disadvantages but laboratory tests using a loud speaker with a 
directional horn and a model fuselage structure in the form 
of a cylinder are recommended for validation of the basic 
analytical model and its use in predictions for add-on noise 
control treatments. Tests on present-day propeller-powered 
aircraft could be performed but they would not be directly 
applicable in terms of validating the predictions for the high- 
speed aircraft powered by the high-speed propfan concept for a 

propeller. 
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2.0 INTRODUCTION 


2.1 Noise Control Study 

As part of an overall program to develop aircraft with Improved 
fuel efficiency, consideration has been given to the design 
of high-speed propeller-driven aircraft. Preliminary studies 
[1, 2, 3] have shown that significant fuel savings can be 
achieved by these aircraft but, at the same time, several 
problem areas were Identified. One potentially severe problem 
is that of high interior noise levels in the passenger compart- 
ments. In order to explore this problem more fully, a study 
was initiated to develop an analytical model for the trans- 
mission of propeller noise into an airplane fuselage. The 
model is then to be used to explore potential noise control 
methods for three study aircraft of different fuselage dia- 
meters. The noise control techniques could be add-on treatments 
which would be applied to current day conventional fuselage 
structures, or advanced noise control methods which involve 
design changes to the primary structure. 

The development and application of such an analytical model is 
presented in this report. A general overview of the study is 
provided in this section, while in Section 3 the development 
of the analytical model is described. The model is applied 
to three study aircraft whose aerodynamic, structural, and 
acoustic characteristics are presented in Section 4. Then the 
analytical idealization used to describe the structure of the 
fuselages for the three aircraft is presented in Section 5. 

The following three sections discuss the application of the 
analytical model to the baseline airplanes (Section 6) and 
explore the potential of add-on (Section 7) and advanced noise 
control treatments (Section 8). Finally, Section 9 considers 
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the feasibility and design of experiments to validate the 
analytical model and the predictions for the baseline airplanes. 

The title of this study refers to high-speed propeller-driven 
aircraft. However, the implicit assumption is that a speci- 
fied propeller, called a "propfan" in the relevant technical 
literature (see, for example, [1]) will be considered as the 
propulsive agent. The name propfan will be used in this 
report to identify the high-speed propeller. 

2.2 General Characteristics of Anal ytical Model 

The basic concept of the analytical model is that of power 
flow, and the cabin noise levels are estimated by means of a 
series of power balance equations which equate transmitted and 
dissipated acoustic energy. The model uses as a starting 
point the analytical study performed for the acoustic environ- 
rnent in the payload bay of the space shuttle orbiter vehicle 
at lift-off [4] and the analytical model considers both 
resonant and non-resonant response in one-third octave fre- 
quency bands. Resonant response occurs when the transmitting 
structural modes, or the receiving cavity acoustic modes have 
resonance frequencies in the frequency band of interest. Non- 
resonant response is associated with modes whose resonance 
frequencies lie outside the frequency band of interest. The 
relative importance of resonant and non-resonant response in 
a given frequency band will depend on a number of factors, one 
of the most important being the number of modes with frequencies 
within that band. For example, if there are few or no acoustic 
modes with resonance frequencies within the band of interest, 
then non-resonant acoustic response may be an important factor 
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in determining the total acoustic power flow in that band. 

Thus it might be anticipated that non-resonant response might 
be of particular importance at low frequencies. On the other 
hand, power flow from non-resonant structural modes may control 
the fuselage noise reduction when the radiation efficiency of 
resonant structural modes is very low. 

The analytical model incorporates various simplifications and 
assumptions in order to provide an engineering prediction tool 
which does not require too extensive computation requirements. 
One such simplification involves the division of structural and 
acoustic response into low and high frequency regimes. At very 
low frequencies there are few modes in a given frequency band 
and individual modes must be considered when calculating struc- 
tural response and coupling between structure and acoustic 
cavity. Then as modes in a given frequency bandwidth accumu- 
late as frequency increases, some approximations can be made 
for the expressions for modal coupling. Finally, at high fre- 
quencies it is impractical to consider modal response on a 
mode-by-mode basis and recourse is had to a modification of 
the method of statistical energy analysis. 


2 . 3 Application of the Model to Propeller Aircraft 

The analytical model has been developed in a general sense so 
that it can be applied to propeller noise fields where the 
excitation is discrete frequency in character and has strong 
spatial decay of the pressure amplitude, as well as to cases 
where the excitation is broad-band and homogeneous, such as 
the space shuttle. However when the model is applied to the 
present study aircraft certain assumptions are made to reduce 
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the computational requirements without causing significant 
adverse effects on the accuracy of the calculations. These 
assumptions will be discussed at the appropriate place in. . 
this report, but some of the more important ones are mentioned 

in this section. 

In the analytical model for noise transmission into a propeller- 
driven airplane, resonant and non-resonant response is con- 
sidered in both low and high frequency regimes. However, as 
will be discussed later, it is found that the acoustic power 
flow is dominated by resonant response. Furthermore, because 
the modal density is high, only the high frequency idealisation 
is used, for power flow calculations. 

One assumption concerns the size of the structural unit to be 
considered as a noise transmitting element. While it is 
possible that the fuselage structure can be considered as a 
single unit, it is feasible to break the fuselage down into 
several smaller units. Parametric studies in Section 4 will 
show that above a certain size, the calculated acoustic power 
flow will be essentially independent of the dimensions of the 
structural unit. This is due, to a large extent, to the rapid 
spatial decay of the amplitude of the excitation pressure 
field. Several advantages are obtained by this segmentation 
of the structure, including the ability to estimate the spatial 
variation of sound level inside the passenger cabin. Utilizing 
the concept of structural units, the fuselage is divided into 
several elements in the longitudinal direction. In the 
circumferential direction it is assumed that noise transmission 
through the floor can be neglected and that the transmitting 
structure can be considered as a floor-to- floor or a floor-to- 
ceiling unit, depending on fuselage diameter. 
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In terms of frequency regime, it is assumed that, when the 
modal bandwidth is greater than the separation between modes 
for the (cabin) acoustic volume, the various high frequency 
approximations can be used. Then the acoustic radiation into 
the volume can be described using a modified statistical energy 
analysis approach. 

The basic analytical model estimates noise transmission through 
a single wall structure into a receiving cavity. However a 
typical airplane fuselage contains a sidewall treatment which 
will provide thermal and acoustic control. This sidewall 
treatment can be idealized as a double wall system, the outer 
panel being the skin of the fuselage structure. The space 
between the two panels is filled with insulation material 
such as glass fiber batts. The assumption is made in this analy 
cal model that the additional transmission loss provided by the 
noise control treatment can be modeled as an extra transmission 
coefficient which is applied directly to the calculated power 
flow through the untreated fuselage structure. Whilst this 
assumption is not strictly valid [5], the errors associated with 
the assumption are not large and the simplifications to the 
analytical model are significant. 

In applying the analytical model to an excitation field such 
as is generated by a propeller, consideration has to be given 
to the uncertainties which arise due to the dominance of dis- 
crete frequency components in the excitation spectrum. This 
is particularly important when the modal densities for the 
structure and receiving volume are low. In such cases the 
accuracy of the predicted resonance frequencies becomes critical 
in determining the accuracy of the predicted noise levels, 
and therefore statistical evaluations of the predicted noise 


reduction are 
dence limits 


introduced. An approach which provides confl- 
for the predictions is presented in the description 


of the analytical model. 


2.4 Study Aircraft Computations 


Following 
propeller 
fuselage , 
aircraft . 


development of the analytical model to predict 
noise transmission into the interior of an airplane 
various computations are made for the three study 


As a first step the external excitation pressure field has to 
be defined and this is accomplished on the basis of existing, 
albeit sparse, analytical and empirical information. This 
information refers to both high speed (propfan) and genera 
aviation propellers. The propfan data were obtained from a 
prediction procedure for sound pressure level provided by the 
Hamilton Standard Division of United Technologies. This 
procedure [6] permits estimation of the spatial distribution 
of the sound pressure level under free field conditions, pro- 
vided that the propfan performance characteristics are known 
or have been previously estimated. Propfan performance charac- 
teristics were estimated using Hamilton Standard procedures [7]. 
While there are several methods for predicting propfan noise, 
s ome of which predict higher levels [1], the Hamilton Standard 
procedure results from continuing investigation of propfan 
characteristics and thus it seemed appropriate to use it for 
this study. These free-field noise levels were converted to 
blocked surface pressures, as described in Section 4.6. General 
aviation propeller noise data [8,9] are used to infer trace 
velocity properties for the pressure field. It is recog- 
nized that the operating conditions for a general aviation 
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propeller are significantly different from those of a propfan 
but, in the absence of alternative data, the general aviation 
propeller measurements can be used as a preliminary model. 

When more data are forthcoming for the propfan, the analytical 
model can be brought up-to-date if need be. 

Prior to calculating the noise transmission into the passenger 
cabin, the fuselage structure is idealized into equivalent 
orthotropic panels by smearing out the assumed structural 
characteristics for the different study aircraft. In addition, 
a transmission loss spectrum is provided for the baseline side- 
wall. Absorption spectra are assumed for the cabin interior. 
These acoustic absorption spectra are based on empirical data 
for conventional interior sidewall and trim configurations. With 
this information the longitudinal spatial variation of the cabin 
sound level can be calculated, and the noise reductions deter- 
mined to reach the goal of a maximum level of 80 dB(A) . 

The next step in the computation process is to utilize the 
analytical model for a variety of add-on noise control methods 
to determine treatments which can provide the required noise 
reductions. These add-on treatments have to be such that they 
can be applied to existing conventional (baseline) fuselage 
structures without modification to the load-bearing character- 
istics of the structures. Prom the different concepts investi- 
gated, only the double wall with a heavy, limp interior trim 
panel appears to have the potential of achieving the large noise 
reductions required at the low frequencies associated with 
propfan excitation. 

In order to perform computations for the advanced noise con- 
trol concepts, where it is assumed that such concepts involve 
modifications to the fuselage primary structure, it is necessary 
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to calculate the structural response for each of the new 
structures. For each new structure, calculations are per- 
formed to determine- the characteristics of the equivalent 
orthotropic panels and these characteristics are used to com- 
pute the noise transmission. Within the scope of the analysis, 
no primary structure was identified which could provide the 
noise reductions required to meet the 80 db(A) criterion. 

2.5 Experimental Verification 

in the development of any analytical model it is always ad- 
visable to perform experimental validation of the model. 

Such validation provides confirmation that the assumptions 
introduced into the analytical model are valid for the speci c 
cases to which the model is applied. The validation also 
provides confidence in the accuracy of the analytical model. 

Several validation experiments were performed on the analytical 
model for the space shuttle payload bay acoustic environmen 
[10 11]. Since that model provided a basis for the presen 

analytical model for propeller-driven aircraft there is some 
confidence that the present model will provide reliable pre- 
dictions. Even so, consideration should be given to obtain ng 
direct verification of the propeller noise transmission mo e 
Model verification did not form a part of the study reported 
herein, but the study did include an assessment of the im- 
portant parameters in the model and recommendations for 
prospective verification experiments. 

The choice of experiment depends to some extent on the actual 
objective of the tests. If the objective is to confirm the 
specific noise reductions predicted for the three study 
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airplanes, then the experiments should be designed to model 
the important parameters of those study aircraft. Alterna- 
tively, if the objective is to validate the analytical model 
under more arbitrary conditions, an existing propeller-driven 
airplane could be selected, the analytical model used to pre- 
dict the noise transmission into that aircraft and an experi- 
ment devised to confirm those predictions. It should be borne 
in mind, however, that if the test conditions differ radically 
from those corresponding to the study aircraft, the validation 
of the analytical model may be performed for a situation un- 
representative of propfan operations. 

The evaluation of potential test procedures discussed in 
Section 9 includes both of the above objectives although the 
emphasis is placed on the validation of the analytical model 
as applied to the study aircraft. This emphasis is chosen 
because it allows the analytical model to be tested under 
conditions appropriate to the subject under consideration — 
high-speed, propeller-driven aircraft. 


3. ANALYTICAL MODEL 

3.1 Power Balance Approach 

The basic concept of the analytical model is that of the balance 
of acoustic power flow from the exterior of a bounding surface to 
an interior-contained volume. This approach has recently been 
developed by Pope and Wilby [4] for application in the acoustic 
modeling of the payload bay of the Space Shuttle orbiter vehicle. 

The first step in the analysis considers the vibrational response 
of the fuselage structure to the external pressure field. The 
model takes into account the external pressure field in terms of 
level distribution over the fuselage surface, spectrum shape, and 
spatial correlation. It describes the structure. in terms of the 
mode shapes of equivalent orthotropic panels, including frequency- 
dependent effects of skin and stringer and. frame stiffnesses 
and masses, pressurization, and curvature. The coupling between 
the external pressure field and the fuselage vibrational response 
is formulated in terms of modal joint acceptances. 

The second step considers the power radiated by the fuselage 
structure into the cabin. Acoustic power is radiated into the 
volume by the structural modes, and the power is accepted by the 
acoustic modes of the volume. Resonant and nonresonant response 
of the structure and volume are considered in the analysis, resonant 
response being that for which the participating modes have resonance 
frequencies within the frequency band of interest. Nonresonant 
structural response can be either mass— controlled (resonance fre- 
quencies below the band of interest) or stiffness— controlled 
(resonance frequencies above the band). However, for present pur- 
poses, the stiffness-controlled response is excluded, as. the 
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radiation efficiencies of stiffness-controlled modes are low 
compared with those of resonant and mass-controlled modes (Eq. (6l)). 

The final step determines the space-averaged sound pressure levels 
in the interior volume by equating the net inflow of acoustic 
power to that which is absorbed within the interior space. The 
model accounts for any transmission of acoustic power from the 
interior back to the exterior of the fuselage. 

The following sections present a detailed derivation of these acoustic 
power flow equations, with particular emphasis placed on the 
following developments which have specific application to propeller 
noise: 

1. Transmission of discrete tones through a structural element 
into a cavity; 

2. Estimates for the mean and standard deviation of the acoustic 
power flow; and 

3. Closed-form solutions for the coupling between a nonhomogeneous 
excitation and the fuselage vibration response. 

3 . 2 General Formulation 

In the general formulation, the transmitting structural system is 
assumed to cover a cavity which has absorbing walls as the interior 
surfaces. The frequency range considered is sufficiently wide to 
assure that the cavity possesses at least one mode which is resonant 
in the frequency band of interest. The transmitting structure is 
assumed to have modes resonant below and within (or simply within) 
the band of interest. 
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The transmitting structural system will respond to excitation by 
an external random pressure field according to the relation 

W|p(x,w) = j G( x,x' ;w) [P^( x T ,w) - P^(x' ,w)]dx' . (D 

G(x x' ;«) is the Green's function for the structure^ P T ( x *■ ) 
respresents the external exciting pressure, and P T (x ) the induced 
interior pressure on the structure. Mathematically W T (x,w), 
P°(x'co), p^x'jW) are Fourier transforms of truncated records (in 
time) of length T and are related to the spectral components of 
the random displacements and pressures. 

The induced Interior field Is related to the displacement response 
of the structure by 

P*(x',w) = -pjW 2 f G p (x,x’ ;w)W T (x,a))dx (2) 

where G p (x,x’ ;w) is the Green's function for the interior space. 


The external driving field is given by 


P^ x ’ 




^ = 




p to 
e 


f G°( x ,x ' ;w)W T (x,w)dx (3) 


where P (x’ ,w) is the transform of the blocked pressure, 

G°(x,x' ju ) 1 ) 1 is the exterior field Green's function (source point 
on the structure), and the integral represents the radiated pressure 
field. The ambient densities in the interior and exterior spaces are 
given by p ± and p e , respectively, and the associated speeds of 

sound are c^ and c g . 
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3.2.1 Structural Response 


These coupled Integral equations can be reduced to a system of 
simultaneous algebraic equations to obtain the modal response £ t 
of the transmitting structure [ 13 : 


[- Y M + co 2 ( p.T^ 
r r i 


P e J rr >]E r T +' 


co 


I Cp;i" ♦ 

s^r 




(4) 

In Eq. (4) external and internal acoustical coupling of structual 
modes is described by the intermodal coupling coefficients defined 
as [13 


J rs (oo) = JJg°(x,x’ ; oo)^^J ^ (x , )^ S (x)dx' dx, ' 
I rs (co) = // G p(^ J ^’ ; w )'4 r (x’ )t S (x)dx' dx. , 


(5) 



(o>) 


is the generalized blocked force 


r pb«, T = / p bj T <x '’“ ),| ' r(x ' >ax '’ (6) 

ip r (x) is the eigenfunction of mode r, is the modal mass, and 
is a receptance given by 


Y r = w 2 [l - (o)/to r ) 2 - in r ]. 


(7) 


where n r and co^ are the modal loss factor and resonance frequency, 
respectively. 
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The matrix form of Eq. (4) is 


^rs-^rT* 


= {-r r } 

1- 1 0 5 > 


pbJl 


T 


and has the solution 


«rT> “ C “rt^- r pbt T } > 


from which the displacement response of the structure is obtained: 


W T (x s u)) = 


r 3 t 


' 1,I ’ <x>a rt r pb«, T ( “ ) 


( 8 ) 


The one-sided cross-power spectral density for the displacement at 
points x,x' is given by the relation 

S (x,x';a>) = lim( 2/T) ( l/2ir)W rT1 ( x)W^( x ' ) , (9) 

T ^CO 1 1 

and, from Eq. (8) 


As “rt^&n 


I vV ,i), 1 <i ' 1 I/ s p (i " i?i ;w) 

r ,t , l ,n 


x i|> t (x ,, )i|> n (x ,M ) dx" dx"' , 


( 10 ) 


where S (x",x m ;w) is the one-sided cross-power spectral density 
of the exterior blocked pressure, i.e. 

S (x" ,x' n ;a>) = lim (2/T)(l/27r)P . (x”,co)P* (x’" ,a>) . (11) 

P t-x» b T DX T 
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3.2.2 Power Radiated into the Cavity 

The one-sided cross-power spectral density between pressure at 
position x and velocity at x' is 

Sp . v(x,x* ;w) = lim ( 2/T)(iu)/2Tr ) P~( x,w )w!( x’ ,(o) . (12) 

1 t-*» 1 1 

Equation (2) yields 

Sp. v(x,x' ;w) = -ip. co 3 f G ( x,x" ;w) S ( x" ,x f ;w)dx" . 

1 1 J p w 

Let x'-> x. Then, the spectral density of power radiated is ob- 
tained by integrating Sp^vCxjX ;w) over the transmitting area; i.e. 

^r ad ( w ) = / s Pi v ^»^ »w)dx. (13) 

This quantity is complex, and the real part is the real power. 

Prom above, 

wj^(u) = - iPj^ 3 / / ®p(^,^’ ’ w )^w(^' ,^> w )dx’ dx. (14) 
Using Eq. (10), this becomes 

W rad ( “ > = ~ lp i“ S J “rt“*n £ f_ V *>*' ;a>)* r (J) 

r,t,£,n x x* 

x ip^(x' )dxdx’ J J S p (x",x' M ;co)^(x") 

? x ,M 

x ^ n (x ,M )dx ,, dx ,M . (15) 
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X* 2/ 

Now introduce the intermodal coupling coefficients . I (-to) to 
obtain 


= - ip^ 3 J “rtV lI ’ lt( “ ) / / S p (x,i';«.)/(x)t n (x’)<a X dx'. (16) 

r,t s &,n 


If the cross terms for the external and internal acoustic coupling 
coefficients for the structural modes can be neglected, Eq. (16) 
reduces to: 


■ 


ip, uf 


I 


T rr, v 
I (co) 


; J J"S p (x,x’ ;w)\Jj r (x)ip r (x , )dxdx* , 

(17) 


r>r> 

MY - w 2 (p J + 
r r e 


Pji 


rr 


)| 


The assumption of negligible coupling of the structural modes by 
the radiated sound field (i.e., I rS = 0 for r / s) is equivalent 
to the assumption of weak coupling between the structure and the 
cavity. Justification for this assumption is given in [ 4 ] , where 
it is pointed out that the assumption does not preclude well- 
coupled modes that occur when acoustic and structural modes have 
resonance frequencies closely spaced relative to the modal band- 
widths. The assumption that the cross terms for the external 
coupling coefficient can be neglected (i.e., J =0 for r r s) 
is adopted in [4] for a homogeneous excitation field. It is 
assumed now that the cross terms will also be negligible for the 
present inhomogeneous pressure field. 


The cavity Green’s function is 


G p (x,x' ;w) 


I 


( x) c[> n ( x * ) 

(k£ - k 2 )///? 2 (x)dV 


(18) 
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where 4> n (x) is the complex mode shape, and the complex eigenvalue 
is 


k - = k - iK , 
n n n 


(19) 


The acoustic surface admittance of the cavity walls is 


S = K - iff = P ± ° ± / z 

where % is the normalized conductance and a is the normalized 
susceptance of the wall. z is the inner surface’s specific acoustic 
impedance, that is, the ratio of pressure to normal velocity at the 
wall, 

z = p/v. 

Analysis is confined to the case |3| <<1. Then 


which is real. With the normalization such that 


f<t > 2 dV = Y/e , 

J Y n n 3 

where V is the volume of the cavity and e n is mode dependent, the 
joint intermodal coefficient in Eq. (18) becomes 


I rr (w) 



V(! 2 - k 2 ) 
v n ' 


ip r (x)ilj r (x' )dx dx’, 


which reduces to the form 


T rr, v A 2 
I (0)) = y 


n 


(k 2 

n K rx 


k 2 ) 


f 2 (n,r) , 


( 20 ) 


where the function f(n,r) is the coupling factor between the 
structure and the cavity. 

. f(n,r) = j f <i> n (x)i|j r (x)dx. 


The spectrum of real power flowing inwardly can now be written. 



2p ± (o 3 A 2 y* e n k n K n V f 2 (n,r) 

= L 1^ - * 2 F L |M Y - V(p J rr + P.I^I 2 

n r 1 r r e 1 


X 


//• 


C (x,x' 
P ’ 


;ix>)ip r (x)^ r (x' )dx dx’ 


( 21 ) 


where C (x,x';w) is the real part of S p (x,x' ;<*)); i.e. the co-spectrum 
or co-spectral density function of the exterior blocked pressure, 
and 


k 2 - k‘ 
n 


2 = (k 2 - 
n 


n 


- k 2 ) 


K k o' 


( 22 ) 


Now k n = u n ^°± and K n = w n n n^ 2c i where n n is the loss factor 
of the acoustic mode with resonance frequency a> n . Thus, 


and 


Vn = %% /Za l 



( 23 ) 
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The band-limited input power. 


for a bandwidth Aw, is 




w Jad ( “ )d “ 



Then, under the assumption that the radiation loss factor is 
relatively small [4], the term in I can be neglected, and 


MY - w 2 (p J 


rr 


r r 


h 1 


rr 


)| 2 « 


M 2 |Y 


M 2 w 4 
r r 




where n r is now defined as the sum of the structure's dissipative 
and external radiation loss factors. 


The band-limited power flow into the volume becomes 


Re 



p ± c 2 w 3 A 2 ^ 

V J 

Aw 

Aw 


V' £ n 1 \i 

^ w 2 [(.l - ^) 2 + n 2 ] 
n n v w 2 n 
n 


f ^njr) JJ Cp(x,x' ;w)'4 r (x)’l ,r (x’ )dxdx* 


M 2 w 4 
r r 


[(■ - Sr)’ * <] 


/ dw . (24) 


It is clear then that, in general, calculation of the band-limited 
power flow into the volume requires estimates of both acoustic and 
structural resonance frequencies and mode shapes and their respec- 
tive loss factors, together with a model for the excitation co- 
spectral density function. 
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3.2.3 Power Absorbed by Cavity Walls 

When calculating the acoustic power absorbed by the interior 
surfaces of the fuselage cabin, it will be assumed that the 
surfaces are locally reacting. Then, the motion at one portion 
of the surface is dependent only on the acoustic pressure incident 
on that portion. This assumption is valid for the limp, absorptive 
surfaces envisaged for the cabin sidewall treatments. 

From [4], the acoustic power absorbed by the walls is given by 
the general relationship 


TT int _ 1 

¥ . = 

abs Pi c i 


2 2 

n w 



A 5 
w w 



> 


where A is the absorbing area of wall w, ? w ls the wal1 conductance, 
<p 2 > is the space-averaged, mean square pressure in the volume V 
due to the response of the nth acoustic mode, and e is a modal 
factor. On an acoustic mode-by-mode basis, this becomes [4] 


W 


int 

abs 


V 


p i c i 


<p 2 > 

*n 


^ * e x- 
a 


+ e 


V + V 


q 


b 


+ e 


? z+ + ^ 


(25) 


n=p,q,r 


where it is assumed that the volume is approximately rectangular 
in shape, with dimensions (a,b,c). The wall conductances are now 
denoted by C x+> £ x _, etc; and e p , e q , e r are modal participation 
factors with values of 1 or 2 depending on mode order. 


Equation (25) requires a knowledge of the wall conductances before 
the absorbed power can be calculated. Such information 
is not readily available in practice, although values of the 


22 - 


absorption coefficients are often known or can be estimated. 
Therefore, it is necessary to convert Eq. (25) into a form 
that incorporates the absorption coefficient a rather than the 
conductance ?. Following the approach used in [4], the conver- 
sion has been accomplished by modeling the low-frequency absorp- 
tion characteristics of the volume in such a manner that the 
relationship maps into the correct form at high frequencies. The 
rate at which the formulation approaches the high frequency value 
is comparable to the accumulation of modes in the volume. The 
resulting equation is 




<rf> 


abs - 5^ + VV r e + (ol y+V + %-\-K + ( V A z+ + V A z- 




X]. 


(26) 


where r x = 

e N n 


n 


^ e q e r 
n H n n 


N n is the number of acoustic modes, and 


the assumption is made that the exterior input and interior coupling 
are complex enough to assure that the modal responses differ little 
from the mean, so that <p^ 2 > = <p^>. "*' n h^-Sh-fnequency regime, 

where there are many o.blique modes, the average conductance is about 
one-eighth of the average absorption coefficient [ 13 . 


3. 3 Representations for the External Pressure Field 

A joint acceptance function j 2 (uj) can be introduced to quantify the 
the coupling between the external pressure field and the rth 
structural mode. The function can be defined as 





S D (x,x’,a))^ r (x)^ r (x' )dx dx' 


(27) 
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where S (u>) is the one-sided power spectral density of the 

exterio/blocked pressure field at a reference location X Q . As 
before, s'(x,x' ;u>) is the cross-power spectral density of the 
exterior blocked pressure, defined by Eq. (11) or by 

00 

1 f , 1 V -ia)T 

S p (x,x ? ;to) = — / Rp( x J x > T ) e dT 

—00 


where R p (x,x’;x) is 
pressure over the exterior 


the average cross-correlation of the blocked 


/ , \ lim 1 I 

Rp ( x > x > T ) t->°° T J 


T/2 


P 


b& 


(x’,f) P b „Cx’,t + T)dt 


-T/2 


3.3.1 Homogeneous Exterior Pressure Field 

For a homogeneous pressure field, such as is normally observed 
beneath a turbulent boundary layer or in the far field of a jet 
(provided the structural element considered Is not too large), the 
following form of S p (x,x' ;«) is found to represent the measured 

data quite accurately. 


S p (x,x T ;w) = 


S p U;u>) = 


S ((o)e' 


■ck 5| Q ik? 


(28) 


where X = x ’~ x > c is the correlation decay parameter, 
excitation trace wavenumber, and reference location Xo 
because of the assumption of homogeneity. 


k is the 
Is arbitrary 


The co-spectrum is 

C (T;a>) = S_(w) ■e" ck ^Uos kf . 
p ' P 


(29) 
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Both the cross-spectral density function and the co-spectrum can 
be assumed to be separable In directions parallel and perpendicular 
to the direction of propagation. 

3.3.2 Non-homogeneous Exterior Pressure Fields 

White p.4] has shown that estimates of the joint acceptance can be 
made for a pressure field whose power spectral density varies 
smoothly over the structural mode shape. The assumptions that the 
co-spectrum Is homogeneous and that the spectral density Is 
moderately uniform over the surface were made. 

In general, strong variations in excitation co-spectrum can occur 
over the surface of the transmitting structure. Consider the 
situation where a periodic point source is located at P, adjacent 
to the (x,y) plane as shown in Figure 1.- An expression is re- 
quired for the cross-power spectral density of the pressure field 
over the two-dimensional (x,y) plane. (]Xb,Yo) is located at a 
distance d from P, such that the line joining (Xo,Yo)to P is 
perpendicular to the (x,y) plane. 

The pressure at points away from P varies inversely with r, the 
radial distance from P; i.e. 

p(x,t + t*) = p(X 0 ,t) 

x 

where t* = (r x - d)/c = £(|x - X 0 | 2 + A 7 ^ 1 ^ 2 - dj/c , 

X * 

T o = (r x» “ d)/c > and 
|x - Xo| 2 = (x - X 0 ) 2 + (y - Y 0 ) 2 . 
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FIGURE 1. 


PERIODIC SOURCE LOCATED ADJACENT 
TO (x , y ) PLANE 




The cross-correlation of the pressures at x and x* can be written 

J72 

R (x,x T ; t) = lim ^ / p(x,t)p(x',t + x)dt 

p m<» 1 J m _ 

T/2 

where p(x,t) = ^ Q " | % + ^)h p( *Q> b " T S } 
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If p(Xo,t) is periodic with period T q and frequency w Q = 2ir/T o , 
the above equations become 
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where q(x,x’) = d 2 /[(|x - Xo | 2 ' + d 2 )(|x’ - X6 1 2 + d 2 )]^ . 


With 


p ( Xo , z ) = ^ 


ibw 0 z 

c e , we find 

m ’ 
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1 f 2 -ibco z 
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Then, on substitution into Eq . (30), 
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R (x,x* ; t) 
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q(x,x’ ) 
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The cross-power spectral density follows, on transforming and 


X X 

substituting for and t q , as 


S (x,x';u>) = 2 q( fe" £ 8 

P ° b=- 
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p„ o 
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where 6(w) 


CO 



The co-spectrum Is 

00 

C (x,x';o>) = q(x,x’)^ cos 

b=l 

x 6(w - bto ) (33) 

o 

where <p (bw )> Is the mean square pressure at Xo attributable 
Xo o 

to the bth harmonic. 

The co-spectrum (Eq. (33)) I s a complicated function of observer 
positions relative to source position but depends on the ampli- 
tude of the pressure field relative to that at refererence position Xo 
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The function q(x,x' ) expresses the spatial decay of the power 
spectral density between the observer positions x and x* and the 
reference point Xo. 

Considerable effort has been expended in the literature to define 
the function q(x,x') and coherence and phase relationships (cor- 
relation functions) for such excitations as jet noise [15,16]. In 
general, an approximation approach has been followed using the 
homogeneous cross power spectral density function of Sec. 3*3.1, where 
S (w) has been obtained by averaging the measured data over the 
structural area of concern and where average values of k and c 
have been selected. This is valid where the spatial variation 
in power spectral density is small, but it is likely to be in- 
creasingly inaccurate as the spatial variation becomes more severe. 

For example, contrast a point excitation with a homogeneous excita- 
tion. The inclusion of the spatial distribution function q(x,x* ) 
permits a more realistic and precise modeling of the coupling 
between the nonhomogeneous excitation and the fuselage response 
to be performed. 

On the other hand, few measurements have been made of the cross 
power spectral density function of tonal nonhomogeneous pressure 
fields. From the available data [8,9], it can be argued that the 
following is an appropriate general formulation of the cross power 
spectral density for use in near- field and far- field applications. 

The selected function is 


. s p(x “ Xo, x* - Xo; bw Q ) = q 1 (x - Xo)q 1 (x’ - Xo)R p (?;b« 0 > 

where q^(x - Xo) represents the spatial variation in tonal rms 
pressure over the structure surface; R p (5;bu) Q ) is a homogeneous 
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correlation function for the excitation [see Section 3.3.1 ] and 
f _ x'- x; defined by 

-ck( bw ) | £ | _ 

R (f;b.u> ) = e cos[k(b« )?]; 

p 5 o u 

and k(bw ) is the excitation trace wavenumber at frequency bw Q , 
averaged over some appropriate spatial extent. 

Often the excitation coherence will be close to unity (i.e. c = 0) , 
as observed in the propagation direction in model tests for plane 
waves propagating over a structure Qj] and on the Aero Commander 
fuselage [8 , 9] measuring circumferentially in the plane of rotation 

of the propeller. 

At present, values of c and k(bw Q ), and the nature of q^x - Xo), 
have to be estimated empirically. Ray acoustics, recent Hamilton 
Standard test data, and Hamilton Standard prediction methodology 
can be used to this end. For example, Hamilton Standard theoreti- 
cal data [6] can be used to estimate the function q ] _(x - Xo): it 

is found that an exponential function describes the predicted 
spatial variation in power spectral density with fair accuracy 
in both the axial and circumferential directions. Then, considering 

for example 


a joint acceptance function for the nonhomogeneous field can be 
defined at frequency bw Q in the form 

j ^ (b.to ) = ~~2 J cos[k(bw o )(x' - I) ]/(I)/(x f )dxdx’ . 

X X* (34) 
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This expression is separable in the appropriate coordinate direc- 
tions and may be evaluated for various mode shape functions. 
Furthermore, Xo may coincide with the edge of a structural member 
or it may be located at an arbitrary position on the structural 
surface. A general derivation is presented in Appendix A. The 
form of joint acceptance and the dependence on coherence, 
amplitude decay rate, and wavenumber matching between excitation 
and response is discussed in Section 3.6. 

We note that, from say Eq. (27) 

oo 

J f C p (x,x^w)/(30i(i r (?)dxdx» = ^ A 2 j^( w ) <p2 o (ba) o )> 6(w - bw o ). (35) 

x x* b=i 

3 . 4 Tonal Power Flow Equation 

The preceding analysis of the power flow into a cavity has placed 
no restrictions on whether the structural or cavity response is 
resonant or nonresonant. Application of the analysis to typical 
baseline aircraft (to be discussed later in Section 6) indicates 
that the modal density of the acoustic cavity formed by the air- 
plane cabin is sufficiently high that the power inflow is dominated 
by resonant response of the cavity. Thus, Eq. (24) for the power 
inflow can be simplified, from a computational standpoint, by con- 
sidering only resonant response of the cavity. This is achieved 
in Eq. (24) by performing the summation for cavity modes of order n 
only over those modes that lie within the bandwith Aw of interest, 
as denoted by the summation limit neAw. Then Eq. (24) becomes, 
with substitition of Eq. (35), 
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where = bW Q , w fe ~ w n since both lie in Aw, and 

K (u L 5 Re K* (u) L- 

The summation over excitation harmonics is also limited to those 
harmonics lying in Aw (i.e. be Aw). 

The bandwidth Aw can be chosen such that only one harmonic of the 
excitation lies within Aw. Then Eq. (36) becomes 
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where w^eAco. 

3. 5 Power Flow Statistics 

The power flow equation presented in Eq. (37) is based on a mode- 
by-mode calculation of the response of the structure and cavity to 
an excitation at a given frequency. This is feasible from a compu- 
tional point of view where there are relatively small numbers of 
structural and acoustic modes in the frequency band of interest, 
but the computational requirements become excessive when the number 
of modes is large* Thus, for the present study the power flow 
equations are considered for three frequency regimes. These regimes 
are defined in terms of modal overlap, which is essentially the 
ratio of modal bandwidth to modal separation for either the structure 
or the cavity as the case may be. 

At low acoustic modal density and loss factor, such that the modal 
overlap is much less than unity, the cavity response depends 
intimately on the cavity resonance frequencies and loss factors 
and on the excitation frequency. The exact result given by Eq. (37) 
should be used for calculation of the power radiated into the volume; 
approximations for the modal admittance should not be made. 
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The high-frequency case is one in which the modal overlap is 
high for both the structure and the cavity. However, there can be 
intermediate frequency regimes where the modal overlap is high 
for either the structure or the cavity but not for both. Consider, 
for example, the case with a high acoustic modal overlap. This 
situation is one in which a tone of frequency drives a group 
of acoustic cavity modes contained in a narrow frequency band Aw. 
Then, with a modal overlap of approximately unity, it becomes 
appropriate, for computational efficiency and because of a lack of 
precision in mode shape and resonance frequency predictions, to 
seek solutions to the power flow analysis which are statistical 
in nature. Of particular interest is the expected value of the 
power radiated into the cavity E[W^ n (w )] and the corresponding 
variance o 2 ] . 

3.5.1 Expected Value of Power Flow into Cavity with High Modal 
Overlap 


The power radiated into the cavity depends on the acoustic modal 
admittances as well as the interior coupling parameters f 2 (n,r), 
mode normalization constants e n , and loss factors ri n . On assuming 
an average acoustic loss factor n and an average product of f 2 (n,r) 
and e , <e f 2 (n,r)>, the expected value of the internally-radiated 
acoustic power may be written as 


E[W in (^)] 



n <e f 2 (n,r)> 
n n ’ 


x E[GU b ,w n ) ] (38) 



MH 

where G (. w b , w n ) = ^ 

neAw 
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Thus, as the modal overlap Increases, one seeks an estimate for the 
modal admittance which Is Independent of the precise modal details, 
so that the Internal coupling becomes a function of modal overlap 
(and hence modal density and average loss factor) only. In the 
case of high acoustic modal overlap, there Is Interest In the 
expected value of the above modal admittance function GCu^ ,u> n ) for 
fixed values of o) n (neAw) as w fe samples In the band Aw. This Is 
equivalent to the expected value of G for fixed w^ , where the 
a) ’ s are distributed in frequency according to their distribution 
functions for an ensemble of (slightly) different cavities. 


Now if 






g( 0), , w ) 
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( 39 ) 


where N . is the number of acoustic modes resonant in Aw. Then, 
neAw 

assuming that w b and w^ are statistically independent, that the 
probability density of w^ in Aw is (Aw) -1 , and that the probability 
density of w n in Aw can be derived from the acoustic modal density/ 
it can be shown that (see Appendix B) 
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This expected value for the acoustic admittance of a group of 
resonant acoustic modes can be substituted into Eq . (38) to give 


E[W (Wb)] = 
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where <e f 2 (n,r)> = — Y e n f 2 (n,r). 

n w neAw 4- 


Thus, the expected value of the band-limited power flow into a 
group of acoustic modes from a tonal excitation depends only on 
the acoustic modal density of the acoustic modes and is independent 
-of the acoustic loss factor. 


It may be shown, see Appendix C , that 

w?V rev 

W <£ n f2(n ’ r)> = ^ (a) b } 

where i 2 (aO is the joint acceptance for the structure exposed 
to a reverberant pressure field. On substitution, Eq. (4l) becomes 


- 36 - 


(42) 


E[W ln (u, b )] 


Aw 


pj/4 A 

2ttc . 


< ^Xo (w b ) ' ^ 




rev 
2 '“ ) 


b^K 


]Yh 


r 


Y (w, ) 
r b 


Additional results, for the power flow from resonant structural 
modes (i.e. modes whose resonant frequencies w^ lie with.in the 
analysis bandwidth Aw) and the power flow from mass-controlled 
modes (w < Aw), may be easily derived. These are presented in 
Section 3.5.3. 

3.5.2 Variance of Power Flow into Cavity with High Modal Overlap 


Calculation of the variance for the internally-radiated acoustic 
power (with high modal overlap in the cavity) requires calculation 
of the average modal acoustic admittance function G(w^,w ), where 
by definition [C.l] 


a 2 [G(w 
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= E[G 2 (w b ,w n )] - E 2 [G(w b ,w n )] 
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and E[G(w b ,w )] has been already derived. 


Now 


E[G 2 (w b ,w n )] = E 
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as w b and w n are statistically independent. 
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Then, from Eqs. (40), (43), and (44) 
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Thus, the variance is inversely proportional to the cube of the 
band-averaged acoustic loss factor and directly proportional to 
the acoustic modal density. 


Of interest is the normalized standard error e r or the coefficient 


of variation, defined as 


a[W. n (w b )] 


'r E[W ln (w b ) J * 


(46) 


Prom Eqs. (38), (40), and (45a), 


e = [ 2ttti n (w H )w H ]~' g 
r n n b b 


where n n (w, )w, is the modal overlap term for the cabin acoustic field, 
n n b b 


Thus, the standard error decreases as acoustic modal density and 
loss factor increase. 
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3.5.3 High Frequency Formulation 

When the modal overlap Is high for both acoustic and structural 
modes, a high-frequency solution can be found for the acoustic 
power Inflow. 

S.5.3.1 High-Frequency Acoustic Power Flow via Resonant Structural 
Modes 

Adapting Eq. (42), the power Inflow from resonant structural modes 
is given by ' 
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E[w ln (, Wb )] = < pj 0 (» b ) > e J. r / <\ 
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With j ^( oj^) and j^(w. H ) set equal to band-averaged values <j*(w H )> 
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and <jj(« ^ respectively, where the subscript P Is used to 

denote that the joint acceptance Is based on a progressive excita- 
tion field. 
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Using an approach similar to that used in Section 3.5.1, 
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where n is the band-averaged structural loss factor and ^(u^) 
is the structural modal density. 


Hence , 
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It may be shown [4] that the band-averaged radiation resistance 


<R^ n ^> can be calculated as 
rad 


2 n 2 


p . to A 
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which leads, on substitution In Eq. (49), to the alternative 
representation for E[W ln (w b )] as 


E r w in ( “b )] ' 5^ ^loK’* W 


<J>P> P < R rad > (51) 

M ^r 


The normalized standard error is found to be the sum of the 
normalized standard errors Inherent In the structural modeling 
and in the coupling to the cabin acoustic field, viz. 

e r = [2iTn r n r (a) b )w b ]- 1 2 + [2-rrn n n n (m b )a) b ] 

where n r n r (u> b )w b is the modal overlap term for the structure. 

3.5. 3. 2 High-Frequency , Nonresonant Acoustic Power Flow 

It is assumed here that, as there is a large number of modes with 
resonance frequencies below the band of interest, only these modes 
will be considered for nonresonant power flow calculation. This 
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restriction is made because, in general, these mass-controlled 
modes are more efficient acoustic radiators than are stiffness 
controlled modes. 


If it is further assumed that for these mass-controlled modes 
w b >> w r* then from Eq. (7) 
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Substituting in Eq. (42) 
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The net acoustic power flowing from outside to inside from non- 
resonant, mass-controlled modes is then 
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where the second term represents the acoustic power flow from the 
cavity to the exterior via nonresonant structural modes. 


3.5.4 Net Power Flow 

The high-frequency approximations developed in Section 3.5.2 and 
3. 5. 3.1 provide estimates of the power flowing from the outside to 
the inside of the cavity. However, in order to obtain an estimate 
of the net inflow of acoustic power, the outward power flow from 
the cavity to the exterior has also to be taken into account. 

It may be shown [4] that the analogous result for the acoustic 
power re-radiated from the fuselage due to excitation from the 
acoustic field inside the cavity is given by 
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where <P?(w b )> is the induced interior space-averaged mean square 
pressure attributable to the bth harmonic. 


The net acoustic power flowing from outside to inside from resonant 
structural modes is then obtained by combining Eqs. (51) and (53) 
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or, similarly, by combining Eqs. (42) and (53). 


3.5.5 Summary of Power Balance Equations 


The preceding analysis has developed the acoustic power flow 
equations for a simple cavity in a rather general manner, although 
emphasis has been placed on high-frequency approximations that are 
appropriate to the baseline airplanes considered in this study. 

In addition, an excitation field with discrete frequency components 
and spatially-decaying amplitude has been included as a representa- 
tion of a propeller noise field. 

Power flow equations have been developed for three frequency 
regimes — a low-frequency regime where the modal overlap is low 
for both the structure and the cavity , a high-frequency regime 
where the modal overlap is high for structure and cavity, and an 
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intermediate regime where the modal overlap is low for the 
structure and high for the cavity. In the low-frequency regime, 
the power flow equation, Eq. (37), represents the net inflow 
into the cavity. However, when approximations are introduced 
for the high acoustic modal overlap case, power outflow has to 
be considered separately, Eq. (53), because Eqs. (42), (51), and 
(52) give estimates of the gross rather than net power inflow. 

As shown in Section 2, the net power inflow has to be balanced 
by the power absorbed in the cavity. Equation (26) has been 
developed to provide an estimate of the power absorbed by the 
surfaces of the cavity (or cabin). In addition, acoustic power 
Is absorbed by non-bounding surfaces such as seats, passengers, 
and service structures. The power absorbed by such items is 

<P ^ V 

W abs ’4" P 1 cY £ a j A j 

j 

where ou is the statistical absorption coefficient associated with 
area and <p*> is the mean square pressure in the cabin interior. 

In a more general situation, the cabin volume might be composed 
of a series of subvolumes, either coupled directly or partitioned 
by barriers. Such a situation was followed in the analysis so 
that noise control methods such as partitioning of the cabin 
volume could be considered. 

Equations (26), (37), (42), and (51) through (54) constitute the 
basic analytical model. However, these must be supplemented with 
information regarding cabin shape, excitation field, interior 
absorption, and sidewall configuration. The joint acceptances, 
interior coupling factors, and structure and; volume resonance 
frequencies; allowed' wavenumbers , structural loss factors, radia- 
tion resistances, and modal densities must be estimated for the 


airplane configurations under study. The following sections 
address these items. 

3.6 Joint Acceptances 

3.6.1 General 

The joint acceptance function describes the coupling between the 
excitation field and the structure, and in the present analysis 
it is defined, for a single harmonic at frequency %, by Eq. (27) 


J>b> = 
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2 <p Xo ( ' a) t 


7 > 


/-/- 

Jyyyi' 


C (x,x' H r (x)i4 r (x’ )dxdx’ 
p b 


where A is the panel area, C p (x,x';w b ) and < pj 0 ( w b ) > 5 respectively 
are the co-spectrum and mean square pressure at Xo for the bth 


harmonic of the nonhomogeneous excitation field, and ip (x) is the 
shape of the rth structural mode. 


In order to evaluate the joint acceptance, it is necessary to 
have representations for the excitation field and the structural 
mode shapes. For the present model, two excitation fields are 
encountered, one being a nonhomogeneous propagating acoustic field 
generated by the propeller and the other a diffuse acoustic field 
associated with the interior acoustic field of the cabin. The 
fuselage structure is composed of curved panels whose boundary 
conditions will probably lie between fixed and simply supported. 

Mode shapes are complicated by variations in skin thickness and 
the presence of stiffeners, windows, etc. As it is not possible 
to account for the detailed forms of the mode shapes, the simpli- 
fying assumption is made that the structural modes, can be represented 
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by sine functions such as are associated with simply-supported 
boundary conditions and rectangular plates. Experimental and 
analytical data for stiffened cylinders [18] and for the Space 
Shuttle [4] support this approach. In any case, the precise mode 
shape used in the representation is probably not too critical. 

For example, joint acceptances have been calculated for rectangular 
panels using mode shapes associated with simply-supported and 
clamped boundary conditions. The excitation consisted of propa- 
gating waves in one case 019] and a subsonic turbulent boundary 
layer in another [20]. At coincidence, the joint acceptance 
was 0.5 to 2.5 dB higher for the simply-supported panel than for 
the clamped panel. Off-coincidence differences between the 
results were less than ±2.5 dB. More significant is the representa- 
tion for the excitation field. 

3.6.2 Nonhomogeneous Propagating Acoustic Field 

As discussed in Section 3* 3- 2 and in Appendix A, the strong spatial 
decay in mean square pressure away from the location of peak 
intensity is represented in this analysis by an exponential ampli- 
tude decay function in the expression for C p ( x,x’ ;to b ) so that 


C ( x , x ’ ; to, ) = 
p * * b 


e -a|x-Xo| e -a|x» 


-Xo| -c | x ' -x | ri , — . .. , . . 

'e cos [k( x - x) ]<pj o (oo b )> 


(55) 

where a is a pressure amplitude decay rate parameter, 
c is a correlation decay factor, and 
k is the wavenumber at frequency to^ . 

It is assumed that C (x,x' ;a) b ) is separable in the longitudinal 
and transverse directions; i.e. 
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where 


C p (x,x' ;u b ) = C p X (x5X ' ; V C py (y ’ y ' ;W b } <P Xo (w b }> 


-a x |x-xd - a x b'-x°! e -°x k * |el M w , 
C (x,x T ; w ) = e e e x 

px b 


C py (y ’ y ' ; “t' ) “ 6 


a y|y- Y °l - a y|y'- Yo l r c y k y |n| 


y 


cos k n , 

V 


K = x’ - x, 


n = y' - y> and 

Xo = ( Xo , Yo ) . 


When a and a„ are zero, the pressure field corresponds to that of 
fully-developed turbulent flow (c x and c y both greater than zero). 
For such a case, the effects of variations in (c x ,c y ) on the joint 
acceptance have been Investigated in detail elsewhere [17]. Of 
concern here is the dependence of the Joint acceptance on varia- 
tions in (a , a , ) , in particular for a coherent excitation, i.e. 

(c c ) cloL to zero. Some simple calculations have been performed 

to demonstrate the influence of ( a x > a y)* 

Figure 2 illustrates the effects of variations in a x on the one- 
dimensional joint acceptance for the fourth mode of a 10m long 
beam when c = 0. The ratio of excitation to structural wave- 
number (k /k ) varies from 10“ 1 (structure wavelength X m << excita- 
tion wavelength X x ) to 10 2 (excitation wavelength << structure 
wavelength). Maximum and minimum values of j* m occur for a x and 
c both zero, corresponding to plane waves propagating over the 
beam length with no amplitude decay. However, as a x takes 
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increasingly higher values, the joint acceptance depends less 
strongly on the matching between structural and acoustic wave- 
lengths: the strong oscillations in j 2 , as k x / k m varies, are damped 

by the presence of significant spatial amplitude decay. It is 
clear that neglect of this spatial amplitude decay may lead to 
significant errors in response predictions. 

Figure 3 shows the dependence of the modal joint acceptance on 
k /k for various mode orders for a 10m long beam and for constant 

x m 

a . When the excitation wavelength is much greater than the 
structural wavelength (k x /k m << 1) in the presence of strong 
amplitude decay, j^ Is inversely proportional to m 2 . When 
k /k >> 1, i 2 is proportional to (k /k ) 4 and inversely propor- 
tional to m 2 . This is in agreement with the limits derived in 
Appendix A. 

For the present study, values for the spatial decay coefficients 
(a ,a ), for the coherence coefficients (c x ,c ), and for the 
average trace wavespeeds (U x ,U ) [or, alternatively, trace wave- 
numbers (k ,k )] were estimated from the Hamilton Standard pre- 
diction procedure [6], Aero-Commander tests [7,8], and ray acoustics, 
as is explained in more detail in Section 5. 

3.6.3 Reverberant Acoustic Fields 

The reverberant field joint acceptance is required both for calcu- 
lation of the acoustic power flow from the reverberant acoustic 
field inside the cavity (cabin) out through the fuselage structure 
(Eq. (53), and for calculation of the high-frequency internal 
radiation ratio [Eq. (52)]. 
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Mode 

Order 

a 

X 

c 

X 

m 

1 

0.5 

0.0 

] 

0.5 

0.01 

— • — 2 

0.5 

0.0 

— 4 

0.5 

0.0 

10 

0.5 

0.0 



Equation (27) is again used to compute the reverberant field 
joint acceptance. The representation for the co— spectrum of 
pressure field, in this case considered to be homogeneous over 
the fuselage surface, is well established [21] and is given by 

0 p (5J,x';(o b ) = C px (5,V = py (l,«> b ) <Pj 0 ( “b )> (56) 

where C (£,to h ) = sin (k£)/k£, 

C (n,w b ) = sin (kn)/kn, 

3? 3^ 

and <p 1 2 (to. ) > is independent of Xo for a homogeneous reverberant 
AO 0 

acoustic field. 

Then, the joint acceptance for reverberant excitation is estimated 
with [22]: 

„ rev 
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Similarly , 
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where I^(n), ^(n), and I^(n) are given by the above equations 
with n replacing m and L replacing L^. Si and Cin are the sine 
and cosine integrals £ 3 ]. 

3.7 Internal Coupling 

At very low frequencies, the acoustic power flow depends intimately 
on the modal characteristics of both fuselage structure and cabin 
volume. Then, detailed modal resonance frequencies and mode 
shapes are required to enable precise calculations to be carried 
out. In practice, computational efficiency requires the use of 
approximations in the description of both the cabin shape and 
its furnishings. For example, in the analytical model, the cabin 
volume is represented as a rectangular parallelepiped whose sur- 
faces can be deformed as shown in Figure 4(a) to represent the 
curvature of the sidewall and ceiling trim panels. 

This approach provides a more realistic model than either a 
cylindrical or an unmodified rectangular model, as can be judged 
by inspection of typical fuselage cabin cross sections (Figure 15). 
The analytical model assumes that the surfaces in the cabin are 
locally reacting, so that response of a given point on a surface 
will depend on the acoustic pressures at that point only. This 
is a reasonable assumption, since many surfaces in an airplane 
cabin are composed of materials such as carpets, thin plastic 
cloth backed by foam, and fiberglass-filled honeycomb with a 
perforated trim cover — which are locally reacting. The basic 
development of the perturbation method for the calculation of the 
cabin modes and resonance frequencies has been developed in [4], 

The resulting power inflow is calculated via the internal coupling 
factor f(n,r), which is discussed in detail below. At high 
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frequencies, where the power inflow depends only on acoustic 
modal density and where the cabin acoustic field tends to be 
diffuse, the internal coupling is formulated in terms of the 
radiation resistances of the fuselage structural elements, as 
demonstrated in Section 3 . 5 . 

3.7.1 Internal Coupling Factor 

The coupling of the fuselage vibration modes is expressed by the 
internal coupling factor f(n,r), defined in Eq. (21) as 

f (n,r) = | f <P n (x)ip r (x)dx 

where ijj r (x) represents the structural mode shape for mode of 
order r, and <J> n (x) is the cabin mode shape for mode of order n. 

The area of the structure that is radiating into the cabin is 
denoted by A. 

The general situation considered in the analytical model is shown 
diagrammatically for the (x,z) plane in Figure 4(b). The figure 
shows a volume of length L exposed to a structural mode that 
extends over a small portion of the structural length L. This 
represents the limited spatial extent of the structure response 
to a localized or nonhomogeneous excitation. 

The associated values of the coupling parameter in the x direction 

are given by the following relations [4], where the structural 

mode shapes are modeled with 4j m (x) = sin k x and the acoustic 

m 

mode shape is modeled by the zeroth order eigenfunction for the 
deformed volume viz. <f> (x) = cos k x.. The acoustic mode shapes 
and the coupling parameter are assumed separable in the x and y 
directions, as has been assumed for the structural mode shapes and 
the joint acceptance function. In particular. 
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f(n,r) = f(pq,mn) = f(p,m) f(q,n) 


where , for the x direction. 
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For the particular case of d = 0, then c^ = 0, 


1 fl - cos TT(m - k) , 1 - cos Tr(m + k) 
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and f pm = 0. Also, k = k p L x /7T and c d = d/L x . 


For the transverse axis. 
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where f and f qn are given by Eqs. (59b) and (59c) for f pm and 
f pm , respectively, with appropriate variable changes [4]. 

The factor f 2 (n,r) is synonymous with the joint acceptance function 
that couples the external pressure field to the structure. Values 
of f 2 (n,r) are nondimensional , positive numbers that are less than 
unity. Typical results are presented in [4]. 


This internal coupling factor is used only at low frequencies. 
When the acoustic modal overlap is much greater than zero, the 
approximation (Section 3*5.1) 


n „ ( *b> 


<e f 2 (n,r)> 


m 2 V 


n 




allows calculation of the interior coupling from band-averaged 
values of the joint acceptance of the structural modes. This is 
always the case for the study airplanes. 
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3.7.2 Radiation Resistances 

The radiation resistance of a structural system is determined 
by the relationship 


R ^ 
rad 


p . c . Acr , . 
i i rad 


( 60 ) 


where a d is the radiation efficiency, and A is the surface area 
of the structure. Using reciprocity, it can be shown that [24] 


a rad 


2A 

IT C ? 


r R 


(61) 


where < > indicates an average over the modes resonant in the 
band of concern. 


3 . 8 Test Conditions on Allowed Structural Wavenumbers 


To calculate the required joint acceptances and internal coupling 
factors, structural bending wavenumbers (or mode numbers) must be 
selected for each analysis frequency band Aw. The selection is 
simple in the case of isotropic flat rectangular panels since, 
for such panels, curves of constant frequency are (approximately) 
curves of constant bending wavenumber. If the structure is 
orthotropic or curved, or has pressurization stresses (as is the 
case for the study airplanes), the determination of allowable 
values of k and k n must be done with a frequency test, either 

w eAw 
mn 

or 


whichever is required. 
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The effects of curvature and pressurization and of the orthotropic 
properties of a structure over a frequency band are translated into 
the model by a set of resonance frequencies and corresponding wave- 
numbers for that band. It is the composite set for each band Am 
which is important. Imprecision, in the sense that a few modes 
may be classified wrongly, is unimportant. The estimate for input 
power will depend on the properties of the whole group rather than 
the specifics of a few individual classifications. If the modes 
can be classified approximately correctly as a group, then the 
desired result will be obtained in the input power calculation. 

The simple joint acceptances and internal coupling terms are 
entirely supported by general classification of modes. Equations 
for the prediction of resonance frequencies for the fuselage 
structures are discussed in detail in Section 3-8.1 below. 

The analysis bandwidth for calculation of resonant and nonresonant 
power flow is chosen to be sufficiently wide that at least several, 
but less than 100, structural modes are resonant within the computa- 
tional bandwidth. At low frequencies, bandwidths are typically 
1/2 or 1/3 octave, while as structural modes accumulate with 
increase of frequency, the bandwidth rapidly decreases (a 1/30 
octave bandwidth is typical at high frequencies). 

3.8.1 Resonance Frequency Equations 

The basic resonance frequency equation used for prediction purposes 
was derived by Mikulas and McElman [25] and is shown in Figure 5 
with the addition of a term due to pressurization. Figure 28 
shows the representation of a typical structural element. Other 
changes to the formulation were made to allow for boundary condi- 
tions other than simply supported. Use of the equation is described 
in Section 4.1., 
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z and z - are the distance to the centroid from skin middle surface for ring and stringer, respectively, 
r s 

X = mrrR/L where m = number of axial half wavelengths. D = Eh 3 /12(l - v 2 ), is the skin bending rigidity. 
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where n = number of circumferential full wavelengths. 


FIGURE 5. MIKULAS’ FREQUENCY EQUATION FOR A SIMPLY-SUPPORTED STIFFENED CYLINDER 
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3.9 Sidewall Representation 
3.9.1 General 

The analytical model predicts the bare fuselage noise reduction 
by equating the net acoustic power inflow from resonant and non- 
resonant structural modes to the power that is absorbed within the 
cabin interior. The model examines the modal response of the fuse- 
lage structure and of the cabin and calculates the coupling between 
the two systems to arrive at the internal noise levels and hence 
the noise reduction for the bare fuselage. In [5], the additional 
transmission loss provided by the sidewall treatment was calculated 
assuming infinite panel conditions and the results showed reasonably 
good agreement with test data for finite panels. Although the 
present analytical model calculates the power flow through finite 
structural panels in order to include both resonant and nonresonant 
structural response, it is assumed, on the basis of the above evi- 
dence and for simplicity, that the infinite panel approach can be 
used to calculate the additional transmission loss provided by the 
sidewall treatments. The infinite panel model for the sidewall 
noise transmission is discussed in this section. 

Figure 6 shows a simplified representation of the sidewall acoustic 
power flow. The external pressure field, characterized by incidence 
angle and intensity, acts on the airplane skin (element 1). The 
acoustic power flow W^ 2 from the skin to the cabin interior 
(element 3) when no sidewall (element 2) exists, is calculated 
from a knowledge of the bare fuselage transmission coefficient, 
and the transmitted power is then equated to the power absorbed 
within the cabin W to find the resulting internal diffuse 
acoustic pressure <p*>. The presence of the add-on sidewall 
acoustic treatment acts to reduce the acoustic power inflow to an 


-59 



FIGURE 6. 


POWER BALANCE REPRESENTATION OF SIDEWALL ACOUSTICS 




amount T g W 12 . In general, the transmission coefficient of the add- 
on sidewall acoustic treatment t cannot be calculated from a con- 

s 

sideration of the sidewall treatment in isolation, but rather a 
knowledge of the bare fuselage impedance is required as well. 


Calculation of the additional noise reduction provided by the 
add-on sidewall structure consequently involves (1) calculation 
of the noise reduction of the bare (infinite panel) fuselage struc- 
ture, (2) calculation of the noise reduction provided by the 
treated fuselage infinite panel representation, and (3) subtraction 
of the bare fuselage noise reduction from that of the treated fuse- 
lage to find the noise reduction for the add-on sidewall structure. 
Allowance must also be made for any changes in absorption coefficient, 
which may result from sidewall changes. 

The noise reduction for the untreated bare fuselage, for a diffuse 
interior acoustic field, can be expressed as [5] 


<p 2 > 

e 


NR = 10 log * 10 log l0 


(Sa> 0 p e°e 


<Pt^ ~1U AT p.C. 

1 P 1 


( 62 ) 


where <p 2 > characterizes the level of the external blocked pressure 
field on the fuselage surface, 

<p?> is the mean-square sound pressure in the cabin, 

(Sa) u is the absorption in the untreated cabin, and 

t is the bare fuselage transmission coefficient. 


Two symbols, S and A, are used for area to distinguish between 
transmitting area A and absorbing area S. 
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From a similar expression for the noise reduction for the treated 
fuselage (including the acoustic sidewall), the additional no s 
reduction provided by the add-on sidewall is [5] 

ANR = 10 log t - 10 log [(Sa)y/(Sa),p] (63) 

where t = t 9 /t is the transmission coefficient of the add-on 
S s p sldewa n acoustic treatment, 

T i S the transmission coefficient of the treated fuselage 
s2 

structure, and 

the subscript T refers to the treated cabin interior. 

Expressions for calculation of the transmission coefficents t s2 . 
and 1 and the interior absorption coefficients are, therefore, 
required . Sections (3.9.2) and (3-9.3) derive relations for T 
and T , respectively, while empirical data are used for absorption 
coefficients a 'although reliable analytical models are aval a t 


3.9.2 Sidewall Transmission 

The transmission coefficient T of a structure can be defined with 
respect to the blocked pressure field as [5] 



p t 

2 p l c l = 

P t 

il 

I — 1 

p i 

P 2 C 2 

P + P r 


p l°l 

P 2 C 2 


(64) 


where for massive structures and incident plane waves, the blocked 
pressure p x is approximately equal to the sum of the inciden P 


reflected (p r ) pressures 
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Figures 7 and 8 present the model of a double-wall sidewall structure 
whose transmission characteristics are of interest in this study. 

The untreated fuselage is shown also. 

Acoustic plane waves are assumed incident on the exterior of the 
structure and reflected and transmitted as determined by the 
various impedances present. The medium on the interior side of the 
structure is assumed to extend to infinity with an acoustic 
impedance of pc. Of interest is the transmission coefficient x 
and the noise attenuation across individual media or layered com- 
binations of such media when aligned in series. From expressions 
for the pressure ratios across boundaries between adjacent media 
and the pressure ratios across the media themselves, it is possible 
to derive expressions for the transmission coefficient across a 
complex series of layered media. Fundamental to this analysis are 
expressions for the characteristic impedance of the various media 
and for the impedance looking into a finite depth medium whose 
propagation and attenuation characteristics may be simply defined. 

It is now possible to derive expressions for the characteristic 
impedance Z 0 and propagation constant b of a porous medium and 
the impedance of a stiffened, pressurized and curved panel. 

Following [5], expressions for Z Q and b for a porous material are 
determined empirically from only the material flow resistivity. For 
example, empirical relationships exist for determining b for 
semi-rigid fibrous materials. Thus 

z o = “ b , and b = a + J 3 (65) 

where, following [-47] 
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EXTERNAL GEOMETRY OF PANEL-PRESSURE 
FIELD CONFIGURATION 


FIGURE 7. 




( 66 ) 


— 0 595 

a = (u)/c)[0.l89(pf7V ^ 

g = (a)/c)[l + 0.0978(917^) 3 

K is the complex compressibility of the gas in the porous 

material ( - atmospheric pressure), 
f is the frequency (u)/2tt), 

Y is the porosity of the porous material, 

c is the speed of sound in the gas in the porous material, 
p is the density of the gas in the porous material, and 
R is the flow resistivity of the porous material 
(0.01.4 pf/R-L - 1 in general). 

The effects of variations in material density are reflected in 
changes in material flow resistivity. Bies and Hansen 026] have 
recently extended the empirical relationships for a and 6 to values 
of pf/R-L much less than 0.01. 

It should be noted that this model does not consider mechanical 
transmission through the material fibers, such flanking of the 
acoustic path providing a limit to the noise reduction of a double- 
wall sidewall. 

Prom results such as [ 27 ], the impedance of an infinite stiffened, 
curved and pressurised fuselage structure for a sound wave incident 
at an angle 6 relative to the normal and at azimuthal angle <J> rela- 
tive to the x-axis, as shown in Figure 7, can be written as 

z (e,<|>) - wv^lxn + -H 1 - ^ ~ x)3 (68) 


where 
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x = O) 2 sin 4 e 2 A/Cc 4 ^) , 

'P = N x sin 2 Q ^ [cos 2 p + 2 sin 2 tj)]/(y 1 c 2 ) , and 
N = Ap R/2 (Ap acts positively outwards). 

X 

A = cos 4 <j) + 2D x y cos 2 (J) sin 2 c|) + D ■ sin 4 tj) 
n is the structural loss factor, 
k = oj/c , 


D x is the bending rigidity of the structure in section per- 
pendicular to the x-axis, D y is the y-axis' bending 

rigidity, and D the cross rigidity given by 
xy 


D 


xy 


(D V + D V + 4D. ) 
xx y y k 


and are the Poisson's ratios for the two axes, and 



Gh 3 

12 


where G is the modulus of rigidity for -shear stresses and 
h is the plate thickness. 


The effect of n is only important near the coincidence frequency. 
Then, the transmission coefficient t (= t^) across the untreated 
panel is given by Eq. (64) as 



PpCp 

Vl 
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Pp c p 

(-AA) 

Pi c i 




where Zp - PpCp/cos 9p and z^ is given by Eq. (68). 
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The transmission coefficient T g2 
given by 


across the treated sidewall is 


1 

T S2 



P 1 p 2 p 3 p 4 ^ 5 . 
P 2 *P3 P4 P5 p 6 


(70) 


The terms in Eqs. ( 69 ) and (70) are given by 
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cosh(b^d cos 0^ + i|>^) p,- z,_ 

p^ cosh ip Hi and P^ = ^ 


where Zg = jwpg + z g , Zg = pgCg/cos 0g, and 0 g = sin-^Cg sin 0^/c^). 


Then, the transmission coefficient of the add-on acoustic treatment 
T s (= t s2 /t ML ) is deflned as 



T s2 

T ML 


'£l\ 

k P 2 /UNTREATED 


'P- 


fP 


n-1 


n /TREATED 


(7D 


where p n = p , and the added transmission loss of the sidewall is 
found from 


ATL = -10 log 1Q x g . 


(72) 


In summary, the treated fuselage acoustic power flow into the 

cabin interior is found by multiplying the bare fuselage acoustic 

power flow into the cabin, calculated from Eqs. (51) and (52), by 

the transmission coefficient of the add-on acoustic treatment t , 

s’ 

computed using the above equations. The flow diagram for computa- 
tion, shown in Sea. 3.11, indicates the stage in the calculation 
procedure when this multiplication is carried out. 


3.10 Interior Acoustics 


3.10.1 Spatial Variation of Interior Noise Levels 

In order to generate a nonhomogeneous structural response to the 
nonhomogeneous excitation pressure field, the computation method 
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divides the fuselage structure into various elements located along 
the fuselage length. Resulting computations of the net acoustic 
power inflow produce a nonhomogeneous distribution of power inflow 
varying in much the same way as the excitation level varies over 
the fuselage exterior surface. Since it is desired to calculate 
the maximum noise level developed inside the airplane cabin, ac- 
count must be taken of this honhomogeneity in structural response 
and power radiation in calculations of the internal noise levels. 

An imaging technique is used to calculate the variation in interior 
noise levels. Each of the radiating structural elements is repre- 
sented as a point source mounted in the fuselage, located at the 
element center, and considered to radiate into the cabin volume. 

The contribution of each point source to the local pressure level 
along the cabin centerline is calculated by adding to the direct 
radiation from the point source itself, the contributions from all 
the image sources that are associated with reflections from the 
various bounding surfaces of the airplane interior (walls, ceiling, 

floor). 

The local intensity is the sum of the local intensities produced 
by each fuselage element (i.e., image-source array). The level 
variation along the fuselage interior is averaged to determine the 
space-average mean square pressure, which is also the basic output 
of analytical model already presented (Section 3-5). The local 
levels are adjusted at each frequency by adding the difference 
between the space-average noise levels calculated by the power 
flow and the image-source array models, respectively; thus, the 
space-average noise levels at each frequency calculated by the 
image array and the power flow model are equal. The differences 
in space-average spectrum levels between the two calculation schemes 
are less than 1.0 dB in most cases and result mainly from a certain 
coarseness in the representation of the image field. 
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3.10.2 Image Array Formulation 


Figure 9 shows the image array used to model the radiation pattern 
from a point source located at S in the sidewall. For model 
simplicity, the cabin fuselage is assumed to be rectangular. For 
example, the image (1,1) shown in Figure 9 represents the con- 
tribution from that ray which travels via the indirect path in- 
volving two reflections, one each from opposite wall and ceiling. 
Equation (73) below expresses the dependence of the local center- 
line intensity on distance from the source, on the absorption 
coefficients of each of the bounding surfaces (the sidewalls are 
considered identical), and on the paths by which the sound waves 
can travel from source to observer position [28], 
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(73) 


is the acoustic intensity at observer position 0, 


<p 2 > 

o 


is the observed mean-square pressure at the observer 
position. 


is the power radiated by the jth structural element 
(considered to radiate as a point source). 


a w ,a c ,a f are absor P tlon coefficient for sidewalls, ceiling, 
and floor, respectively, 

R is the axial distance between observer 0 and jth 

source location S, 

H,W are the mean height and mean width of the cabin. 
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p 

j 

r 

j 

r 

•' P,q,r are indices determining the path taken by acoustic 

ray from source to receiver, defined by 

' ( 23 - 1 for g > 0 q = (-1 ) Y (y + l)/2 

P = ) for y odd 

r 1-23 for 3 < 0 r = (q - 1) 

I . 

1 q = r = y/2 for y even or zero. 



r - The analytical model was computerized to automate the calculations 

i of the noise levels developed in prop fan-powered study airplanes 

An outline of the overall computational procedure, with details of 
j how the various components of the analytical model are assembled 

together within the computer programs, is presented in this section. 

r*- 

') 

f 

Figure 10 presents a flow diagram of the computational processes. 
f" The various computational phases and the associated input require- 

i ments are linked together; references to appropriate sections of 

r the report text and/or relevant equations are included; locations 

) where important decisions must be made are indicated. 

j The computing prediction system in fact consists of five inter- 

dependent computer programs. Four of these organize acoustic and 
j" structural data and perform preliminary calculations, the results 

of which are used as input to a main computer program for final 
p calculation of the cabin interior noise levels. 

r- Tt should be noted that, while the effects of changes to the struc- 

[ ture require that all programs be run, changes to the sidewall trans- 

mission characteristics and to the cabin interior can use data stored 
on existing input files generated by previous (baseline) computer 
runs. 
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INPUTS 

(Computed or Empirical) 


OPERATION 


TEXT . REFERENCES 


Structural Properties and Damping 
Characteristics 
as per Figure 6 


Excitation Characteristics: 

Spatial Variation 
Trace Velocities 
Coherence 

(Inboard and Outboard Propellers) 


Select Structural Representation: 
Fuselage Configuration 
Structural Segmentation 


Sections 4.5. 1, 4.5.2, 4.5.3 


— Section 3.8.1 



Frequency-Dependent Model for Structure: 
Calculate Structural Resonance 
Frequencies & Classify Modes 






Calculate Band-Averaged 
Response Functions: 

Joint Acceptances, Modal Densities 
Radiation Efficiencies, Loss Factors 


Sections 3.8, 5. 1 


. — — Sections 4.3, 4.6 


Sections 3.6, 5.5, 5.6 


Cabin Properties and 
Dimensions 
(Tables 1 and 3) 


Determination High and 
Low Frequency Regimes 
Using Modal Overlap Criterion 


Section 3.7 


Calculation Cabin Acoustic Resonance 
Frequencies and Internal Coupling 
Factors for Each Frequency 



1 . 

Sidewall Transmission Loss 
(Figure 20) 


Acoustic Power Flow Calculations*. 
Resonant and Non-Resonant, 
Spatial Variation 
Along Cabin Length 


1 







Cabin Absorption 
(Figure 21) 



Formulate Power Balance Equation: 
Equate W nef and W Qbs 






Solve for Internal 
Space-Averaged SPL 


Compute Axial 
Variation in SPL 


Compute Space-Average & Axial 
Variation in A-Weighted SPL 


Section 3.5.3 
(Equations (51) and (52)) 


Section 3.5.5 

(Equations (52), (54) and (2< 


Section 3.10 
(Equation (73)) 


FIGURE 10. FLOW DIAGRAM FOR COMPUTATION OF INTERNAL NOISE LEVELS FOR STUDY AIRPLANES 
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4.0 BASELINE AIRCRAFT 


The previous section described development of the analytical 
model which is to be applied to the design of fuselage side- 
walls for the baseline prop fan-powered study aircraft such that 
the interior noise levels will meet an 80-dB criterion for 
the average A-weighted sound level during cruise. The fuselage 
sidewalls of the study aircraft are to be designed to have 
minimum surface weight density and high sound transmission 
loss. Two concepts are to be considered for reducing low- 
frequency propfan noise through modification of the fuselage 
sidewall design: (1) add-on noise control features wherein 
nonstructural noise-control elements are to be added to the 
conventional sidewall structure, and (2) advanced noise 
control designs wherein the materials and the basic structural 
configuration of the baseline aircraft could be modified to 
obtain higher sound transmission loss and hence lower interior 
noise levels than could be obtained by add-on designs, possibly, 
at a lower total weight penalty. 

Because the effectiveness of add-on or advanced acoustical 
treatment depends on the characteristics of the fuselage 
structure to which the additional acoustical treatment is 
applied, it was necessary to select baseline aircraft that 
would allow meaningful assessments of the application of the 
analytical method. The study required that three baseline 
aircraft be selected, namely a wide-body aircraft, a narrow- 
body aircraft, and a small-diameter aircraft. The fuselage 
diameters were to be as follows: 

Airplane Fuselage Diameter 

4.88 to 6.10 m ( 1 6 to 20 ft). 

3.05 to 3.96 m (10 to 13 ft) 

1.83 to 2.44 m ( 6 to 8 ft) 


Wide-body 

Narrow-body 

Small-diameter 
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It was considered desirable, as well as acceptable, that the 
baseline study aircraft should be derived from aircraft that 
were already designed or flying. This section presents the 
general requirements for the existing-design aircraft as well 
as the baseline study aircraft. The process of designing the 
propfans to power the baseline study aircraft is then describe . 
This is followed by a description of the relevant structura 
characteristics of the baseline study aircraft, and a dis- 
cussion of the acoustical characteristics of the fuselage 
sidewall and interior. Finally, the properties of the ex- 
terior noise field are presented. 

4. 1 Study Requirements 

In addition to the requirement to have fuselage diameters in 
the range described above, the baseline aircraft were to have 

•wing-mounted engines, 

•an initial-cruise altitude of 9000 m ( 30,000 ft) and 
•a cruise Mach number of 0.8. 

Also, the wide-body and narrow-body aircraft were to be able 
to carry at least 50 passengers over a range of at least 805 

(50 0 mi). 

The propfans were to be designed using data developed by the 

Hamilton Standard Division of United Technologies Corporation, 

[ 7 ]. The fuselages of each of the three baseline study air- 
craft were to have structural characteristics compatible wi ^ 
cruising at the specified cruise altitude. The structures o 
the three baseline study aircraft were expected to have some 
differences because of the different design approaches used y 
the manufacturers of the existing design aircraft to achieve 
f a ii-safe structural requirements. 
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4. 2 Choice of Ex i s ti ng-Des i gn Aircraft 


The following three aircraft were selected, as being able to 
meet the general study requirements outlined above. They were 
the McDonnell Douglas DC-10-30 for the wide-body aircraft, the 
Boeing 707-320B for the narrow-body aircraft, and an aircraft 
having the fuselage of a Grumman Gulfstream I, but modified 
for high-speed subsonic cruise with a swept wing similar to 
that on a Gulfstream II instead of the unswept wing on the 
Gulfstream I. 

Each of the three existing-design aircraft meet all the re- 
quirements for initial cruise altitude, cruise Mach number, 
payload, and range [48,49], Also, each aircraft has a fuselage 
structure capable of carrying all design loads, which there- 
fore forms a realistic fuselage design for a baseline study 
aircraft . 

Each of the existing-design aircraft was modified, conceptually, 
to accept the propfan propulsion system. The DC-10 is powered 
by three turbofan engines - one on each wing, supported below 
the wing from a pylon, and one mounted through the vertical 
stabilizer. The aft fuselage of the DC-10 was modified to 
eliminate the center engine. The two wing engines were re- 
placed by propfans, and an extra two propfans were added to 
the wings. The 707 had the four turbofan engines replaced by 
four propfans. The small-diameter airplane used two propfans 
located approximately at the spanwise wing station where the 
propellers are mounted on the straight wings used with the 
Grumman Gulfstream I aircraft. All aircraft are low wing 
designs. The propfans were installed such that the center of 
thrust was above the wing rather than below the wing as it is 
for the 707 and the wing-mounted engines on the DC-10. 


:: rr — 

body aircraiu ^ would have reduced the 

have used more than four engines. That would ^ englne 

diameter of each of the propfans, and if rance would 

iocation had remained fixed, the prope ler t P cle 

j m hpn the exterior noise field impinging 

practical. 

4 . 3 Propfan Sizing Analysis, 

, .. n 2 3] have considered the selection of 

Previous ® u ^ a ’ glven aircraft mission. In those 

studies 8 considerations were given to aerodynamic, mechanical 
„d cosmic factors in producing an optimum propfan design 
f or a specified mission. For the study here, however, it 

not feasible to conduct such an analysis since 
was not feas conceptually, on existing-design 

fans were to be installed, r propfans 

aircraft. Therefore, the approach used to si . 

for the three haseUne study aircraft proceeded from a defini 
tion of the total thrust required to propel the aircraft 

the desired cruise Mach number and initial cruise altitud 

values for total installed thrust per airplane, obtained 
the basis of previous analyses, are listed be ow. 


Airplane 


Wide-body 

Narrow-body 

Small-diameter 


Installed Net Cruise Thrust 
Per Airplane 

133.45 kN (30,000 lb) 
71.17 kN (16,000 lb) 
22.24 kN ( 5,000 lb) 


- 78 - 


These values of total net thrust per airplane, coupled with 
the number of propfans per airplane, defined the net thrust 
to be produced by each propfan. 

The Hamilton Standard data in [ 7 ] describe performance char- 
acteristics for propfans having 6, 8, or 10 blades. For a 
propfan producing a given net thrust, a study of the perform- 
ance data in [ 7 ] indicated that the smallest-diameter prop- 
fan would be one having the largest number of blades. Since 
data were only available for 6, 8, or 10 blades, propfans 
designed for the three baseline study aircraft were based on 
propfans having 10 blades. 

The next design parameter that was selected was the operating 
tip speed of the propfan. The tip rotational speed V rot , was 
chosen to be 243.84 m/s (800 ft/s) . This tip speed would be 
used during takeoff to minimize the noise level at the Appen- 
dix C noise-measurement location for noise certification under 
Part 36 of the Federal Aviation Regulations. That same tip 
speed would be used during cruise. 

With specifications for.(l) the installed net thrust per 
propfan, (2) the number of propfans per airplane, (.3) the cruise 
tip speed, and (4) the number of blades per propfan, it was 
necessary only to select one more design parameter in order 
to calculate the rest of the propfan characteristics. 

That design parameter was the cruise power loading P/D 2 where 
P is the power supplied to the propfan and D is the propfan 
diameter. From the studies that had been previously con- 
ducted of propfan installations (see Refs. 1 , 2,3), and from 
discussions with personnel at Hamilton Standard, it was 
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decided to select a P/D* ratio of 280 kW/m* or 35 SHP/ft* 
where SHP is the propfan shaft horsepower. The same value of 
power loading was assumed to apply to the propfans for each 
of the three baseline study aircraft; it would permit fuel- 
efficient cruise at the specified cruise altitude and cruise 

Mach number. 

It was further assumed that the atmospheric conditions at the 
design cruise altitude would be those of an International 
Standard Atmosphere - namely a density ratio of p Q /p - 2.6686 
and a speed of sound c = 303.2m/s where Pq is the density of 
the air at sea level and p is the density at 90.14 . 4 m. For a 
cruise Mach number of 0.8, this assumption gave an airspeed, 

V of 242.58 m/s. The helical tip Mach number was then 


M 


hel 



1.13 


Propfans on each study aircraft were powered by turbofan 
engines operated as turboshaft engines to provide shaft power 
through a reduction gearbox. The turboshaft engines pro- 
duced some jet thrust to augment the thrust produced by the 
propfans. Previous analyses had indicated that the jet thrust 
was about 10 percent of the total propfan thrust. Those anal- 
yses had also indicated that the sum of the drag resulting 
from flow around the nacelle and the gearbox , and the drag 
resulting from interference effects between the wing and the 
nacelle and gearbox, was about 6 to 8 percent of the. total 
propfan thrust. Therefore, it was assumed that the jet thrust 
was approximately equal to the sum of the nacelle' drag and 
the wing/nacelle interference drag. With this assumption, the 
net installed thrust from the propfan equalled the net 
uninstalled thrust, F n ^ u » 
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The propfan performance analyses in [ 7 ] are based on the 
generalized non-dimensional parameters used for analyzing the 
performance of propellers. The parameters are: the power 

coefficient, C p , the advance ratio J, and the thrust co- 
efficient C p . A net efficiency parameter, n net > is also 
defined. 

The propfan design process consisted of the following steps: 

1. Determine the advance ratio J; 

2. Determine the power coefficient C p ; 

3. Prom C p and J, find the thrust coefficient 
Cp and the net efficiency n net ; 

4. Prom the thrust coefficient and the net thrust, 
find propfan diameter; and 

5. Find propfan rotational speed from the ND product 
and propfan diameter. 

The analysis was carried out only for the design cruise Mach 
number of 0.8 at the design cruise altitude, for propfans 
having 10 blades. 

Advance ratio is defined as 

j = v a /( v ro t // ' rr ^ = 3,13 for P resent conditions; 

Propfan power coefficient C p is defined as 
C p = power/pn 3 D 5 , 

and can be related to propfan design parameters by 
C p = (P/D 2 ) (p Q /p) / [5.671 (ND/1000) 3 ] 
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where P Is in kW, D is in m, N is in rev/min, and e Q is 1.225 
kg/m 3 . 

The ND product is, 

ND = 60V A = 4657 m/min. for present conditions, 
rot 

Using P/D* = 280 kW/m*, the power coefficient for the propfans 
is then 

C p = 1.30. 

Propfan performance data from Table VII of [7(0] were used to 
develop the design chart shown in Figure 11 which relates 
power coefficient, advance ratio, thrust coefficient, and net 
propulsive efficiency. Entering the design chart at a power 
coefficient of 1.30 and interpolating at the advance ratio 
of 3.13 give's an uninstalled net thrust coefficient of 0.355 
with a net propulsive efficiency of 0.85, the efficiency para- 
meter being defined as 

n net = J C F /C P 
n n,u 

An efficiency parameter of 0.85 is compatible with values that 
were used in previous propfan design studies and is close to, 
but does not equal, the maximum-possible efficiency for the 
design advance ratio 3.13. The maximum net efficiency accor - 
ing to the data in [ 7 ] would be slightly higher than 0.85. 

To achieve more efficiency, however, would require operating 
the propfan at a lower power coefficient which would mean a 
lower thrust coefficient and hence a larger propfan diameter 
in order to provide the same net thrust. 
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Advance Ratio, J 



/ 


^net " /^p 



Ref: Table VII of [7(c)] 





Propfan diameter was determined from the net thrust coefficient 

with 

F 

n 

0 = net thrust/pn 2 D 4 

^n 

where net thrust is in newtons, and p, n, and D have the units 
stated above. 

The equation for C p can now be written as ■ 

n 

C = (F ) (p n /p)/[340. 28 (ND/1000 ) 2 D 2 ] 

F n 0 

n 

where F is the net thrust in newtons required per propfan, 
and other parameters have units specified earlier. 

From previously given net thrusts per aircraft, the installed 
net thrusts per propfan are: 


Wide-body: 33,362 N (7500 lb) 

Narrow-body: 17,793 N (4000 lb) 

Small-diameter: 11,120 N (2500 lb) 

Then, on substituting values for p q /p and the ND product, the 
propfan diameter D can be calculated for each value of net 
thrust. Rotational speed N is found from ND/D = 4657/D. 

4 . 4 Characteristics of Baseline Air craft 

The baseline study aircraft were designed by placing the 
propfans on the wings at appropriate locations. For wide-body 
aircraft, the inboard propfans were placed at the same wing 
station at which the wing engines are currently installed on 
the DC-10-30. The outboard engines were placed at an 
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appropriate location farther outboard. For the narrow-body 
aircraft, the propfans were placed at the same stations as the 
existing locations for the wing engines. For the small- 
diameter aircraft, engines were installed approximately at 
the same wing stations that the propellers are installed for 
the Grumman Gulf stream I aircraft, but on the swept wing. 

Results of the propfan sizing study for the three baseline 
study aircraft are shown in Table 1. Each of the aircraft 
meets the general study requirements described earlier. The 
data in Table 1 indicate that the ratio of propfan diameter 
to fuselage diameter is approximately equal to, or somewhat 
greater than, 1. ; 

On the two larger airplanes, tip clearances range from 0.4 
propfan diameters for the inboard propfan of the wide-body 
aircraft to 2.7 propfan diameters for the outboard propfan on 
the narrow-body airplane. Tip clearance for the small dia- 
meter airplane is 0.3 propfan diameters. 

The sweep of the wings for the wide-body and narrow-body air- 
planes means that the outboard propfans, even though they are 
significantly farther outboard, could cause the exterior noise 
field to be spread over an extensive region of the fuselage 
and hence eliminate the option of using localized noise control 
methods in a particular area of the cabin (as is discussed in 
Section 7.0). Localized procedures used in the past on pro- 
peller-powered airplanes have located coat closets in the 
regions of the plane of the propeller in order to put a com- 
patible cabin furnishing in the region of highest noise levels. 
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Table 1. 

PERFORMANCE PARAMETERS FOR 
PROPFAN-POWERED BASELINE 
STUDY AIRPLANES 



Study Airplane 

Parameter 

Wide-body 

Narrow-body 

Small- diameter 

Design range, km (n, mi) 

7400 (4000) 

6900 ( 3700 ) 

6670 (3600) 

Maximum takeoff weight, kg 

252,000 

148,000 

30,000 

(lb) 

(555,000) 

(327,000) 

(66,000) 

Typical number of psgr seats 

270 

145 

18 

Fuselage diameter, <j> m (ft) 

6.02 ( 19 . 75 ) 

3.76 (12.33) 

2.44 (8.0) 

Fuselage length, m (ft) 

52.0 (170.5) 

45.7 (150.0) 

21.7 (71.3) 

Number of propfans 

4 

4 

2 

Net cruise thrust per propfan, 

F , N (lb) 
n 

33,362 (7500). 

17,793 (4000) 

11,120 ( 2500 ) 

Propfan diameter, D, m (ft) 

5.83 (19.1) ■ 

4.25 (14.0) 

3.37 (11.1) 

Propfan diam. /fuselage diam. 

0.97 

1.13 

1.38 

Propfan rotational speed, N, rpm 

799 

1096 

1384 

Power required per propfan. 

' 



P, kW (hp) 

9517 (12,762) 

5058 (6782) 

3172 (.4254) 

Blade-passage frequency, f ., Hz 

133 

183 

231 

Inboard tip clearance, r/D 

0.4 

1.3 

0.3 

Outboard tip clearance, r/D 

CO 

i — 1 

2.7 

— 








_U 




J 


j i 


J 






The data in Table. 1 also show the power required per propfan and 
indicate that it ranges from approximately 10 megawatts per 
propfan for the wide-body airplane to 3*2 megawatts for the 
small-diameter airplane. The table also lists the maximum 
takeoff gross weights and design ranges for the three airplanes. 
Typical numbers of passenger seats are given for mixed-class 
configurations for the wide-body and narrow-body airplanes; 
the small-diameter airplane would have 18 seats in a typical 
business/executive-j et arrangement. 

Table 2 lists the values of performance parameters that are 
common to each of the three baseline study airplanes. The 
values in Table 2 were either derived as indicated above, or 
were specified as part of the study requirements. Figures 
12, 13 , and 14 show the three baseline study airplanes and 
illustrate the relative fuselage and propfan diameters, the 
placement of the propfans on the wings and the clearances 
between the fuselage and the propfans. Note that, although 
the design of the propfans and the installation on the air- 
craft is considered to be reasonably realistic, the airplane 
designs were not developed as a result of a detailed design 
study and thus are somewhat conceptual. However, the designs 
are considered realistic enough to provide suitable baseline 
aircraft for development and initial application of the 
analytical model. 
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PERFORMANCE PARAMETERS COMMON 
TO EACH BASELINE STUDY AIRPLANE 


Initial cruise altitude, m (ft) 

9014.4 (. 30 , 000 ) 

Cruise Mach number, M 

0.80 

Cruise airspeed, V a ,m/s (knots) 

242.58 ( 471 . 6 ) 

Speed of sound, c, m/s (knots) 

303.23 ( 589 . 5 ) 

Density ratio, p Q /p 

2.6686 

Number of propfan blades 

10 

Propfan tip rotational speed. 

243.84 ( 800 ) 

V rot’ m/s 

Propfan helical tip Mach number, M^ el 

1.13 

Propfan power loading, kW/D 2 (SHP/D 2 ) 

280 ( 35 ) 

Propfan advance ratio, J 

3.13 

Propfan power coefficient, C^ 

1.30 

Propfan uninstalled net thrust 

0.355 

coefficient, C F 

n,u 

0.85 

Propfan net propulsive efficiency, n net 



1 



FIGURE 







W-BODIED AIRPLANE 


.J - J .__.J J -- 


J 




Figure 15 illustrates typical seat arrangements and cabin 
interior contours for the three study aircraft. The wide— body 
and narrow-body seating arrangements are fairly easy to define 
as there is little variation among commercial aircraft. The 
wide-body aircraft has two aisles in a first-class and coach- 
class configuration whereas the narrow-body aircraft has a 
single aisle. 

Seating arrangements for the small-diameter fuselage are less 
easy to define as there are many configurations used in 
current business/executive jet aircraft. However, there will 
be a single aisle down the center of the aircraft, as shown 
in the fuselage cross-section. 

4.5 Baseline Fuselage Designs 
4.5.1 General Considerations 

The fuselage structures of the three baseline aircraft are of 
conventional skin-stringer-frame construction, an example of 
which is given in Figure 16 . The figure shows longitudinal 
stiffeners, or stringers, with hat sections and circumferential 
stiffeners, or frames, with zee sections. The frames are 
coupled to the skin by shear ties, and to the stringers by 
clips. Fail-safe straps or doublers are bonded to the skin 
at the frame locations to provide local reinforcing of the 
skin. These four structure elements — skin, doubler, 
stringer, and frame — are all taken into account in the 
descriptions of the baseline structures. 

Other parameters to be considered are fuselage radius, cabin 
windows, and in-plane stresses in the skin. These stresses 
result from the cabin pressurization and fuselage bending 
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r 


FIGURE 15. TYPICAL PASSENGER CABIN ARRANGEMENTS FOR 
THE STUDY AIRPLANES 







FIGURE 16 TYPICAL SK I N-ST Rl NG ER- FRAM E CONSTRUCTION 
FOR CONVENTIONAL FUSELAGES 
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loads, and form a complicated distribution with tensile 
stresses on the upper portion of the fuselage and compressional 
stresses on the lower portion. An example of such a stress 
distribution Is given In [29] for the Boeing 7^7. 

4.5.2 Baseline Structures 

Baseline structures for the three study airplanes were defined 
on the basis of existing technology and typical characteristics 
of current-day aircraft of comparable size. Dimensions and 
shapes for the stringers, frames and doublers for the study 
aircraft are shown in Figures 17 , 18 , and 19 for the wide-body, 
narrow-body, and small-diameter fuselages respectively. All 
airplanes are assumed to have hat-section stringers. The two 
larger diameter fuselages have zee-section frames, and the 
small-diameter fuselage has cee-section frames. All fuselages 
are assumed to have doublers beneath the frames, although 
shear ties are included only in the wide-body fuselage. Skin 
thickness is assumed to be constant for the wide-body and 
small-diameter fuselages, but, since some aircraft have ^ a 
skin whose thickness varies along the length of the cabin, a 
varying skin thickness was assumed for the narrow-body baseline 

structure. 

Characteristics of the baseline fuselages, including skin 
thickness, frame and stringer spacing, cabin pressure differ- 
ential, and material properties of the structure are listed 

V- 

in Table 3 . 

In order to reduce the complexity of the analytical model, the 
structure is assumed to be uniform in construction throughout 
the region above the cabin floor (acoustic power flow through 
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(a) FORWARD OF OUTBOARD PROPELLER 


SKIN THICKNESS: 1.78 mm 


FRAME 


STRINGER 


20.7 

K — H 




(b) AFT OF OUTBOARO PROPELLER 


SKIN THICKNESS: 1.78 mm 


STRINGER 


20.7 


i».< 

1.78 


All Dimensions in mm 


FRAME 


4.1 




FIGURE 17. STRUCTURAL DETAILS FOR BASELINE 
WIDE-BODY FUSELAGE 
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Distance From Forward End of Cabin (%) 


STRINGER FRAM£ 

25.4 

H — H 



All Dimensions in mm 


FIGURE 18. STRUCTURAL DETAILS FOR BASELINE NARROW-BODY 
FUSELAG E 



SKIN THICKNESS = 1.02 mm 


FRAME 


STRINGER 
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FIGURE 19. 


STRUCTURAL DETAILS FOR BASELINE 
SMALL DIAMETER FUSELAGE 
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TABLE 3 

CHARACTERISTICS OF BASELINE FUSELAGES 


AIRPLANE 

WIDE 

BODY 

NARROW 

BODY 

SMALL 

DIAMETER 

FUSELAGE DIAMETER 

(m) 

6.02 

3-76 

2.44 

CABIN LENGTH 

(m) 

41.45 

34.00 

10.06 

CABIN WIDTH 

(m) 

5-72 

3.55 

2.18 

CABIN HEIGHT 

(m) 

2.44 

2.31 

1.76 

SKIN THICKNESS 

(mm) 

1.78 

1.02 to 3 • 05 

1.02 

FRAME SPACING 

(m) 

0.508 

0.508 

0.330 

STRINGER SPACING 

(m) 

0.183 

0.229 

0.127 

CABIN PRESSURE DIFFERENTIAL 

(N/m 2 ) 

5.515 x 10 5 

5.515 x 10 5 

5.515 x 10 5 

MATERIAL 


Aluminum 

Aluminum 

Aluminum 

YOUNG'S MODULUS E 

(N/m 2 ) 

7.2 x 10 10 

7.2 x 10 10 

7.2 x 10 10 

DENSITY 

(kg/m 3 ) 

2770 

2770 

2770 




the structure beneath the floor Is assumed to be negligible 
In the analytical model). This means that the analytical model 
does not take Into account circumferential variations In 
skin thickness, stringer spacing, or frame depth. Furthermore, 
the model excludes windows and window cut-outs. A more- 
detailed representation of the sidewall structure in the analy- 
tical model is not feasible within this study. 

4.5.3 In-plane loads 

The fuselage structure is subjected to in-plane tensile/com- 
pressional loads N x and N y , in the axial and circumferential 
directions respectively, and to shear loads N xy . '■^ e ^ w0 
main contributors to these loads are pressurization and body- 
bending forces. If the fuselage radius is R and the pressure 
differential is Ap, then the cabin pressurization loads intro- 
duce contributions of 0.5 R.Ap and R.Ap to N x and N y respect- 
ively, the shear load contribution being zero. In all calcu- 
lations the fuselage pressure differential was chosen to be 
55 . lkPa (8 p.s.i.). Body-bending loads vary with airplane 
flight conditions and with axial and circumferential location 
on the fuselage. Also, they will be affected by the weight of 
the aircraft. At some locations the loads will be positive 
(tensile) whereas at other locations they will be negative 
( compressional) . 

Because of the variation in the magnitude of the body-bending 
loads, the influence of these loads on structural response 
has been neglected in the present analytical model. Pressuriza 
tion loads have been included as being indicative of an average 
value throughout the fuselage structure. In any case, results 
presented in Section 6.0 indicate that in-plane loads are not 
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important at low frequencies, which are critical from a noise 
control point-of-view, because of the effects of the frames 
and stringers. 

Circumferential pressurization stresses in the fuselage skin 
will be relieved to some extent by the frames and doublers. 
Typical average values for the percentage of nominal loop stress 
accepted by the skin lie in the range of 75 to 80 % . However, 
the analytical model assumes that the full pressurization loads 
are carried by the skin since the effect on calculated panel 
resonance frequencies is small. 

4.5.4 Noise Control Treatment 

Two properties of the noise control treatment have to be 
defined for the baseline airplanes. The first is the trans- 
mission coefficient for the sidewall treatment and the second is 
the acoustic absorption in the cabin. 

The baseline transmission loss spectrum was obtained from a 
review of empirical data for sidewall treatments typical of 
those in present-day use. These sidewalls consist mainly of 
glass fiber batts and a trim panel. In some cases damping 
tape may be used and in other cases lead- impregnated vinyl 
sheets may be installed, but such items are usually introduced 
for specific problems. These special treatments are excluded 
from the baseline sidewall. 

Much of the suitable transmission loss data is associated with 
tests in laboratory facilities where the fuselage structure is 
represented by a stiffened flat or curved panel, and where 
the transmission loss of interest is the insertion loss provided 
by the add-on treatment. It is this insertion loss which is of 
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interest to the present study because the noise transmission 
of the basic fuselage structure is calculated by the analytical 

model . 

The baseline sidewall is shown diagrammatically in Figure 20. 

It is assumed that acoustic insulation material is installed 
over the stringers and between the frames. Then a secon 
layer of the material is placed as a continuous blanket 
over the top of the frames. Finally, the trim panel is in- 
stalled such that there is a gap between the trim and the frame 
cap, the gap being filled by the glass-fiber batt. The glass 
fiber material is assumed to have a density of 9-61 kg/m ^ 

(0.6 lb/ft 3 ) and the trim panel a surface density of 1.76 kg/m 
(0.3 lb/ft 2 ). The total thickness of the batts is assumed to 
be 10.2 cm (4 inches) for all three study airplanes. In the 
case of the narrow-body and small-diameter designs, the region 
between fuselage skin panel and trim panel will be completely 
filled with the glass fiber material. The wide-body airplane 
will have an airgap between the batts. Installations of this 
type are described in [30] for a small-diameter (Gulf stream II) 
airplane, in [31] for a narrow-body (Boeing 727) airplane 
and in [32] for a wide-body (DC-10) airplane. The associated 
insertion loss, or additional transmission loss, associated 
with the baseline installation is shown in Figure 20. The 
curve is based on data such as is given in [ 30 ] and [33]. 

Acoustic absorption in the passenger cabin is provided by the 
sidewall trim panels (including the ceiling), carpet, seats, 
and other surfaces. Absorption data for materials used in air- 


Transmission Loss re Bare Panel dB 



Frequency Hz 


. BASELINE TRANSMISSION LOSS FOR ADD-ON SIDEWALL 


FIGURE 20 



craft interiors are not readily available but estimates can be 
made of average absorption coefficients, or corresponding 
absorption units, provided by the different items. Estimated 
absorption coefficient spectra for typical sidewall, carpets, 
and seats are given in [34], where the coefficients for the 
seats are based on typical seat distributions as a function of 
wall area. The absorption coefficients in [34] were obtained 
from laboratory test data for the individual items rather than 
for measurements in furnished aircraft cabins. However, there 
are some unpublished data for measured space-averaged absorption 
coefficients in furnished cabins and these data provide the 
band of data presented in Figure 21. 

Using the data in [34 ] and the band of measured data in 
Figure 21, a series of component absorption coefficient spectra 
were constructed and these are shown in Figure 22 . The com- 
ponent spectra were chosen so that the space-averaged absorp- 
tion coefficient a, for a typical airplane cabin, would lie 
roughly at the upper boundary of the measured data in Figure . 
The calculated a spectrum is shown in Figure 21 for comparison 
with the measured values. The main difference between the 
data in Figure 22 and the data in [34], is in the assumed ab- 
sorption coefficient for the sidewall trim. This spectrum 
was modified to incorporate results of Goss [ 3 0] and to bring 
the composite spectrum into agreement with the measured values 

O f a . 

Also shown in Figure 21 Is an absorption coefficient spectrum 
based on reverberation measurements in a bare fuselage [ 35 ]. 
Since the absorption coefficient was obtained in this. case 
using conventional reverberation equations the value includes 
the effect of noise transmission out of the fuselage interior. 
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Average Absorption Coefficient a 



Frequency Hz 


FIGURE 21. COMPARISON OF MEASURED AND MODEL ABSORPTION COEFFICIENTS 



Absorption Coefficient 




A typical value for the absorption coefficient for bare 
aluminum, excluding the effects of noise transmission is 0.05, 
as was used in [ 4 ] . 

4.6 Propeller Noise Field 

In the development of the analytical model it was shown, in 
Section 3.3, that the inhomogeneous pressure field generated 
by the propeller can be described in terms of the spatial 
distribution of the pressure amplitude or level, the spatial 
decay of the coherence function, and the trace velocities of the 
pressure field in the axial and circumferential directions. 

This section will describe briefly the characteristics of the 
acoustic fields used to represent the exterior pressure fields 
on the three study aircraft- 

The basis for calculation of the spatial distribution of the 
pressure level is the prediction procedure developed by 
Hamilton Standard [6]. Typical free-field directivity curves 
[6] for the near-field pressures are shown in Figure 23 for 
a range of tip clearances similar to those for the study air- 
planes. The available data refer to an 8-blade propeller with 
a tip rotational Mach number of 0.7, but the general trend of 
the data can be used for the selected 10-blade propeller with 
a tip rotational Mach number of 0.8. 

Predicted variation of the overall sound level along the 
fuselage exteriors of the three baseline aircraft are shown 
in Figure 24. In calculating these surface pressure distribu- 
tions, it was assumed that pressure doubling (+6dB) occurred in 
the neighborhood of the plane of rotation of the propellers and 
that the pressure increase was +3dB elsewhere. The curves in 
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Fore and Aft Location in Fan Diameters 


RECTIVITY AS A FUNCTION OF TIP CLEARANCE 




Overall Sound Pressure Level, dB re 20 >t/P< 



FIGURE 24. VARIATION OF EXTERNAL SOUND PRESSURE LEVEL 
ALONG THE CABIN 
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Figure 24 for the wide-body and narrow-body airplanes show 
the double-hump characteristic associated with the two pro- 
pellers on each side of the fuselage. The small-diameter 
design shows only a single peak as there is only one propeller 
on each wing,. The predicted peak pressure level is a function 
of propeller efficiency, tip rotational speed, cruise Mach 
number, tip clearance, altitude and number of blades. Since 
all these parameters, with the exception of tip clearance, 
are the same for all the study airplanes, the peak overall 
sound pressure levels shown in Figure 24 differs from airplane 
to airplane as a function of tip clearance only. The peak 
level is highest (l46dB) on the small-diameter fuselage, 
where the tip clearance is smallest (see Table 1 ) and is 

lowest (137dB) on the narrow-body airplane where the tip 
clearance is greatest. 

The data in Figure 24 show the rapid decrease m sound level 
as distance increases from a region of peak level. This 
spatial decay provides values for amplitude decay parameter a x , 
these values varying with position along the fuselage. 

The spatial decay of sound pressure level in the circumferential 
direction was estimated by first calculating the free-field dis- 
tribution using Figure 25(b). Then the influence of the fuse- 
lage structure was taken into account by adding 6 dB at the 
position of minimum tip clearance (location B m Figure 25(a)) 
decreasing to 3 dB at locations on the vertical centerline, for 
example, location D in Figure 25(a) . 

The analytical study was required to consider three different 
spectrum shapes for the excitation field on each of the study 
airplanes. These spectra, which refer to the relative levels 
of the different harmonic components of the propeller noise 
field, are shown in Figure 26. The first spectrum identified 
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as the Hamilton Standard SP15 Spectrum is obtained from [6], 

The other two spectra show harmonic levels falling off much 
more slowly as harmonic order increases; in one case the 
change in harmonic level being -3dB per harmonic and in the 
other case all the harmonic levels are equal. For present 
purposes it was assumed that the constant harmonic level 
specification was applied to the ten lowest order harmonics 
only. 

Pressure field coherence and trace (or phase) velocity were 
not covered by the Hamilton Standard prediction procedure [ 6 ], 
and other means had to be used to obtain the required informa- 
tion. One source of information was the data analysis [8, 9] 
performed on measurements from a test program on an Aero 
Commander airplane. The analysis included coherence and trace 
velocity results, although the test conditions (tip clearance 
approximately 0.05D, forward speed less than 40 m/s (120 ft/sec 
differed considerably from those associated with the study 
airplanes. 

The Aero Commander data [9] show that, where there is forward 
velocity of the airplane, the decay of the pressure field co- 
herence for the low order harmonics is very slow in both the 
longitudinal and circumferential directions. This is to be 
expected when the excitation consists of a series of discrete 
frequency components. Thus, a similar slow decay of press.ure 
field coherence was assumed for the study airplanes, the co- 
herence decay parameters c and c of Eq. (29) (c = (° x j c )) 

k / <y 

each being given a value of 0.01. 

The trace velocity in the longitudinal direction was estimated 
on the basis of acoustic wave propagation from an effective 


source location. The source 'was located at 10 % of the pro- 
peller radius, in the effective plane of rotation of the 
propeller, as shown in Figure 2 7. The effective plane of 
rotation was used to take into account the forward motion of 
the airplane. Figure 27 shows the construction of the trace 
velocity relationship for the inboard propeller of the wide- 
body airplane. The net trace velocity for the pressure field 
includes the influence of flight velocity so that the trace 
velocity is subsonic on the forward region of the fuselage 
and about twice the speed of sound on the aft region. Close 
to the plane of rotation of the propeller, the acoustic field 
is modeled as plane waves incident normally to the fuselage, 
so that the axial trace velocity tends to an infinite value. 

In the circumferential direction the trace velocity is obtained 
by assuming that the pressure field rotates with the propeller. 
Thus, for a constant tip rotational speed, the greater the tip 
clearance, the higher will be the circumferential trace velocity. 
The value of the trace velocity is supersonic for all conditions 
of Interest, and, because of the model chosen, the circumfer- 
ential velocity increases with distance from the propeller 
plane R, as is shown in Figure 25(c) for the inboard propeller 
of the wide-body airplane. The circumferential trace velocity 
in the region of the propeller plane (peak levels) is predicted 

from 

U y = 7rDf Q ( l+2r/D) 

where r/D is the tip clearance in fan diameters and f Q the 
blade passage frequency. Away from the propeller plane, Uy is 

calculated from 

U y = 2rrf q R 


where R is identified in Figure 27 for the 7th element. 
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Trace Velocity (m/sec) 




5.0 




MODEL OF STRUCTURE 

Calculation of the structural response to the propeller 
excitation pressure field requires the calculation or deter- 
mination of several functions such as joint acceptance, modal 
density, radiation ratio and structural damping. These func 
tions can, in turn, be determined only after a model has been 
constructed for the fuselage structure. As was discussed m 
Section 4, one of the basic assumptions about the baseline 
fuselage structure was that the construction would be the 
conventional skin-stringer-frame type used in current day 
aircraft. This section describes the steps taken to describe 
such a structure in a form suitable for computation in the 
analytical acoustic model presented in Section 3. The dis- 
cussion centers on the resonant structural response, since, 
as will be shown in Section 6, the power flow associated with 
resonant response is significantly greater than that due to 
non-resonant response. 

5,1 General Representation 

The analytical formulation in Section 3.2.2 gives the total 
power flow from the external excitation field through the 
sidewall structure into the interior, and, in the process, 
predicts a homogeneous, average level for the structural 
response. To allow for non-homogeneous spatial distributions 
of the external pressure, the structural response and the 
interior noise levels, it is necessary to divide the fuselage 
structure into several segments along the fuselage length. 
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Then the total resonant power flow equation for N structural 
elements can be written in the approximated form 


N 


W 


m 


U b ) 


■ 1 


27TAJP.C. 


n. 


1-1 


<P («- b ) > j 

XO D J 


(“h) 


n r 


< J r ( o b ) > P j < a rad 

( 74 ) 


Eq. (74) is obtained by substitution of Eq. (60) into Eq . (51) 
and by replacement of the modal mass Mr by the resulting value 
appropriate to sinusoidal structural modes, 

i.e. . M r = yA (75) 

“4 

where A is the structural area and’ y is the surface density. 
Eq. (74) applies to the frequency regime where structure 
and receiving cavity have sufficiently high values for the 
modal overlap. The summation is performed over structural 
elements j . 

The choice of the size and number of the structural elements is 
influenced by the following points. 

(i) The element dimensions should be large enough that the 
amplitude of the external pressure field decays by at 
least 10 dB within the element. 

(ii) The elements chosen should have structural wave number- 
frequency distributions characteristic of the detailed 
overall fuselage modes. 
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(iii) The total power flow into the Interior must not depend 
on the sizes of the elements chosen. 

The fuselage structure Is Idealized so that each structural 
element Is assumed to be part of a pressurized cylinder, I.e. 
a cylindrical plate of uniform thickness with uniformly spaced 
ring stiffeners (or frames) and stringers (Figure 28). It Is 
assumed that no power flows Into the interior from below the 
cabin floor so that the element is not a complete cylinder. 
Furthermore, when calculating resonance frequencies. It Is 
assumed that the structural element has clamped boundaries. 

The structural parameters used to describe the element will 
vary according to the wavelengths and frequencies of the modes 
of vibration [36]. 

At low frequencies, where the structural modal wavelengths, 
both axially and circumferentially, are greater than the frame 
and stringer separation distances respectively, both frames 
and stringers may be ’ smeared-out * o.rthotropically over the 
skin surface. The frequency equation, derived by Mikulas and 
McElman [25] and presented in Figure 5 of Section 3.8.1 (with 
the addition of a term due to pressurization), was found during 
study calculations, to give reasonable estimates for the reson- 
ance frequencies of a simply-supported ring stiffened cylinder [36], 
Mikulas* equation was then used for the resonance frequencies 
of an element of a stiffened' cylinder with edges clamped by 
using the following changes and approximations. 

( i) For axial mode order m, the simply— supported wavelength 
is replaced by the smaller clamped-clamped wavelength 
[37], and m is replaced by (m + 0.5)/1.05 
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A = Frame Pitch 
A = Stringer Pitch 
R = Radius of Curvature 


FIGURE 28. REPRESENTATION OF FUSELAGE STRUCTURAL 
ELEMENT 



For a complete cylinder, n is the number of circum- 
ferential full wavelengths. In this analysis, 
n Is the number of circumferential half wave- 
lengths, and 

n is replaced by j (n + 0.5)/l*05 

The ratio of circumferential to axial wavenumbers 
becomes 


k m R L y (m + 0.5) 

For modes of low order m in the axial direction and low order 
n in the circumferential direction, the resonance frequency 
equation is dominated by the membrane stiffness term due to 
curvature. Then, as n increases, the ring bending stiffness 
term becomes dominant. Mikulas' equation with frames and 
stringers represented as ' smeared*-out ' masses and stiffness, 
therefore covers regions I and II defined by Hu et al £36] 

The modal density is dependent on both bending and membrane 
stiffnesses and is not very sensitive to changes in frame 
bending stiffness alone. 

As frequency increases, the axial half-wavelength of the modes 
becomes less than the spacing between the frames, and the 
frame motion becomes small compared to the motion of the panel 
between the frames. Inclusion of the frame stiffness now 
results in overestimates of the panel resonance frequencies. 
Thus, in this mid-frequency region, the frame stiffness and 
mass are omitted from the analytical model, and only the 
stringers are ' smeared-out » over the surface. Mikulas' equa- 
tions are again used to calculate the modal frequencies . 


(ii) 


( iii) 


In the high frequency region, both the axial and circumfer- 
ential half-wavelengths are less than the spacing between 
frames and stringers respectively, and subpanel resonances 
occur. Only the skin stiffness and mass: are included now in 
the analytical model, as the panel motion is now effectively 
independent of the frames and stringers. 

The frequencies at which the structural parameters are changed 
are given by 

(i) The minimum value of the resonance frequency of 
the subelement, of size (A x L y ) , clamped along two 
adjacent frames, for any circumferential mode 
order n with half-wavelengths greater than the 
stringer pitch, using (stringer + skin) structural 
parameters, and 

(ii) The fundamental frequency of subpanel with clamped 
edges, whose size (£ x £ y ) ls given by (frame pitch 
x stringer pitch) , using skin structural para- 
meters only. 

The first of the above frequencies defines the boundary between 
low and intermediate frequency regimes in the structural model- 
ing, and the second frequency is the boundary between the 
intermediate and high frequency regimes. It should be noted 
that these panel sizes are used only in defining the frequency 
regimes for the structural modeling and are not used in calcu- 
lating structural response. In calculating panel response, 
panel dimensions L L are always used. 

x y 


5 . 2 Structural Represent ations, 

The factors to be considered in selecting the size of 
structural elements were given In Section 5.1. The discussion 
in this section Identifies those parameters which are Influenced 
by element area while Sections 5.3 to 5.5 Investigate these 
parameters In more detail with the objective of choosing a 
minimum element size while maintaining the accuracy of the 
response calculation. 

Empirical evidence suggests that, at low frequencies, the ^ 
element size should be much larger than the frame and stringer 
spacings. For example, vibration data for a Boeing 737 fuse- 
lage structure [38] exposed to boundary layer and jet noise 
excitation show that acceleration coherence can remain high 
over several adjacent panels (see Figure 29). This is true 
even though the coherence coefficient of the excitation pressure 
field decays much more rapidly [38] than it does for the prop- 
eller pressure field considered in the present study. 

In the circumferential direction, the coherence coefficient 
for propeller noise maintains a high value over large separa- 
tion distances, and the amplitude decays little with distance 
[Q], The structural element should, therefore, extend at 
least from the cabin floor to the ceiling centerline or, 
preferably , from floor to floor. For a ^asel no a - 

plane, the structural element would include 15 to 30 stringer . 
This idealization in the circumferential direction cannot give 
a spatial variation of sound level across the width of t e 
cabin, but gives, instead, a local average value of structure 
response and interior noise level. This is considered adequate 
because measurements at appropriate frequencies m presen 
day aircraft indicate that the lateral variation in ca in 
sound level is small (less than ±2dB on an octave band basis). 
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Longitudinal Separation Distance (m) 


598 Hz 


Circumferential Separation Distance (m) 


FIGURE 29. ACCELERATION COHERENCE MEASURED ON ADJACENT 
PANELS OF BOEING 737 FUSELAGE (M = 0.78, JET 
AND BOUNDARY LAYER NOISE) 






In the longitudinal direction, coherence is again high at low 
frequencies, but the pressure amplitude decays rapidly with 
distance. For example, the amplitude of the pressure field 
predicted using [6] decreases by 20 dB or more over a distance 
of 10 frame spacings on the larger baseline aircraft. 

For a non-homogeneous excitation, the acoustic power flow 
equation for N structural elements is given by Eq. ( 7 * 0 . 

Assume that the elements are identical, with equal values 
of generalized mass M rj and loss factor n r j • Furthermore for 
excitations involving strong spatial amplitude decay over the 
length of the structural element. 


c PXo (to b^ 


s PXo (u b )> J^l U 


1 , 2 , 3 - • •) 


where <p y 2 >, is the maximum pressure amplitude applied at 
the edge of ilement j . Therefore only the term involving 
element 1 will contribute significantly to the power inflow, 
if other factors such as excitation coherence, amplitude 
decay rate and trace velocity remain constant or vary only 
slightly. Hence the acoustic power flow equation may be 
written in the approximate form:- 


W± n U b ) ~ 


2lT Pj c l <P X o(“ b ) > j A .1 nr » 3 — b) < j? (a) b )> P j <arad> J ; J 1 

w b Pj' n r,j (76) 


The effect of element size then reduces to a consideration as 


to whether 


Ai - A. n r,j < ^r V. <0 rad j 

^ J 


for j = 1 


( 77 ) 


is independent of A. 
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The term <-p 2 > ,• _ is constant for different values of A., ’ 
A 0 J - i J 

since as area A. increases, the change is made such that the 
—• J 

location X Q of the maximum pressure always lies on one edge of 
the element j=l. For convenience, Eq. (77) may now be written 
in the alternative form 


<CT rad > j <J 'r > p.j) 


(78) 


The various terms on the right hand side of this equation 
will now be considered in more detail in Sections 5.3-5. 5. 


5*3 Modal Density 

For a flat plate, the modal density is given by 


Kp c l 

where K p is the radius of gyration of the panel. This equation 
will be approximately valid even for pressurized and stiffened 
panels, provided that the value of K p is determined for the 
equivalent ( smeared-out) panel. From Eq. (78), it is required 
that n rj j/Aj for fuselage panels of various Aj , but constant 
surface density and stiffener arrangement, should be independent 
of Aj , at least when the modal density is sufficiently high 
so that several modes exist in the frequency band of interest. 

Figures 30 and 31 present calculated values for n(oi)/Aj for 
typical narrow-body and wide-body fuselage elements. The curves 
in both figures, are associated with the low-frequency struc- 
tural model which includes frames, stringers and skin. For the 


n r ,j U b ) 
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VARIATION OF MODAL DENSITY WITH PANEL ELEMENT AREA • 
LOW FREQUENCY STRUCTURAL MODELS 


FIGURE 31. 






narrow-body fuselage (Figure 30), when L y 7.0m (i.e. element 
extends from floor to floor), the ratio n(co)/A Is indepen en 
element length L x when L, is greater than approximately 4m (8 
frames) for frequencies 100Hz and above. When L x ' 2 ' 5 
n(w) /A shows little change at frequencies above 200Hz, but 
lower frequencies the differences can be as much as ±5d . 

For the wide-body fuselage, an element of width L y = 5.22m 
(i.e. extending from floor to ceiling centerline) and a length 
L of at least 2.5m. (5 frames) is sufficient to give a value 
of n(w)/A which is independent of area for frequencies of 125Hz 
and above (Figure 31). An element of dimensions L x = 3.17m and 
L = 5.66m is satisfactory for the small-diameter fuselage 

y 

(Figure 31) • 


The modal density of the structural elements varies with 
frequency (as shown in Figures 30 and 31), but it also epen 
on the structural representations used in the analysis. 

Figure 32 shows the modal density for a typical wide-body 
fuselage element calculated using the three different struc- 
tural representations discussed in Section 5.1* 
that there is a large difference in modal density between the 
low-frequency structural representation of * smeared-out * 
frames and stringers and the high-frequency structural repre- 
sentation where only the skin is considered. Differences of 
this type have been observed in practice [36] . 
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Modal Density/Unit Area, -i— J Modes rad” 



Typical Wide Body Fuselage Element 
Size: L x = 5.08 m (10 Frames) 

L = 5.22 m (28 Stringers) 

y 

FIGURE 32. VARIATION OF MODAL DENSITY WITH FREQUENCY 
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5.4 


Joint Acceptances 


The effect of element size on the product Aj <j r > p is shown 
ln Figure 33 for typical narrow-body and in Figure 34 for wide- 
body and small diameter fuselage elements, the calculations 
having been performed for the low frequency structural models, 
which include frames, stringers and skin. For the narrow-body 
element lengths L x greater than 4m, Aj 2 <j 2 > p is independent 
of A.. When the elements are much smaller than 4m, the approxi 
mate 0 ' equation (Eq. (76) ) ' for the power flow is no longer valid 
and the contribution from adjacent elements cannot be ignored. 
More importantly. A 2 <j 2 > p is significantly underestimated 
in these situations due to the preclusion of long wavelength, 
well-coupled structural modes. 


Figure 34 shows that A 2 <J r 2 > p is independent of Aj , for the 
wide-body fuselage, for an element of width L y = 5- 22m and 
length L greater than 5.08m, and for the small diameter fuse- 
lage fo/an element L = 5.66m and length L x greater than 3.17m 

These conclusions regarding the effect of element size on the 
progressive wave joint acceptance are valid only for a strong 
decaying excitation pressure amplitude. If the excitation, 
is a homogeneous, reverberant acoustic field* the progressive 
wave joint acceptance is replaced by the reverberant field 
joint acceptance <j 2 > R . From Eq. (61), 


A. 


<J r >R 


77 c i °rad /2a) 


( 61 ) 


Now, as will be shown in Section 5.5, the structural element 

can be chosen such that a rad is independent of Ay Then 

A <i 2 >-n will also be independent of A.. This is shown m 
j J r R J 

Figure 3:5 • 

* required for acoustic power flow calculations (See Section 


3.6.3) 
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Band-Averaged Progressive Wave Joint Acceptance A \j 
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FIGURE 33. 


EFFECT OF PANEL AREA ON BAND-AVERAGED 
PROGRESSIVE WAVE JOINT ACCEPTANCE 
(TYPICAL NARROW BODY FUSELAGE ELEMENT, 
LOW FREQUENCY STRUCTURAL MODEL) 
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FIGURE 34. EFFECT OF PANEL AREA ON BAND-AVERAGED PROGRESSIVE WAVE 
JOINT ACCEPTANCE: LOW FREQUENCY STRUCTURAL MODELS 
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Band-Averaged Reverberant Field Joint Acceptance A <^j r 


L 

L 

A 

X 

X 

. 2. 

(m) 

(m) 

(m ) 

2.54 

2.29 

5.8 

2.54 

7.0 

17.8 

5.08 

7.0 

35.6 

10.2 

7.0 

71.4 


A / *2\ _ 7r c ? rad c /,,* 

'Jr 'R ~ 2 ,,2 Ec l- ( 61 ) 



Frequency (Hz) 

(Note: Except at Low Frequencies, All Curves 
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FIGURE 35. EFFECT OF PANEL AREA ON BAND-AVERAGE 
REVERBERANT FIELD JOINT ACCEPTANCE 
(TYPICAL NARROW BODY FUSELAGE ELEMENT 
LOW FREQUENCY STRUCTURAL MODEL) 





5.5 


Radiation Efficiency 


At low frequencies, with structural elements of sizes pre- 
viously discussed (L x “ 2,5 m, L - 5m) so that they include 

y 

several frames and stringers, the radiation efficiency o ra ^ 
is unity for all modes, and hence is independent of the 
element area Aj . This occurs because at low frequencies the 
modes are acoustically fast (and hence efficient radiators) 
due to the effects of fuselage curvature, the ' smeared-out ' 
stiffness of the frames and stringers and, to a lesser extent, 
cabin pressurization. Modal wavelengths at low frequencies 
are large compared to typical frame and stringer pitch, and 
the corresponding mode shapes are well represented by simple 
sinusoidal functions, even though the element boundaries are 
assumed clamped for resonance frequency calculation. The effect 
of stiffeners is to modulate the. basic mode shape with a wave- 
length equal to twice the stiffener pitch. As frequency in- 
creases however, the stiffeners participate less in the fuselage 
motion. Then, mode shapes are determined more by subpanel prop- 
erties with stiffeners acting as restraining edges. The radia- 
tion properties of such 'ribbed' structures may be estimated 
following Maidanik [39] or Barger et al [ 4 0 ] : the latter ap- 
proach was adopted, from which it may be shown that, for an 
isotropic free vibration field 


a . _ 8 
rad - - 


L_ 

A 


1 + (1 + 


jm_) 2 ' 


Apk 


P J 


( 80 ) 


where L is the sum of the total length of rib attachment to the 
skin surface and one-half the total length of ribs attached to 
the area edge. A is the plate area, A c = 2ir kc^/c^ is the 
skin critical wavelength, k is the skin radius of gyration, 
h is the skin thickness, A r is the cross-sectional area of 
the stiffeners, k p is the structural wavenumber at frequency 
co, and e = sin~ (k c ^/kp). Since L/A is constant for a uni- 
form grid of arbitrary area the radiation ratio is independent 
of the area. 


- 134 - 


5.6 


Loss Factors 


For lightly damped aluminum structures, only structural damp- 
ing and acoustic radiation will contribute significantly to 
the decay of vibrational energy levels. Structural energy 
dissipation will be associated both with material (internal) 
damping mechanisms, and with energy losses occurring at the 
structural boundaries; of the latter mechanisms, viscous 
damping arising from the pumping of gas within the joints can 
be important in non-rigid joints such as those occurring in 
bolted built-up structures [41]. Acoustic radiation damping 
can become a significant means of energy dissipation for 
lightweight structures whose radiation ratios are high. 

The overall modal loss factor n r may be expressed by the 
relation 


’V + n j + n rad 11 s true + 11 rad (8l) 

where is the internal loss factor associated primarily with 
material damping (which for aluminum can be considered to be 
essentially independent of frequency) , nj is the loss factor 
associated with dissipation at structural joints and n rdd is 
the modal loss factor resulting from acoustic radiation 
damping. n rad is given by 

"rad = R rad /u,Ap = pC ' , rad / “ ,J (82) 

where w is the circular frequency, y and A are the surface den- 
sity and area of the structure, R rdd ls the radiation resistance, 
p and c are respectively the density of, and speed of sound in, 
the fluid surrounding the structure and c rad is the radiation 
efficiency or radiation ratio. 
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For lightly damped structures, nj_<10“ 3 , while for rigid Joints 
allowing no relative motion between the component parts, 
n j <ri i* For lightweight structures and supersonic modes, 
o - 1, so that n rad can be greater than both and nj. This 
can be particularly important for cylindrical structures as 
they have supersonic modes at frequencies well below the 
acoustic coincidence frequency of a flat plate with the same 
thickness as the skin. The same statements are assumed valid 
for band-averaged loss factors, n r , which represent the arith- 
metic mean of the modal loss factors in the frequency band of 
interest. The structural component n g ^ ruc of the loss factor 
for an undamped bare structure, with no applied damping treat- 
ments, is assumed to be 0.01 at 100 Hz and to vary inversely 
with frequency [4], i.e. 

n , „ = 0.01 t2.7r.lOO/&>) (83) 


When the radiation loss factor (Eq. (82)) is included, it can be 
seen (see Figure 36) that the structural loss factor for an un- 
damped structure is small compared with the radiation loss 
factor at most frequencies of current interest. Limited experi- 
mental data is available [42] to support the calculated values 
as seen in Figure 36. Also shown are the total loss factors 
for a highly-damped structure, assuming that n gtruc has a con- 
stant value of 0.05, a value considered to approach the maximum 
achievable value in practical fuselage structures at low fre- 
quencies using proven techniques [42], In this case the radia- 
tion loss factor is no longer the most significant contribution 
to the total loss factor at low frequencies. 


- 136 - 



Loss Factor 










5. 7 Structural Idealization for Computation Purposes 

The preceding discussion has described how the fuselage struc- 
ture is represented as a series of structural elements (see, 
for example. Figure 28). The elements extend from floor-to- 
floor or floor-to-ceiling in the circumferential direction, 
and divide the fuselage into several segments in the longitudi- 
nal direction. The selection of these elements for the three 
study airplanes is discussed in this section. In addition, 
the dimensions and structural characteristics of the elements, 
and the excitation characteristics associated with each element 
are presented in tabular form for reference. 

The division of the fuselage structure into a series of elements 
was determined to some extent by the amplitude decay of the 
excitation pressure field for each study airplane in both the 
longitudinal and circumferential direction. Element sizes 
were selected for each airplane on the basis of sensitivity 
calculations similar to those presented in Section 5-1 - 5.5 
for the wide-body airplane. 

Location of the boundaries between adjacent segments was 
selected in the first place on the basis of the location of 
peak sound level from the inboard propeller. For computation 
reasons it was desired to have this peak level occur on the 
boundary between two structure elements, as is indicated in 
Figures 37 through 39 which show schematic sideviews of the 
structural segmentation. Thus, for example, the peak sound 
level associated with the inboard propeller of the narrow-body 
airplane (Figure 38) occurs on the boundary between structure 
elements 3 and 4. Similarly, a second boundary between ele- 
ments was selected to correspond with the peak sound levels 
associated with the outboard propeller (Figures 37 and 38). 
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Location of Peak Excitation Pressure 
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FIGURE 38. 


SKETCH SHOWING 
AND PROPELLERS: 


LOCATIONS OF STRUCTURAL ELEMENTS 
NARROW BODY AIRPLANE 
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FIGURE 39. SKETCH SHOWING LOCATION OF STRUCTURAL 
ELEMENTS AND PROPELLER: SMALL DIAMETER 

AIRPLANE 
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Having defined these boundaries, and, having determined 
acceptable element sizes from parametric analyses in Sections 
_ 5 . 5 9 the segmentation of the fuselage could be completed 
The resulting divisions for the three airplanes are shown in 
Figures 37 through 39 and associated segment locations, dimen- 
sions and excitation characteristics are given in Tables 4 
through 6. For convenience, element locations are given in 
terms of the distance of the upstream edge of the element from 
the forward pressure bulkhead of the cabin. 

The selection of the location of peak sound, level as the factor 
determining element boundaries allows the excitation field 
for a given element to be described. in terms of a monotoni- 
cally-decaying pressure amplitude with a mean trace velocity 
either upstream or downstream. This. representation reduces 
the complexity of the analytical model. 

Structural data for the three baseline study airplanes, as 
required for the calculation of acoustic power flow into the 
fuselage interior, are presented in Table 7. These data are 
also applicable to the baseline airplanes with add-on noise 

control . 
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Table 4 


Structural Segmentation Details: Wide-body 


Element 



Element No. 



Descriptors 

1 

2 

3 

4 

5 

6 

7 

Location of up 
stream edge of ele- 
ment, re forward 
bulkhead (m) 

o„o 

7.11 

14.21 

L9.29 

24.37 

29.45 

34.53 

Element dimensions 
L (m) 

X 

L (m) 

y 

7.105 

7.105 

5.08 

5.08 

5.08 

5.08 

6 . 92 

5.227 

5.227 

5.227 

5.227 

5.227 

5.227 

5.227 

Location of maxi- 
mum excitation 
level re forward 
bulkhead (m) 



19 . 

29 




Amplitude Decay 
Rates 

a L 

X X 

5.77 

5.87 

3.77 

3.56 

3.50 

3.50 

4.77 

a L 

y y 

1.15 

1,15 

1.15 

1.15 

1,15 

1.15 

1.15 

Trace Velocities 
(m/s) 

U v 

68 . 

83 . 

260 . 

745 . 

574 . 

556 . 

551 . 

X 

u 

V 

1390 . 

851 . 

498 . 

498 . 

780 . 

1153. 

1626 . 

Coherence Decays 
c =c* k 

X XX 

c =c' k 

y y y 

0.01 

0.01 

0.01 

0.01 

0 ,01 

0.01 

0.01 

0.01 

0.01 

0.01 

0,01 

0.01 

0.01 

0.01 


- 143 - 




Table 5 


Structural Segmentation Details: Narrow-body 


Element 

Element No. 

Descriptors 

1 

2 ' 

3 

4 

5 

6 

Location of upstream 
edge of element, re 
forward bulkhead (m) 

0.0 

5.08 

10.16 

15.24 

20.32 

25.^0 

Element dimensions 

L x (m) 

L y (m> 

5.08 

7.00 

5.08 

7.00 

5.08 

7.00 

5.08 

7.00 

5.08 

7.00 

8.60 

7.00 

Location of maximum 
excitation level re 
forward bulkhead (m) 

15 

.24 



Amplitude Decay Rates 
a L 

X X 

a L 

y y 

3.12 

0.70 

3.12 

0.70 

2.51 

0.70 

1.98 

0.70 

1.98 

0.70 

3.35 

0.70 

Trace Velocities (m/s) 
U x 

u y 

95. 

1703. 

148. 

1241. 

553. 

926. 

1038. 

926. 

633. 

L24l. 

572. 

1880. 

Coherence Decays 
c =c ' k 

X XX 

c =c 1 k 

y y y 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 

» 

0.01 

0.01 
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Tabl e 6 


Structural Segmentation Details: Small diameter 


Element 

Element No. 

Descriptors 

1 

2 

3 

Location of upstream edge 
of element, re forward 
bulkhead (m) 

0.0 

3.73 

6.90 

Element dimensions 




L x (m) 
L y (m) 

3.73 

3.165 

3.165 

5 . 66 

5.66 

5 . 66 

Location of maximum ex- 
citation level re forward 
bulkhead (m) 

3.73 

Amplitude Decay Rates 




a L 

XX 

5.71 

4.50 

4.71 

i a L 

; y y 

1.79 

1.79 

1.79 

Trace Velocities (m/s) 




u y 

145 . 

473 . 

658 . 

451 . 

560 . 

791 . 

Coherence Decays 




c =c' k 

X XX 

0.01 

0.01 

0.01 

c =c' k 

y y y 

0.01 

0.01 

0.01 
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i Table 7 

Structural Data for Typical Element of 
Baseline Airplanes 


Baseline Aircraft 

Element No. (in propeller plane) 
Frame Spacing % m 
Stringer Spacing m 

Skin 

Thickness h mm 

Mass/unit area ph kg/m 

Bending stiffness, D > D X > D V > N,m 

Torsional stiffness, D, , N.m 

Wave speed, c L m/sec 

Stringer (smeared out over 
width % ) 

| Equivalent thickness, A g /&y 
Mass/unit area 

p A / 

[Bending stiffness ,E g X s /Jt 
Torsional stiffness, 

G J /* N - m 

, , ^ . s s y 

Centroid to skin 

middle surface z g 
[Bending stiffness, D x N.m 

1 (skin & stringers) 

Frame ( smeared out over width SL^) 


mm 

kg/m 2 
N .m 


Equivalent thickness A^/Jt 
Mass/Unit area pA /Jl x 
Bending stiffness, 

E 1 /% 

Torsional stiffness 

G J /a 

I 2? T* X 

|centroid to skin 

middle surface 

|Bending stiffness D 
(skin & frames) y 


m 

kg/m 2 


N.m 


Stiffener Cross-sectional Area m 2 

rib = <Vx + Vy )/( h + *y> 


A 


Wide 

Narrow 

Small 

4 

. 5080 
.1829 

4 

. 5080 
. 2286 

4 

. 3302 
.1270 

1.7780 

4.8006 

38.133 

12.584 

5491.1 ! 

3.0 

8.100 

183.18 

60.45 

5491.1 

1.016 

2.7432 

7.H51 

2.3480 

5491.1 

.9882 

.4006 

.4064 

2.6680 

6510.05 

1.0817 

4338.9 

1.0973 

687.23 

19.388 

2.3164 

3.8956 


-.01583 

18,5^9 


■1.0448 

2.8209 

174,872 

24.510 

-.05468 

323,420 

2734x10”. 3 


-.01855 

13,410 


.6012 

1.6233 

53,956 

8.1609 

-.03921 

110,780 


-.00744 

1,895.2 


.3716 

1.0034 

20,213 

1.2314 

- . 03911 
51,237 


.1579x10” 3 .7136xl0” 4 
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6.0 NOISE LEVELS IN BASELINE AIRPLANES 

As a first step in the analysis of the interior noise levels 
in the three study airplanes, interior noise levels were calcu^- 
lated for the three airplanes with bare interiors and with 
baseline sidewall treatment. The results of the calculations 
are presented in this section for the three different excita- 
tion spectral shapes described in Section 4.6. Subsequent 
sections will present the results of applying the analytical 
model to add-on and advanced noise control concepts . The 
information presented in this section shows axial distribu- 
tions of fuselage acoustic power in-flow, axial distributions 
of interior overall noise levels, resonant and non-resonant 
contributions to the total acoustic power flow, and low-order 
harmonic levels. Finally the required noise reductions for 
the baseline study airplane to achieve the criterion of 80 dB 
maximum interior A-wei.ght.ed -noise level are presented. 

6.1 Hamilton Standard Excitation Spectrum 

The major part of the discussion in this section will be 
concerned with analytical results obtained using the excita- 
tion spectrum defined in the Hamilton Standard prediction 
procedure [6], since this spectrum is associated with the 
highest predicted interior sound levels. 

6.1.1 Acoustic Power Flow into Fuselage 

The calculated acoustic power flow into the baseline airplane 
interior for the two lowest order harmonics and for each study 
airplane is shown in Figure '40: here the airplane fuselages 
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are each lined with the baseline sidewall treatments. It is 
clear that, while maximum transmission occurs in the region 
of maximum excitation level, the power flow is distributed 
over a considerable portion of the fuselage length. This 
is due, for the wide-body and narrow-body fuselages, to the 
influence of the outboard propellers whose directivities have 
less pronounced effects than do the inboard propellers. The 
contributions of the inboard and outboard propellers are in- 
dicated separately. In general it is seen that fuselage 
sections extending from just forward of the inboard propeller 
plane of maximum excitation level to well-rear of the outboard 
propeller plane of maximum excitation level will require 
acoustic treatment to achieve significant reductions in 
acoustic power flow. 

Further it will not be possible to confine the acoustic power 
flow to a small element of the enclosed cabin by partitioning 
the cabin space. This conclusion follows from the fact that 
the calculated power flow distribution is a direct indication 
of the noise transmission distribution through the sidewall 
irrespective of any subsequent distribution of energy within 
the cabin. 

Spectra of the net calculated power flow through the baseline 
fuselage of each of the study airplanes (Eqs. (51) and (52)) 
are shown in Figure 4l. Also shown in the figure are the contri 
buttons associated with resonant and non-resonant transmission 
through the fuselage structure. At higher frequencies where the 
effects of frames are not included in the structural representa- 
tions, the structural resonant modes become less efficient 
radiators so that the resonant power approaches that of the 
mass-controlled non-resonant modes. The calculations clearly 
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indicate however, that over the frequency range of interest, 
resonant transmission controls the acoustic power flow. 

6.1.2 Interior Noise Levels 

Having calculated the power inflow, through each of the fuse- 
lage sections, it is now possible, by means of the imaging 
procedure outlined in Section 3*10 to calculate the variation 
of sound level along the length of the cabin. 

Prior to considering each particular study airplane, the 
effect of variations in interior absorption coefficient on 
the axial variation in interior sound levels is examined. 

Figure 42 shows, for the narrow-body airplane, the axial 
sound level variation for two different values of absorption 
coefficient for the case where a single panel element (number 4) 
is radiating into the cabin interior. a= 0.05 corresponds 
essentially to a bare fuselage structure while a = 0.75 would 
be considered the maximum achievable value in the study air- 
planes. The smaller the value of a, the more diffuse the 
internal sound field. Some unpublished empirical data is 
available to support these calculations. For example. Figure 43 
presents a comparison of calculated and measured variations in 
cabin sound levels for a small diameter propeller-driven air- 
craft: good agreement is found between measured data and 
calculations using the image source array model. 

Figure 44 shows the calculated axial variation of cabin A-weighted 
interior noise levels for the study airplanes, where the baseline 
fuselage structures have been used. Results for the bare fuse- 
lage and baseline add-on sidewall treatment for the Hamilton 
Standard excitation spectra [6] are presented. Maximum noise 
levels of the treated baseline fuselage are between 5 and 10 dB 


Sound Pressure Level 


FIGURE 



42. AXIAL VARIATION IN SOUND LEVEL FOR 
DIFFERENT INTERIOR ABSORPTION 
COEFFICIENTS (NARROW BODY) 


Change in Noise 
Level dB 


Mach No. = 0.49 
Altitude = 6 , 100 m 


O 53 Hz 
A 106 Hz 
V 212 Hz 
□ 450 Hz 



0 1 2 3 4 5 6 

1 l l l - I _J I 
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Data Points are the Average for 5 Microphone Positions at Each Longitudinal Location. 


FIGURE 43„ COMPARISON OF CALCULATED AND MEASURED 
VARIATION OF SOUND LEVEL ALONG CABIN 
OF A SMALL DIAMETER PROPELLER-DRIVEN 
AIRCRAFT (Data From BBN Files) 


-Weighted iouna Level 


WIDE BODY 


Bare Fuselage 
Baseline Interior 



Design Criteria 


Distance Along Cabin (m) 


NARROW BODY 


Design Criteria 


10 20 30 

Distance Along Cabin (m) 



_ 

— 




^ ««* 

\ 

\ 





SMAL 

r 

L DIAMETER 


FIGURE 44. 


jbMAL L uimac* | Desig - Criteria 

0 2 4 6 8 10 

Distance Along Cabin (m) 
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lower than for the bare fuselage. That the effect of the 
baseline sidewall and associated Increased Interior absorption 
Is so small results because the A-weighted Interior noise 
levels are In general controlled by the excitation level at 
the propeller blade passage frequency, where the transmission 
loss provided by the baseline sidewall Is small, particularly 
for the larger diameter fuselage. The axial variation for the 
bare fuselage Is notably less than that for the baseline 
treated fuselage, a result of the decreased absorption present 
In the bare fuselage. 

It Is clear that conventional airplane sidewall and fuselage 
designs are Inadequate to achieve the design criteria (i.e., 

80 dB(A) maximum Internal noise level) In the study airplanes. 
Therefore substantial Improvement will be required In the 
fuselage/sidewall noise reduction. 

Cabin noise spectra at the location of maximum calculated dB(A) 

i 

levels for each study airplane for bare and baseline fuselages, 
and for the Hamilton Standard excitation spectra [6] are pre- 
sented In Figure M5. For each airplane, the propeller funda- 
mental dominates the Internal A-welghted noise level, as well 
as the Interior noise spectra. Baseline sidewall treatments 
provide minimal noise reduction for these low-frequency ex- 
citations, but are increasingly effective as harmonic order 
Increases . 

Table 8 presents a summary of the minimum noise reductions 
required for the baseline structures to achieve the design 
goal of 80 dB(A) : . these minima have been assessed at the 
propeller fundamental frequencies since It is at these low 
frequencies that it is most difficult to achieve high noise 
reductions. 


Sound Pressure Level, dB re 20^ Pa 



CABIN SPECTRA AT LOCATION OF 
FOR BARE AND BASELINE TREATED 


MAXIMUM dBA 
FUSELAG ES 


LEVELS 


FIGURE 45. 
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1 


1 


TABLE 8. 

MINIMUM ADDITIONAL NOISE REDUCTION REQUIRED FOR BASELINE STRUCTURE 


Baseline Structure 

W i d e - B o dy 

Narrow-Body 

Small Diameter 

Maximum Cabin dBA level 

105.5 

104.9 

108.1 

Frequency of Fundamental, Hz 

133 

183 

231 

Sound pressure level at fundamental 

120.8 

116.5 

117 • 4 

frequency, dB 

95.5 

91.8 


Sound pressure level at fundamental 

89.5 

frequency to meet 80 dBA goal, dB 




Minimum Noise Reduction required 
at Fundamental Frequency, dB 

25.3 

24.7 

27-9 










Limited experimental evidence from current commercial aviation 
airplanes exists to support the above conclusions. In Figure 
46 are shown comparisons between calculated baseline noise 
reductions for each of the study airplanes with empirical data 
for narrow-body fuselages, where noise reduction is taken to. 
be the difference between peak excitation levels and the maxi- 
mum internal noise levels. The current predictions agree well 
with the range of data, lending some confidence to the analyti 

cal predictions. 

6.1.3 Error Analysi s 

For the study airplanes of interest. Table 9 and Figure ^47 
show the expected error limits involved in the computation; 
these estimates result only from modeling of the fuselage 
structures and the coupling with the cabin volume, and do not 
include such errors as may result from inaccuracies in the 
description of the excitation pressure field. 

Following the approach outlined in Section 3.5*3, the normal 
ized standard error in calculations of the net acoustic power 

flow' can be derived as 

e r = [ 2lT ^r n r (w b )u) b]” / + [ 27TT1 n n n ( “b } “b ] 

The upper and lower bound 2a limits for 95* confidence in 
power flow calculations are calculated as 10 logio (l+2e r ) 
and 10 logic (l-2e r ) above and below the expected value of 
the power flow. These error limits should be interpreted as 
follows: at any time, the actual value of the acoustic 

power flow is expected to fall within the noted limits 95* 
of the time such a check is made. 
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FIGURE 46. COMPARISON OF CALCULATED BASELINE NOISE 

REDUCTIONS WITH EMPIRICAL DATA FOR NARROW 
BODIED FUSELAGES 
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FIGURE 47. 
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POWER FLOW INTO THE 
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AIRPLANE CABINS 


6 


Table 9. 


Calculated error limits for the 
acoustic power flow into airplane cabins 



Upper Limit 
(10 log(l+2E))*dB 

Lower Limit 
(10 log (l-2e))dB 

Harmonic Order 

1 

2 

— 

3 

1 

2 

3 

Wide-Body 

+ 2.7 

+ 2.1 

i — 1 

• 

i — 1 
+ 

1 

CO 

• 

vn 

-4.2 

-1.5 

Narrow-Body 

+2.3 

+ 0.9 

+0.9 

-4.9 

-1.1 

-1.2 

Small-Diameter 

+ 2.3 

CO 

• 

o . 

+ 

+0 . 6 

C\J 

LPv 

1 

O 

i — 1 

1 

1 

O 

~<3 


*2e corresponds to 2a limits on the acoustic power flow. 


6. 2 Alternative Excitation Spectra 

The calculated axial variation in cabin A-weighted interior 
noise levels for the study airplanes, with baseline fuselage 
and sidewall structures, but with alternative excitation 
spectra, are presented in Figure 48, and the corresponding 
cabin noise spectra at the location of maximum dB(A) level' are 
shown in Figure 49. . The differences in maximum dB(A) noise 
level between different spectra for each of the airplanes 
(Figure 48) result from changes in the level of the' funda- 
mental relative to the SP15 baseline excitation, as shown 
directly in Figure 49. Errors in estimation of the funda- 
mental of the excitation spectrum are directly reflected as 
errors in the A-weighted interior noise levels. 

Since the Hamilton Standard spectrum represents the state-of- 
the-art prediction and also leads to the most conservative 
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CABIN SPECTRA AT LOCATION OF MAXIMUM smiKin 
LEVEL FOR DIFFERENT EXCITATION SPECTRA S0UND 
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(i.e., highest) interior noise esWma * > ctrum model for 
calculations will be performed usxng this spectrum 

the excitation. 


6.3 Concl usi ons 

Clearly, substantial increases in the ^“^^^“^'/"^“study 
fuselage/sidewall structure are -"or^ ^ Qf 80 dB 
airplanes to achieve the g The requlred 

maximum anywhere within t re<au ction vary between 

rrr. « irr — :^i :::»ation r - 

^ than the noise 

. . „n fncpiaee structures, typicax 
"s! 1 “many for the larger d— aircraft which 

w l areer diameter propellers, 

are driven by slower, large* 
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7. ADD-ON NOISE REDUCTION METHODS 


7 . 1 Genera 1 

The range of add-on noise control treatments available for con- 
sideration in the analytical model includes variations of the 
sidewall transmission coefficient, variations to the internal 
absorption coefficient, increased damping and stiffness of the 
primary structure by add-on means, and optimizing the location 
of interior furnishings such as galleys and lavatories. Following a 
series of preliminary calculations of the structural character- 
istics and of the baseline interior noise levels, it is apparent 
that certain approaches would not provide significant increases 
in the fuselage noise reduction. For example, at low frequency, 
the "smeared" properties of the fuselage structure are dominated 
by the major stiffeners (frames, stringers) so that the adhesion 
of honeycomb panels to subpanel elements would produce little 
change in "smeared" stiffness and mass characteristics of the 
structure. Consequently, the associated changes in fuselage 
noise reduction would be small. Partitioning of the cabin volume 
and, similarly, placement of galleys, lavatories, etc. at 
particular fuselage locations close to the areas of peak excita- 
tion pressure levels, are also of little benefit for noise control 
because of the. large extent of the fuselage structural area aft of 
forward propeller plane that transmits significant levels of 
acoustic power into the cabin interior. 

Thus, efforts have been concentrated on those add-on methods of 
noise control for which it is clear that substantial increases in 
fuselage noise reduction would be possible. These include 
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1. optimizing the low-frequency absorption coefficients; 

2. increasing the structural loss factor of the fuselage at 

low frequencies -by the addition of patches of damping material 
to the skin and stiffeners; and 

3. upgrading the noise reduction . of add-on sidewall treatment 
by construction of a limp— mass double wall system. 

The results of computations for these add-on treatments are now 
presented. 

7 . 2 Maximizing the Interior Absorption Coefficients 

Baseline absorption coefficients for the various study airframes 
are discussed in Section 4 . 5 . At low frequencies, the absorption 
coefficients are generally less than 0.3 so that there is some 
potential for increasing the absorption at these frequencies. The 
improvements would be practical only on sidewall and ceiling sur- 
faces where significant cavity depth exists, although the use of 
relatively massive trim panels in a double-wall configuration will 
tend to limit the increases in sidewall absorption coefficients. 

For modeling purposes, however, changes in interior absorption 
are assumed to occur only on the sidewall surfaces. The maximum 
absorption considered practical is shown in Figure 50 together with 
baseline data. The effects of this change in cabin absorption on 
cabin A— weighted sound levels are shown in Figure 51 for each of 
the study aircraft. At the location of maximum internal sound 
level, reductions of only 1 to 2 dB are produced, while away 
from this location the benefits of increased sidewall absorption 
are somewhat greater with changes in space-average A-weighted 
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FIGURE 51. EFFECT OF CABIN ABSORPTION ON CABIN 
A"W El G HT E D SOUND LEVELS 
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The reason that 


noise levels varying between 2 and 4 dB. 
the peak levels change by such small amounts is due to the 
dominance of the direct contribution of sidewall radiation. 

The reverberant level contributions to local A-weighted sound 
levels lie some 3 to 4 dB below the direct field levels. 1 

It is concluded that maximizing interior absorption coeffi- 
cients will lead to only small decreases in interior noise 
levels, insufficient to achieve the 80 dB A-weighted interior 
sound level criterion. 

7.3 Increased Structural Damping 

It has been shown, in Section 5.6 and, in particular. Figure 
36, that the low frequency structural total loss factors for the 
baseline aircraft are mainly determined by radiation losses. 

The application of damping treatments to the fuselage struc- 
ture must therefore be quite extensive to be at all effective 
in increasing the total low frequency structural loss factors. 

In practice this is not an easily accomplished task, particu- 
larly since at low frequencies the add-on treatment must be 
applied to stiffening members rather than to the skin alone , 

When damping treatments are applied to the fuselage skin, 
increased structural loss factors and increased fuselage noise 
reduction are produced only at frequencies of the order of, 
and above, the fundamental resonance frequencies of the sub- 
panel elements (bounded by frames and stringers), as shown 
in Figure 52 taken from [35] . 

When damping treatments are applied to stringers in stringer- 
panel structures, using constrained layer techniques, the 
structure component, t s .j- ruc , °f the total loss factor can be 
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increased, by conventional techniques, to values in the range 
0.05 to 0.10 [42]. If similar treatments are applied to the 
fuselage frames as well, it is probable that values of n 
of 0.05 will be possible at low frequencies, where the 
frame stiffness influences the fuselage response and acoustic 
power flow. Taking this value n struc =0.05 for all frequen- 
cies, the total loss factor for the fuselage structure has 
been calculated for the three study airplanes and the result- 
ing values are given by the upper curves in Figure 36. 

These values of total structural loss factor were then used 
to estimate the probable maximum effect of add-on damping 
treatment on interior sound levels for each of the study 
airplanes. A typical spectral result is shown in Figure 53 
for the narrow-body airplane, where results are compared for 
the baseline airplane and the airplane with anticipated maxi- 
mum structural loss factor. Results showing the effect of 
increased structural damping on the maximum A-weighted levels 
for each of the study airplanes are also shown in Figure 53. 

It is clear that, at best, add-on damping treatments will 
provide only small noise reductions (3 'to 5 dB) relative to 
baseline levels, and that add-on damping treatments of them- 
selves will not provide the noise reduction required to 
achieve the design interior sound level criterion. Further, 
the combination of add-on damping treatments and tuned in- 
terior absorption, whose effects are essentially additive, will 
achieve interior sound levels only 5 to 7 dB lower than those 
of the baseline airplanes. 

7.4 Double-Wall Sidewall Configurations 

A preliminary investigation of add-on sidewall treatments by use 
of the analytical model indicated that large decreases in cabin 
sound level can be achieved by means of double-wall structures 
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which have fundamental resonance frequencies well below the 
fundamental excitation frequency. The lowering of the double- 
wall resonance frequency can be accomplished either by increasing 
the cavity depth & or by increasing the, trim panel surface 
density p 2 (see Section 3.9 for nomenclature definitions). The 
presence of porous material within the cavity of the double wall 
is important to damp both the double-wall resonance and the 
airgap resonances which occur at frequencies where the airgap 
depth is approximately equal to an acoustic half-wavelength. On 
the other hand, the use of a single layer of glass fibers, 
either alone or in combination with light-weight septa were 
shown in C 53 to provide insignificant increases in sidewall trans- 
mission loss at the low frequencies characteristic of the prop- 
fan fundamental blade passage frequencies. 

Sample calculations, using the model presented in Section 3.9, 
have been performed to show the estimated transmission losses 
associated with different double wall systems where a limp 
trim panel was added to each of the baseline fuselage struc- 
tures. The calculations were performed both for an empty air 
space between the two walls and for a space filled with porous 
material. The depth of the airgap was taken to be 127mm (5 in), 
a value which is typical of the maximum depth acceptable in 
current day commercial aircraft. Resulting transmission loss 
spectra for the wide-body fuselage are compared in Figure 54 
for two values of the trim panel surface density y 2 (2.0 and 
8.0 kg/m , or 0.4l and 1.64 lb/ft 2 ) and with porous material 
having characteristics similar to those of PF105AA glass fibers 
and a bulk density of 9.6 kg/m 3 . The results for other study 
airplane fuselage structures are essentially identical with those 
presented. The curves in Figure 54 show that the effect of increas- 
ing y 2 from 2.0 to 8.0 kg/m 2 is to increase the sidewall transmission 
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Increase in Sidewall Transmission Loss - -10 log (r $ ) d& 








loss by about 12 dB, or by 6 dB per doubling of y 2 . The 
double-wall resonance frequency decreases in accordance with 
the frequency equation given in Figure 5 4. At frequencies 
just above the double-wall resonance frequency, the effect of 
the porous material in the air space is small, at least for 
moderate values of flow resistance R 1# As Rj increases, the 
analysis predicts substantial increases in sidewall trans- 
mission loss, even at the double-wall resonance frequency. 
However, it is possible that at these high values of Rj the 
assumption that there is no structureborne transmission 
through the material fibers is no longer valid. In all cases 
the calculations assume that the trim panel is highly damped, 
or limp, and that there is no mechanical flanking path 
through the trim panel supports. 

Figure 5k also contains the sidewall empirical transmission 
loss spectrum (Figure 20) assumed for the baseline treatment. 
The baseline properties (A = 102 mm, p 2 = 1.76 kg/m 2 ) are 
similar to those for the lighter weight of the two analytical 
model sidewalls, yet the transmission losses are significantly 
lower. This difference between laboratory data and analytical 
predictions has been attributed [ 5 ] to the effects of flanking 
paths but there is, as yet, no confirmatory evidence. Little 
information is available on the vibration transmission through 
presently-used trim panel mounts, but, on the basis of [ 34 ], 
mounts with very low transmissibility will be required if the 
80 dB(A) goal is to be achieved. 

It has been shown in Section 6 that, in order to achieve the 
80 dB(A) criteria, noise reductions of about 25 dB relative to 
the baseline sidewall treatment will be required at the funda- 
mental blade passage frequency for each airplane. This increase 
in transmission loss is greater than the values shown in 
Figure 54 at corresponding frequencies. In fact, calculated 
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TABLE 10 


Predicted Weight Increase, relative to Baseline 


Airplane 

Airplane 

Weight 

(kg) 

Payload 

Weight 

(kg) 

Additional Sidewall Weight 

(kg) 

% Take-Off 
Weight 

% Payload 1 
Weight 

Wide-Body 

252,000 

47,500 

1,890 

0.75 

4.0 

Narrow-Body 

148,000 

24,500 

1,090 

0.74 

4.5 

Small-Diameter 

30,000 

2,350 

450 

1.51 

19.0 


Prom this analysis it is apparent that, at least for the ideal- 
ized double-wall system, the required reductions in cabin noise 
level can be achieved. Ho-wever, in order to achieve these 
exceptionally high transmission losses at low frequencies it is 
imperative that there be no structural flanking paths across the 
trim panel mounts or through the porous material in the cavity. 
The practicality of achieving this vibration-isolated double- 
wall system in an airplane environment has to be proven. 

7 . 5 Summary 

Add-on treatments considered in this study included double walls 
(with and without glass fiber batts), damping, skin panel stiff- 
ness, absorption in the cabin, and cabin partitions. Only the 
double wall with batts can provide the required noise reductions 
at low frequencies. Cabin partitions have no benefit because 
of the spatial extent of the power inflow, and add-on stiffness 
to the panel has little effect at low frequencies because the 
power inflow is determined by stringer and frame characteristics. 
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8.0 


ADVANCED NOISE REDUCTION METHODS 


In the present context, advanced noise reduction methods are 
those methods which involve modification to the fuselage 
primary structure. It is unlikely that such changes could be 
made as retrofit modifications to existing structures. Instead, 
they would be incorporated in new fuselage designs, and, as 
such, they would have to undergo the usual pre-certification 
testing before being placed in passenger service. The testing 
could be extensive if the new structures are significantly 
different from present-day aluminum skin-stringer-frame con- 
figurations. Because of the anticipated long lead time asso- 
ciated with the introduction of advanced noise reduction 
methods, they have been given less emphasis in the present 
study than have the add-on noise reduction methods , 

Since the fuselage frames determine the fuselage stiffness at 
low frequencies, modification to the frames was one of the 
important advanced noise reduction methods explored in the 
study. Also changes were made to the skin thickness. Varia- 
tions to the stringer stiffness were not considered in this 
phase of the study as the effects would lie between those due 
to the frames and the skin. However, some indication of the 
effect of stringer stiffness can be seen in Figure 86 in 
Section 9.5. Factors considered in the study of advanced 
noise control included frame spacing and stiffness, skin thick- 
ness, and the use of honeycomb panels. The calculations were 
performed using the baseline acoustic sidewall system as pre- 
sented in Section 4.5.4. 


8.1 


Variations in Frame ^Parameters 


8.1.1 Frame Stiffness 

frpmp bending stiffness parameter. 
The effects of doubling the frame oenamg 

n in Pieure 57. (This parameter gives the 
Erl/Eh, are shewn in Ergur 57 the skln parameters.) 

frame stiffness normalized with reap leve is 

The figure compares harmonic levels for the caD 

at tie location of maximum A-welghted level, using baseline 
at the iocati „ each study airplane show 

calculations as reference. “ ffnes3 results ln reductions in 
that the increased frame stiffness res 
cabin sound level, but the changes are sma 

The structural changes represented by this doubling ° f ^ame 
stiffness could be achieved either by the use of composite 

„ the frame construction or by redesign of the con- 
"al frame cross-section to increase the frame modulus^! 
inertia I r> while keeping the radius of gyra 1 r 
sectional area A r constant. 

Structural dimensions and fuselage parameters -o^d 
the doubling of frame stiffness are listed m 
13 for each of the study airplanes. 


n 


8.1.2 Frame Spaci ng 

A second approach considered was that of maintaining the frame 

stiffness equal to that of the baseline configuration, but 

doubling the number of frames by halving the frame 

structural dimensions and fuselage parameters --"f 

. rpnhlPF! 11 through 13 lor rne 

this modification are given in Tables ±1 

three study airplanes, and the changes in harmonic sound level 
arHlotled in Eigure 5 7 for the three lowest order harmonics. 
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Table 11 

FUSELAGE DETAILS FOR ADVANCED NOISE CONTROL STRUCTURES. 


WIDE BODY 


Descri pti on 
Element Number 
Frame Spacing 
Stringer Spacing 

Skin 


(m) 

Zy (m) 


Thickness (mm) 

Mass /Unit Area (kg/m 2 ) 
Bending Stiffness D ( N-m) 
Torsional Stiffness D^/N^m) 


Stringer 
Mass /Unit Area 


Bending Parameter 


p A /I (Wm 2 ) 

s y 

E sV D1 y 

Membrane Parameter E S A S / Ehl y 
Torsional Parameter G s J s / D1 y 
Centroid Z- (m) 


Frame 

Mass /Unit Area 
Bending Parameter 
| Membrane . Parameter 
Torsional Parameter 


pA /l (kg/m 2 ) 

E I /D1 
r r x 

E rV Ehl x 

Vr /D1 x 


Centroid Z 

! A .. x 10 3 
rib 


Basel i ne 

2x 

Frame 

E 

1/2 

Frame 

Pitch 

2x 

Skin 

Thi ckness 

4x 

Skin 

Thi ckness 

Honeycomb 

Skin 

4 

4 

4 

4 

4 

4 

0.5080 

0.5080 

0.2540 

0.5080 

0.5080 

0.5080 

0.1829 

0.1829 

0.1829 

0.1829 

0.1829 


1.778 

1.778 

1.778 

3.556 

7.112 

34.408* 

4.8006 

4.8006 

4. 800 6 

9.6012 

19-202 

5.7472 

38.133 

38.133 

38.133 

305.06 

2440.5 

13.012 

12.584 

12.584 

12.584 

100.67 

805.36 

17.474 

2.6680 

2.6680 

2.6680 

2.6680 

2.6680 


170.72 

170.72 

170.72 

21.340 

2.6675 

— ; 

0.5558 

0.5558 

0.5558 

0.2779 

0.1389 

— 

0.5084 

0.5084 

0.50 84 

0 . 06355 

0.007944 

— 

-0.01583 

-0.01583 

-0.01583 

-0.01672 

-0.01850 


2.8209 

2.8209 

5.6417 

2.8209 

2.8209 

1.3807 

4585.9 

9171.7 

9171.7 

573-23 

71.654 

4.1835 

0.5876 

1.1752 

1.1752 

0.2938 

0.1469 

0.3157 

0.6427 

1.2855 

1.2855 

0.08034 

0.010043 

0.00067 

-0.05468 

-0.005468 

-0.05468 

-0.05557 

-0.05735 

-0.03817 

0.2734 

0.2734 

0.3273 

0.2734 

0.2734 

— 


*Equivalent skin thickness based on radius of gyration. 


J 
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FUSELAGE DETAILS FOR ADVANCED 


Descri pti on 

Element Number 
Frame Spacing % (m) 

X 

Stringer Spacing l 

7 


Basel ine 

4 


(m) 


0.5080 

0.2286 


Skin 

Thickness (mm) 

Mass /Unit Area (kg/m 2 ) 
Bending Stiffness D (N-m) 
Torsional Stiffness (N-m) 


3.0 

8.100 

183.18 

60.45 


Stringer 

Mass/Unit Area pA /l (kg/m 2 ) 

s y 

Bending Parameter E I /D1 

s s: y 

Membrane Parameter E A /Ehl 

s s y 

Torsional Parameter G J /D1 
_ s s y 

Centroid Z (m) 
s ' 


1.0817 

23.687 

0.1335 

0.01264 

-O.OI855 


Frame 

Mass/Unit Area P A r /l x (kg/m 2 ) 

Bending Parameter E I /D1 

r r x 


Membrane Parameter 
Torsional Parameter 
Centroid Z Cm) 

A rib x 10 3 r (m 2 ) 


E r A r /Ehl x 

G r J r /D1 x 


1.6233 

294.55 

0.2004 

0.04455 

-0.03921 

0.1579 


*Equi valent skin thickness based on radius 


Table 12 

NOISE CONTROL STRUCTURES: NARROW BODY 


2 x 

Frame 

E 

1/2 
Frame 
Pi tch 

2 x 

Skin 

Thi ckness 

4x 

Skin 

Thi ckness 

Honeycomb 

Skin 

4 

4 

4 

4 

4 

0.5080 

0.2540 

0.5080 

0.5080 

0.5080 

0.2286 

0.2286 

0.2286 

0.2286 

- 

3.0 

3.0 

6.0 

12.0 

23.326* 

8.100 

8.100 

16.200 

32.400 

5.0195 

183.18 

183.18 

1465.4 

11,723 

5610.7 

60.45 

60.45 

483.6 

3868.7 

6.4379 

1.0817 

1.0817 

1.0817 

1.0817 


23.687 

23.687 

2.9608 

0.3701 

_ 

0.1335 

0.1335 

0.06677 

0.03338 

_ 

0.01264 

0.01264 

0.001581 

0.000198 


-0.01855 

-0.01855 

- 0.02005 

-0.02305 

- 

1.6233 

3.2467 

1.6233 

1.6233 

1.2389 

589.10 

589.10 

36.819 

4.6023 

6.6599 

0.4008 

0.4008 

0.1002 

0.05010 

O .3019 

0.08910 

0.8910 

0.005569 

0.000696 

0.001397 

-0.03921 

-0.03921 

-0.04071 

-0.04371 

-0.03698 

0.1579 

0.1929 

0.1579 

0.1579 

- 


of gyration. 
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FUSELAGE DETAILS FOR 


Descri ption 

Element Number 

Frame Spacing & x (m) 

Stringer Spacing ( m ) 

Skin 

Thickness (mm) 

Mass /Unit Area (kg/m 2 ) 
Bending Stiffness D (N-m) 
Torsional Stiffness D^. (N-m) 

Stringer 2 

Mass/Unit Area pA g /l y (kg/m ) 

Bending Parameter E s I s ' /D1 y 

Membrane Parameter E s A g ./Ehly 

Torsional Parameter G s J s / D1 y 

Centroid Z ( m ) 

Frame , 

Mass /Unit Area P A r /3 - x 

Bending Parameter E r I r / D1 x 

Membrane Parameter E r A r / Ehl x 

Torsional Parameter G^J^/Dl^^ 


Table 13 

ADVANCED NOISE CONTROL STRUCTURES: SMALL 


DIAMETER 


Centroid (m) 

A x 10 3 (m 2 ) 

rib 


Basel i ne 

2x 

Frame 

E 

1/2 

Frame 

Pitch 

2 

0.3302 

0.3302 

0.1651 

0.1270 

0.1270 

0.1270 

1.016 

1.016 

j 1.016 

2.7432 

2.7432 

2.7432 

7.1151 

7.H51 

7.H51 

2.3480 ! 

2.3480 

2.3480 

1.0973 

1.0973 

1.0973 

96.588 

96.588 

96.588 


0.4000 

0.4947 

-0.00744 


2x 

Skin 

Thi ckness 

4x 

Skin 

Thi ckness 

Honeycomb 

Skin 

0.3302 

0.3302 

0.3302 

0.1270 

0.1270 


2.032 

4.064 

22.886* 

5.4864 

10.9728 

3.6695 

56.921 

455.37 

3624.5 

18.784 

150.27 

5 .0663 


0.4000 

0.4947 


0.4000 ! 

0.4947 - 


1.0973 

12.074 

0.2000 

0.06184 


-0.00744 -0.00744 -0.00795 


*Equivalent skin thickness has 


1.0034 1.0034 2.0068 

2840.9 5681.7 5681.7 

0.3658 0.7315 0.7315 

0.4066 O.8132 0.8132 

- 0.03911 - 0.03911 0.03911 

I 0.07136 0.07136 0 .08252 

ed on radius of gyration. 


1.0034 

355-11 

0.1829 

0.05083 

-0.03962 

0.07136 


1.0973 

1.5092 

0.1000 

0.007730 

-0.00897 

1.0034 
44.388 
0 .09144 

0 .006353 

-0.04064 

0.07136 


J J 


It is seen that this method produces only small changes in 
harmonic level, and in some cases the sound level increases 
rather than decreases. 

8.2 Variation in Fuselage Skin Thickness 

Two situations were considered with regard to skin thickness. 
First, the skin thickness was doubled with respect to the 
baseline configuration. Then the thickness was doubled again 
so that it was four times the baseline value. The resulting 
fuselage details for the three study airplanes are tabulated 
in Tables 11 through 13. 

Changes in harmonic sound level associated with the doubled and 
quadrupled skin thicknesses are shown in Figure 58 for the 
three lowest-order harmonics and the three study airplanes . 

The spectra correspond to the location in the cabin with the 
highest A-weighted sound level. 

In general an increase in skin thickness results in a decrease 
in harmonic level for the cabin sound spectra. However, only 
in the case of the narrow-bodied airplane, which has a rela- 
tively thick skin in the baseline configuration, does the noise 
reduction exceed 3 dB for either the first order (fundamental) 
harmonic or the A-weighted sound level. As the skin thickness 

increases within the range considered, it is found that the 
structural modal density at low frequencies increases but, at 
the same time, the band-averaged joint acceptance decreases. 

The increase in modal density occurs because, at these fre- 
quencies, the fuselage stiffness is dominated by the frame 
characteristics and the increase in skin thickness is observed 
mainly as an increase in structural mass. 
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These changes in modal density and joint acceptance essenti- 
ally cancel each other for the wide-body and small diameter 
cases and the resulting noise reduction is small. For the 
narrow-body airplane, however, the baseline skin mass is rela- 
tively large so that the cabin noise levels are more sensitive 
to further increases in skin thickness. 

8 . 3 Honeycomb Skin Panels 

During recent years there have been several investigations 
concerned with the use of honeycomb panels for fuselage primary 
structures, although no such structures have yet been intro- 
duced for passenger-carrying aircraft. One advantage of honey- 
comb construction is that fuselage stiffness characteristics 
can be maintained whilst reducing the total weight of the 
fuselage. However, the noise transmission characteristics of 
the honeycomb panels, when exposed to typical fuselage exterior 
pressure fields, have not been explored. Consequently honey- 
comb skin panels were included in this analysis of advanced 
noise control methods. 

It was not possible, within the scope of this program to de- 
velop an analytical model for noise transmission, to derive 
general design characteristics for honeycomb fuselage struc- 
tures. Thus, reference was made to. available structural in- 
vestigations. In particular honeycomb structures were designed 
on the basis of fuselage configurations presented in [29]. The 
honeycomb structures have ring frames with the same pitch as 
for the corresponding baseline study aircraft, but no stringers 
are used. A typical panel configuration is shown in Figure 59, 
and dimensions are listed in Table 14. Structural character- 
istics are presented in Tables 11 through 13. 
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TABLE 14: 


DIMENSIONS FOR HONEYCOMB STRUCTURES 


uminHiiiiiiimiiilii 


WIDE-BODY 

FUSELAGE 


NARROW-BODY SMALL-DIAMETER 
FUSELAGE FUSELAGE 


CORE THICKNESS 

h 

c 

(mm) 

FACE PLATE THICKNESS 

h 

P 

(mm) 

FRAME THICKNESS 

h f 

(mm) 

FRAME DEPTH 

a l 

(mm) 

FLANGE LENGTH 

a 2 

(mm) 

FLANGE LENGTH 

a 3 

(mm) 

WEB-HONEYCOMB SEPARATION 

a 4 

(mm) 

FRAME SPACING 


(m) 


12.70 

0.76 

1.40 

88.9 

25.40 

25.4 

6.35 

0.508 


12.70 

0.51 

1.02 

63-5 

20.32 

20.32 

6.35 

0.330 








FIGURE 59. REPRESENTATION OF FRAME-STIFFENED HONEYCOMB 
FUSELAGE SHELL 


9 


. , oH with the cabin maximum A-weighted 
Predicted spectra associated with pne 

sound level, for each study airplane with honeyed structure 
are presented in Figure 58. It is seen that, for the se e 
structures, there is a general, although small \ pect lon 

sound level relative to the haseline configuration Inspec ^ 
Of the analysis indicates that the noise increas ^ 

decrease in smeared-out mass of the structure. 

of light-weight fuselage structures with bending and membrane 

® . , . j_ 0 -t-bose of conventional struc- 

stiffness parameters equivalent • , 

tures and with similar fuselage loss factors will lead in 
general to increases in interior sound -levels. Optim za ion 
honeycomb structures with respect to noise transmission will 
require extensive use of the analytical model and development 
dynamic models more suitable for honeycomb structures. 


8.4 Summary 

The results of the study of advanced noise control methods 
presented in this section are summarized in Table 15 and 
Figure 60 in terms of changes in A-weighted sound level and 
the surface mass density of the structural element which trans- 
mits the major fraction of the acoustic power inflow. The 
total mass per unit area is quoted in the table, and the change 
in surface density, relative to the baseline structure, s 
shown in the figure. In all cases it is assumed that the 
sidewall treatment is the same as that for the corresponding 
baseline study airplane. 

In general the decrease in the A-weighted interior noise level 
is proportional to the increase in the mass surface density 
of the structure. A 3 dB decrease in noise level is produced 
approximately by a 3 dB (doubling) increase in mass density o 
the fuselage structure independent of whether the mass addit on 


n 
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TABLE 15 


SUMMARY OF NOISE REDUCTIONS ACHIEVED WITH 
ADVANCED CONCEPTS 


Structure 

Wide 

Body 

Narrow Body 

Small Diam. I 

Mass/ 
Unit Area 

Max. Cabin 
dBA 

Mass/ 
Unit Area 

Max. Cabin 
dBA 


Max;. Cabin 
dBA 

Baseline 

9.46 

105-5 

10.29 

104.9 

4.84 

108.1 

2 x Frame Stiffness 

9.46 


10.29 

104.7 

4.84 

104.2 

Half Frame Pitch 

11.05 


13.11 

106.0 

5.85 

106.4 

2 x Skin Thickness 

16.2 

104.4 

15.09 

102.2 

7.59 

106.6 

4 x Skin Thickness 

29.7 

101.9 

24.69 

97.6 

13.07 

104.2 

Honeycomb Skin 

6.26 

108.1 

7.13 

108.8 

4.76 

107.8 


Mass/Unit Area (kg/m 2 ) Includes Skin, Stringers and Frames For 
the Structural Elements For Which the Input Power is a Maximum. 
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occurs as a result of Increasing the number of stiffening ring 
frames, or as a result of increasing the skin thickness. There 
are, of course, some exceptions to this general trend. For 
example, a 4 dB reduction in noise level is predicted for the 
small diameter airplane when the frame stiffness is doubled, 
yet this method does not increase the structural mass. 

From this analysis of advanced noise control concepts it is 
apparent that it will be difficult to achieve the required 
noise reductions by means of these methods alone. Whilst a 
more extensive parametric analysis is required than was possible 
under this study, it is anticipated that add-on noise reduction 
methods will have to be used with the advanced concepts in 
order to meet the 80 dB(A) goal for interior noise levels. An 
optimization procedure will be required to perform the analysis. 
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9 0 REQUIREMENTS FOR EXPERIMENTAL VERIFICATION OF 
ANALYTICAL MODEL 

9.1 Introduction 

The analytical model developed in Section 3 has been used t 
evaluate the various add-on and advanced noise control "etho 
described in Sections 7 and 8. However, it is adV1 “ blb 
the analytical model undergo some form of experimental va 
tion just as was done for a similar analytical model of the 
acoustic environment in the Space Shuttle orbiter payloa 
bay [10, It ] . The key questions to be answered in the va l 
tion of the present analytical model are: 

(i) Does the analytical model adequately predict 
the response of the fuselage structure to 
convected pressure field with a rapidly decay- 
ing amplitude? 

(ii) Does the analytical model adequately predict 

the sound levels in the interior of a fuselage 
with a sidewall treatment typical of present- 
day turbofan aircraft? 

(iii) Can the add-on sidewall treatments provide the 
high noise reductions predicted by the analyti- 
cal model at low frequencies? 
and 

(iv) Does the analytical model adequately predict 
the noise reductions provided by different 
fuselage structure designs? 

Practical questions concerned with the design of validation 
tests revolve around several key steps: construction o a 
model fuselage which adequately simulates the full-scale 
propfan airplane; simulation of the propfan excitation; 
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and definition of the experimental measurements and procedures 
to be followed to answer the above questions. 

In this section, criteria for adequate simulation of typical 
propfan airplane configurations are discussed in detail. In 
Section 9.2, sensitivity studies are conducted to demonstrate the 
Importance of excitation, structural, and cabin variables on the 
measured noise reduction. Alternative excitation simulation 
methods are reviewed in Section 9.3, bearing in mind the criteria 
developed in Section 9.2 for achieving the above test program ob- 
jectives. Section 9 A presents detailed recommendations, resulting 
from the preceding discussion, for the proposed test program. 
Excitation and fuselage models, and measurement and data analysis 
procedures are outlined. The relationship between experimenter 
and analyst is discussed in detail. Validation using flight tests 
with an existing propeller-driven airplane is discussed in 
Section 9.5. 

9.2 Sensitivity Studies 

The planning of a model test study requires the selection of an 
appropriately constructed and scaled test model. Fundamental to 
correct model selection is an understanding of the sensitivity 
of^the calculated noise reduction (and the parameters such as 
<j r > p , ff, n p , n r , t s ) to variations of the characteristics of the 
excitation and the structural and sidewall configurations from the 
prototype propfan airplane designs. In this section, scaling laws 
are presented which show how the characteristics of the excitation, 
structure, sidewall acoustic treatment , and interior acoustics 
depend on the size and material selection in the model fuselage 
relative to the propfan prototype. Approximate limits within 
which critical parameters must be maintained, are presented for 
guidance in the practical: selection of a scale model experiment. 
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9.2.1 Model Scaling 


The Mikulas resonance-frequency equation (Section 3.8.1), 
which smears out 'the mass and stiffness characteristics of 
frames and stringers over the skin surface, has been used to 
predict the wavenumber-frequency distribution for the struc- 
tures of the various study airplanes. The equations were 
modified to include the effects of cabin pressurization. 

The contributions of the skin and stiffener masses and stiff- 
ness are clearly represented in the equations (Figure 5 ) in 
a series of non-dimensional groups which can be readily varied 
independently to quantify the influence of each component in 
the prediction of fuselage wavenumber-frequency distributions. 

Correct scaling implies that the test model fuselage has the 
same non-dimensional wavenumber-frequency distribution as the 
prototype fuselage. The resonance frequencies v mn of the 
model fuselage, when non-dimensionalized by the theoretical 
cylindrical fuselage ring frequency, must remain constant for 
all vibration modes. The Mikulas equation (Figure 5), both 
sides of which are non-dimensional, can be written m abbre- 
viated form, as 


2 

V 

mn 

4 

n 


h l 
e 


' 3 2 

h R 


function 


E I 
r r. 


E A 
r__r , 


M x E* x h 



(84) 


(all non-dimensional groups are listed in Figure 5.) . Here 
h is the skin thickness, h g is the thickness of the equiva- 
lent orthotropic panel (i.e. skin plus smeared-out stiffeners) 
and v = w /to where u> = c L /R and c is the longitudinal 
wavespeed in the skin. It is assumed in the following sections 
that the skin and stiffeners are of the same material. 
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For v mn to be const ant for model and full-scale fuselages 
requires a dimensionally-scaled model where h/l, h/R and h /R 
are maintained constant. Then the non-dimensional groups on 
the right hand side o'f Eq.(84) must be maintained constant, 
which requires that all dimensions are scaled to maintain the 
same ratio with the fuselage radius. Dimensionally-scaled 
models, in which all dimensions are changed in the same ratio, 
will maintain the same non-dimensional wavenumber-frequency ■ 
distribution according as 

w mns /w mnf = R f //R s 

using the same material or, more generally, a material with the 
same compressional wavespeed c L for both model and full-scale 
fuselages, and where is the (m,n)th resonance frequency of 

the model cylindrical fuselage of radius R g • subscript f refers 
to the full-scale fuselage. Thus, for example, use of a half- 
scale test model requires that measurements be carried out at 
twice the excitation frequency of the full-scale fuselage. 

By analogy with uniform axisymetric circular cylinders, the 
non-dimensional modal density N(v) is given by N(v) = 
n (^ b ) Wr/(Lx/ic) , when k is the effective radius of gyration of 
the stiffened cylindrical fuselage. L x /k is constant for a 
dimensionally-scaled model fuselage so that nU)u> r must also 
be constant for model and full-scale fuselages when evaluated 
at the scaled excitation frequency if N(v) is to be constant. 
Figure 6l shows that n(o)^)/R is constant for model and 
full-scale structures when correctly scaled dimensionally, 
minor variations occur due to band-width computation proce- 
dures at low frequencies. In performing this analysis, it 
is assumed that the model and full-scale fuselage structures 
are constructed of the same material, more particularly, mat- 
erials with the same value of Detailed examination of the 

. -199- 


200 


Full scale Ap=55.2kPa 

Half scale Ap=55.2kPa 

Half scale Ap= 0*0kPa 


L /R = 1 .69 L /R = 1 .74 
x y 

h/R = 5.9x10* 


a L = 3.56 c = 0.01 

V y -’- ls v 0 - 01 

U x = 745m/s U y = 498m/s 


10 ' 


-Q 10' 


FIGURE 61 


Harmonic Order 

EFFECT OF STRUCTURAL SCALING ON <j*>, n(w b )/R 5 and <0> 
Wide Body Airplane: Full scale f f = 290Hz ; 


• scale 


l 




r 

r 

resonant modes used in computations show that an essentially 
!~ identical wavenumber-frequency distribution is produced for 

j 

1 a correctly scaled model fuselage when evaluated at frequency 

r 




r 


i 


9.2.2 Noise Reduction Scaling 

The high frequency formulation (Section 3*5.4) is used ex- 
clusively for all study airplanes. For baseline aircraft 
where W t (u>) can be neglected in comparison with W abs '(u), 
the fuselage noise reduction can be calculated, from Eqs . 
(26) and (54), as 


NR(co b ) 


10 log <p Xo ( ' a) b )> 

io r, r~ 
<Pi 


(85) 


where 

< PjCo (M b )> 

< Pl ( “b )> 


8ir(p . c . ) 

l l 


A 


2 -, 


co a) y aS 

D r 


2 

[n ( u). ) u) ]<j (di, ) >t, 
r b r r b P 


<a t I 

_ rad s } 


( 86 ) 


The excitation frequency is denoted in Eq.(86) by , where 
it is implicitly assumed, from Section 9.2.1, that the ratio 
(w b A> ) is the same for model and full-scale conditions. Thus, 
when it is stated that some function, say $(a> b ), is the same 
for model and full-scale conditions, what is really implied 
is that $(a> b /u) r ) remains constant. An equivalent approach 
is to specify that 0 ( 0 ^) is constant for constant harmonic 
order b. 
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As will be shown, parameters contained within the { } brackets 
are all non-dimensional, depending only on non-dimensional 
groups which are, or can be, maintained constant in scale 
model testing. The terms preceding the { } brackets deter- 

mine the dependence of the noise reduction on dimensional material , 
and frequency changes: these, if any, must be clearly defined if 

model acoustic tests are to be properly interpreted and used 
for full-scale noise reduction predictions. The variables on 
which the important parameters determining the structural 
response of, and acoustic radiation from, the fuselage depend, 
are briefly reviewed. 

For simplicity, the following discussion will present numerical 
examples based only on the wide-body fuselage structure, but 
similar conclusions can be drawn for other aircraft. Also, 
the discussion is limited to the three lowest order harmonics 
since these dominate the interior noise levels. 


9> 2. 2.1 Scaling of the external a c oustic excitation 

The excitation is described by its cross spectral density 
function which depends on the pressure amplitude decay rate 
parameters (a L , a L ) , the excitation coherence parameters 
(c c ) and t£e non-dimensional excitation trace wavenumber 
(k X R, y k R). Each must be maintained constant between model 
and full-scale fuselage to maintain correct excitation scaling. 
For example, for a half-scale model fuselage, the pressure 
amplitude decay (in dB/m) for the fuselage model will be twice 
that for the full-scale fuselage since the model is half as 


a xf L xf /L xs 


long as the full-scale fuselage: i.e. a xa 
coherence decay parameters (c x ,c ) can be written as 
. = c > (“I) where c ' is a constant. The Strouhal number 

X IT ^ — • i TT J — . rt (- Ol 

w5/U 


The 


x 


will X be constant if the trace velocity U x is constant. 
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since, as determined earlier, length scales (such as £) vary 
inversely with frequency. (The selection of c is somewhat 

X 

arbitrary on account of a lack of experimental data.) Further, 

as will be shown, the sensitivity of the joint acceptance to 

variations in (c , c ) is weak when(c ,c ) are small. For 

x y x y 

constant trace velocities (U ,U ) in the axial and circumfer- 

x y 

ential directions respectively, (k R,k R) will be constant for 

x y 

a scale model since co.^R will be maintained constant. Further- 
more it is imperative from the point-of-view of excitation 

source directivity, that the acoustic wavenumber to R/c be 

b e 

maintained constant. These various requirements are consistent 
with each other and with the structural scaling requirements. 


9. 2. 2. 2 Band-averaged joint acceptance <j£(.u) b )> p : 

The band-averaged joint acceptance for a progressive wave 
excitation is calculated for a group of resonant structural 
modes by summing the individual joint acceptance functions 
at the analysis-band center frequency and then dividing by 
the number of modes. The variables, on which the modal joint 
acceptances depend, are given by 


jmm^b^ = f unctlon 


where 


k x R 

k m R 
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- b -x 
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k k R » 
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k R 
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(m,n) , (a L ,a L ) , (c ,c ) 
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GO. L 
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y 


(m,n) is the struc- 


tural mode order, (k m ,k n ) the structural modal wave numbers. 


If (aL ,a L ), 
x x» y y * 


and (k ,k ) the excitation wavenumbers 

. «y * ^ j j 

( c x> c y) and (U x 3 U y ) remain constant, as required to correctly 

simulate the excitation characteristics on a test model, and 

the model is also dimensionally scaled, (m,n) , (k R/k R) and 

x m 
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(k R/k R) Will also be maintained constant at the scaled 
excitation fluencies and correct simulation of the joint 
acceptances of the full-scale propfan will be achieve 
Figure 61 presents calculations of the band-averaged joint 
a Iptance for full-scale and half-scale fuselage elements 
for the wide body study airplane for in-flight pressurization 
conditions. The model and excitation have been scale i 
accordance with the previous discussion: it is clear 
<■,£( )> p is constant for constant harmonic order b. 

9.2. 2. 3 RanH-averaqe d radiation efficiency 

The band-averaged radiation efficiency 

<a > depends on the individual modal radiation 

which themselves are functions of mode order relative struc- 

and structure dimensions, i.o, 
tural and acoustic wavespeeds, and struct 

<0 > = function £(m,n), ^Gi^^mn^' ^ci 5 ci 

rad 

a w r - „ R/c . A correctly scaled model will preserve 

th^same^non-dimensional wavenumber-frequency distribution 

for model and full-scale structures while k cl L and k cj R will 
also be constant at the new model-scale test frequencies 
Thus <„ „> will remain constant for correctly-scaled model 

fuselagfstructures. This is also shown in Figure 61. 

9. 2.2. 4 Re nd-Averaq e d structural Loss Fact or.ii: 

The band-averaged structural loss factor 

tural dissipation associated with factors such as joint 
tion, gas pumping, and applied damping treatments, as well as 

on radiation damping forces. Thus 

t = n struc + n rad 
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where, as discussed in Section 5.6 , ^struc ^ e P en ^ s essenti- 
ally on the structural mode order, provided the joint details 
are similar and the frequency change is not too great, while 
the acoustic radiation loss factor ti^ d depends on both struc- 
tural surface density and radiation ratio: i.e.. 


'rad 


p . c . <0 , > 

1 1 rad 


7 v 


where y is the surface density of the fuselage structure. Since 

dimensionally-scaled models will maintain (m,n) , u>.b . v and 

s * b mn 

<a rad > cons ta-nt at the new modal, test frequency then n ra( j, ^ s ^ ruc 

and consequently n will all be maintained constant. 

v » v 

9. 2. 2. 5 Sidewall Transmission Coefficient t : 

, _s_ 

The sidewall transmission coefficient r g of a double-wall 
structure, is essentially a function of the acoustic and 
double-wall resonance wavenumbers, non-dimensionalized with 
respect to cavity depth i , as well as the attenuation and 
wavenumber parameters for the medium in the double-wall 
cavity. Thus, 

t g = function [k cl £ = , k dw £, aZ, 3*] 

If each of these non-dimensional parameters is maintained 

constant between model and full-scale fuselages, r will be 

s 

constant. 

The non-dimensional acoustic wavenumber k C j_£ will be held con- 
stant for constant since u£ is constant for a dimensionally- 
scaled model fuselage. 
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The double-wall, non-dimensional resonance wavenumber k dw t is 


given, approximately, by 


k dw^ 


1 / 2 

/ pt(y 1 +P 2 )\ 

' II 11 . ' 


y 1 V 2 


(87) 


J \ Pi + fil Y 

\ P l h l P 2 h 2 / 


/ : 


where y ls Ua are the material densities for the fuselage and 
trim panels respectively: h l# h 2 are the corresponding panel 
thicknesses, p is the density of the cavity medium. Then if 
the materials are the same for model and fullscale panels, 
the model and fullscale values of k dw & will be the same when 


attenuation'parameter at where a is attenuation per unit dis- 
tance (dB/m) in the porous medium (see Section 3.9.2), may be 
written non-dimensionally as 

R ’*- •'-°- 595 l (88 ) 


and h are all scaled in proportion to radius R. 


The 


at — k t 
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r r 1 1 i 

„56H L k V("l^)J 


while the porous material wavenumber parameter W, where 
g . 2 ,/( acoustic wavelength In the material), can be written 

non-dimensionally as 

-0.700 


g t = k t 


1 R t .. 

1 + 0.35 [ k«/(5T")J 


( 89 ) 


Maintaining the acoustic wavenumber M and the non-dimensional 
flow resistance (R^/pc) constant will ensure and $«. 
are properly scaled and maintained constant. 

Consequently the double-wall sidewall transmission coefficient 
and the add-on sidewall transmission loss will be the same for 
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model and full-scale fuselage configurations. The requirement 
that the flow resistance remain constant for the model 
sidewall Implies a change In flow resistivity R x : this can be 
achieved by changing either the material density p or the 
material fiber diameter d since R 1 «p m 1 - 53 /d 2 £26]. Changes 
In p m are most easily arranged by compressing standard porous 
material samples, although there will be a small increase in 
the total weight of the material. Such an increase will be 
small compared with the double— wall mass, and will not cause 
a significant change in the noise reduction provided by the 
wall panels, provided that the mechanical stiffness of the 
layer is not so high that it invalidates the acoustic modeling. 
Detailed calculations, as seen in Figure 62 for full-( curve (a)) 
and half-scale (curve (b)) models where T g (k£) remains constant 
for and R-^£ constant, verify the scaling arguments. If 

R 1 A ls not scaled correctly, r g (k£) will not be correctly 
scaled as shown by curve (c) in Figure 62. 

9. 2. 2. 6 Interior Acoustics : 

The acoustic modal overlap of the cabin volume, expressed by 

, n ( “b ) = "n n n ( “b J “b 
a(ui b ) S k b 2 


must be kept constant to ensure correct scaling of the side- 
wall radiation characteristics. In a dimensionally-scaled 
model, k b S <* (k^R)(k^L) will remain constant. Then, to keep 

the modal overlap unchanged, the absorption coefficient of 
the various surfaces must be the same for the model and pro- 
totype fuselages, in the above equation, k b = 
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Add-on Sidewall Transmission Loss, dB 



One-Third Octave Band Center Frequencies in Hz 


GURE «. PREDICTED ADDITIONAl ^SI^WALL ^NSM^SS.ON 




The absorption coefficient depends essentially on sidewall 
characteristics since the bulk of the cabin absorption is 
provided by the bounding surfaces: for example, 

a(to^) = function 

all terms of which are non-dimensional and remain constant 
for a dimensionally-scaled model (y 2 is the interior trim- 
panel surface density). Thus a will be constant if remains 
constant, since the same non-dimensional parameters are in- 
volved in both cases. 

In practice, the cabin length L may be shorter than dimensional 
scaling would allow: additional absorption can then be added 
to non-radiating surfaces so that ¥ (w b ) remains constant. 


V> V-“ bP ? 

p i°i p l°i I 


9. 2. 2. 7 Summary : 

Referring to Eq. (86) for the fuselage noise reduction, it has 
been shown that, for a dimensionally-scaled model and ex- 
citation, the terms inside the { } brackets will be equal to 
those for the full-scale fuselage. The change in noise reduction 
from full-scale to model-scale then reduces to 

2 o 
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(90) 


for a dimensionally-sea:iqd' model, where, according to the 


- 209 - 


analytical model, subscripts s and f refer to model-scale 
and fullscale conditions, respectively. 

Thus, the measured noise reduction on a dimensionally-scaled 
test model will he essentially equal to that measured on a 
full-scale fuselage at the same non-dimensional frequency 
provided the excitation has been properly scaled. 


9.2.3 


Ff fp C ts of variations in exc ltatjon^^ 
characteristics. 


In a practical model test situation it may not be feasible 
and/or economic to construct a fully-scaled test mode . or 
example, exact simulation of stiffener details would be a 
relatively expensive task. Thus it is necessary to n 
magnitude of the errors produced in noise reduction calcu - 
tions or, more directly, in the parameters shown in Eq. ) 
on which the noise reduction depends, when other than comp 
simulation exists. The most important of these parameters 
are the band-averaged Joint acceptance, the structural modal 
density and radiation resistance, and the sidewall trans- 
mission coefficient. 

Consider first the effects of pressurization loads. Figures 
61 and 63 show that for the wide body fuselage and for re 
auencies below the subpanel resonance frequency, <J r ( > P » 

nU )/R and ”<* ^pend »lgnificantly on whether 

full-scale pressurization loads are simulated m the model 
scale test. This is true because the frames and stringers 

,. d-Mffness relative to the induced in-plane 
have high bending stiffness reiaui fmmd 

(pressurization) membrane stiffness. Similar resu 

for the other study aircraft, as seen in Figures 64 and 65 . 
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Thus there is little need to pressurize a test model fuselage 
if one is concerned only with low frequency noise reduction, 
for example, at the lower three harmonics which are of present 

concern. 

Of interest now is the accuracy with which the excitation 
characteristics, (a x ,a ), (U X ,U ) and (c x ,c y ) must he modeled. 
Consider the dependence of the band-averaged joint acceptance 
on the amplitude decay rate parameter (a x ,a y ) . Figure 66 
shows, for the wide-body airplane (baseline structure, ele- 
ment 4), the dependence of <j^> p on variations in (a x »a y ) for 
a range of harmonic order, (U x ,U y ) and (c x ,c y ) being kept 
constant. The dependence is marked, as predicted by Eq . (A. 9) 
of Appendix A. The values of <j*> p relative to the baseline 
estimate are shown in Figure 67 for different values of (a x ,a y ). 
It is clear that close simulation of (a x ,a y ) is most important; 
for an accuracy of ±1 dB in <j r >p and s hence, in structural 
response and noise reduction simulation, ( a x s a y) mus ^ 
modeled with an accuracy of ±10#. However, as shown in 
Figure 66 for situations involving strong decay in the pressure 
amplitude, simulation of (c x ,c y ) need not be very precise. 

Figure 68 shows the dependence of <J r > p on variations in trace 
velocities (U ,U ) for the wide body aircraft: it is clear 
that close simulation of U x , and U y is imperative, and limits 
of ±10# are suggested as reasonable. A similar strong 
dependence on the variation in U x ,U y is found for the other 
study aircraft, although the errors do not follow any general 
pattern since basic structural differences exist between 
different study aircraft. 
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FIGURE 68. EFFECT OF VARIATIONS IN TRACE WAVE SPEED ON 
BAND-AVERAGED JOINT ACCEPTANCE: WIDE BODY 
AIRPLANE - ELEMENT 4: LOW FREQUENCY STRUCTURAL 
MODEL 



The effects of incorrect!* ^; i ^r;it^t B i»“ 

parameters are now discussed. Consider 

— - cnr" n ::;’s=rc: LLi*. r 

- -- — re " ai ; r - r. :l:x 1. 

variations in <c rad >» n( * b ) a " ? basellne fra me 

plane which result from variations > ocour sln ce the 

bending stiffness. Small variations rad 

m odes are still supersonic for only 

of the modal density ben41ng stiffness is 

doubled "Tt seels that the frame stiffness and spacing should 

be fairly closely scaled, within ± 20 ? of ^ "lona ^ 

On the other hand, as seen in Fig 
scaled values. the frame membrane 

errors introduced by improperly modeling of 

— t o ( f ^ - 

b f :i:“:: 0 :; :::-sr::r;he combined and 

torsional stiffness parameter both appear to have on y mm 
influences on either parameter. 

Variations in stringer stiffnesses for the widebody airplane 
are shown in Figures 71 and 72 . The stringer bending stiff- 
ness, membrane stiffness and radius of gyration are most im- . 
portant in the intermediate or mid- frequency structural regime 
(where only the stringers and skin properties are considered). 
However, from Figure 72 , it is clear that stringers must be 
included in the test model for good simulation of n(» b ), even 
at low frequencies, although it appears that precise dimensional 

scaling is not necessary. 

Variations in the wide-body skin thickness produce substantial 
changes in response parameters in the mid- and high frequency 
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regimes where only stringer and skin properties are included 
in the fuselage model. Increases in < ^ rad > » nU b >tand <J r > p 
occur due to increases in skin thickness, as seen in Figure 73* 
Similar ■ trends occur for the other study aircraft. These 
increases are counterbalanced to some extent by increases in 
the ’smeared out’ fuselage surface density, but increases in 
skin thickness can lead to decreases in the fuselage noise 
reduction (see Figures 58 and 60 ). Since the fuselage noise 
reduction, at least for the excitation spectrum given by the 
Hamilton Standard prediction procedure- [6 ], depends on the 
lower order harmonics and since at these frequencies the 
noise reduction of the study airplanes is controlled to a 
large extent by stiffener properties, precise simulation of 
the full-scale skin thickness is not necessary: variations 
of ±25$ from dimensionally-scaled values would seem to be 
acceptable. 

The discussion in Section 5.1 shows that panel curvature con- 
trols the wavenumber-frequency distribution of the fuselage 
structure at low frequencies. Figure 7^ shows the effect 
on the various parameters of interest of using a flat panel 
to simulate a curved panel. Thus, it is clear that flat test 
panels are not appropriate as model structures at low fre- 
quencies . 

The area of the test fuselage should be large enough that 
strong decay in the pressure amplitude occurs between the 
peak location and the panel boundaries. As shown in Figures 
30 and 33 , and as discussed in Section 5.2, significant 
underestimates of the panel modal density parameter, n(u> m )/A 
and joint acceptance parameter, result from using 

too small a panel area. The criterion suggested by Eq . (A. 9) 
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of Appendix A, that L x >a x 1 and ^ , must be used in 

selecting appropriate panel dimensions both circumferentially 
and axially over the fuselage surface. Thus, while the analy- 
tical model does not of itself require a cylindrical test 
section, it does require that the panel elements have the 
correct curvature and have dimensions which satisfy the decay 
rate criterion. Correct simulation of the fuselage structural 
response and associated noise reduction will require in general 
a model fuselage which extends circumferentially from floor to 
floor. 

9.2.4 Summary of Practical Model Simulation Requirements 

The parametric studies in Sections 9.2.2 and 9.2.3 have pro- 
duced requirements for the design of validation tests for the 
analytical model. These requirements can now be summarized 
for easy reference. It is also important at the same time to 
reiterate the assumptions which formed the basis for the 
requirements . 

Firstly, it was assumed that the tests will attempt to model 
as closely as possible a fullscale baseline fuselage structure 
and sidewall treatment, so that it will not be necessary to 
perform extensive analytical work to complement the measure- 
ments. Secondly, since the noise control problems are most 
severe for the low order harmonics, emphasis is given to the 
simulation of the fullscale conditions at the model scale 
frequencies associated with the low order harmonics. 

Requirements for the model are as follows :- 

Excitation: Trace velocities and amplitude decay rates 

should be simulated closely (within an accuracy 
of ±10%). Coherence decay can be simulated less 
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accurately. (Note that, without superimposing 
a mean airflow, the trace velocity will be 
supersonic, l.e. it will represent only that 
portion of the pressure field which Is In 
the plane of rotation or downstream of the 
propeller. This is the region of high acoustic 
transmission) . 


Structure: The structure should be curved and have frames 

and stringers. The structure should extend from 
floor to floor, and be at least two fuselage 
diameters long. The stiffness, mass and spacing 
of the frames and stringers should be simulated 
closely, (within an accuracy of ±20$). Simula- 
tion of skin mass is less important (an accuracy 
of ±25$). In-plane loads due to fuselage 
pressurization need not be simulated. 


Sidewall Simulation of the add-on sidewall treatment is 
Treatment, most important. Structural flanking must be avoided. 


Cabin 
Spaces : 


The model overlap of the acoustic space should 
be simulated closely. However, details of the 
cabin space, such as location of absorption, and 
shape of furnishings and passengers, are of low 
importance: the length need not be scaled provided 
that the absorption is scaled. 


9.3 Candidate Noise Sources 

The analysis performed in Section 9«2 indicated that the noise 
transmission through the fuselage structure was strongly depend- 
ent on the trace velocity and spatial decay of the excitation 
field. Thus the selected model sound source should be capable 
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of representing these properties fairly accurately. It should 
be noted that the analytical model does not distinguish between 
a trace velocity associated with the propagation of acoustic 
waves and a trace velocity associated with forward motion of 
the fuselage. Thus it is not necessary to introduce mean flow 
over the model structure, although if the mean flow is zero it 
will not be possible to model the subsonic trace velocities 
which occur over the forward regions of the baseline airplanes. 
These regions are, however, of lesser importance in terms of 
noise transmission. 

Several alternative noise sources have potential for use in 
simulation of the excitation. Propeller noise sources include 
the various model propfans, and full or model-scale propellers, 
either operated statically or in the presence of air flow, 
such as provided by a wind-tunnel or by in-flight conditions. 
Electro-acoustic sources, such as a single or an array of 
suitable loudspeakers, have found application as simulation 
devices in some verification programs [ 10 3 . The arguments for 
candidate sources are reviewed herein. 

9.3.1 Model Propfan 

A model propfan has been used in NASA Lewis/Hamilton 
Standard test programs. The model has the correct geometry 
and number of blades (at least as far ar the current state of 
technology is concerned) and would be operated at scaled rpm 
to maintain correct tip Mach number. It could be operated 
alongside a small model fuselage in a wind tunnel, although to 
reach a flow Mach number of 0.8 it is probable that a test 
section without acoustic treatment (and with attendant data 
interpretation problems) would have to be used. Tunnels with 
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acoustically treated test sections generally operate at Mach 

numbers no greater than 0.35. It should be noted that 
the simulation of propfan performance by overspeeding the model 
propfan to achieve the helical tip Mach number associated with 
cruise conditions does not necessarily simulate propfan 
acoustics correctly. It has been claimed that such a situation 
does not correctly simulate pressure amplitude, and that it 
will not simulate directivity and trace velocities either. 

Firstly, the main disadvantage of the model propfan is that, 
to maintain the correct ratio of fuselage diameter to propeller 
diameter, the model fuselage would be very small thereby making 
it very difficult to scale the fuselage structure and acoustic 
treatment. For example, considering the narrow-body airplane, 
(D/<j> = 1.13 and r/D = 1.8) use of the SR-3 model propfan £43] 
with D = 0.622m requires that the model fuselage be 0.55m or 
21.5 inches in diameter. Then, maintaining the ratio of add-on 
sidewall treatment depth l to fuselage diameter, a/* = 0.0*12, 
the model fuselage would require a sidewall depth of 6.3mm and 
an associated flow resistance for the porous infill in excess 
of 10^ mks rayls/m. Such scaling of the sidewall is impractical. 
Furthermore, effects of the measurement processes, such as the 
presence of a microphone array in the cabin volume, or (light- 
weight) accelerometers on the model skin, might lead to 
erroneous results. 

Secondly, even though the propeller and fuselage would be small, 
it would still be difficult to install both of them in a wind 
tunnel test section (open or closed) with a representative 
distance between propeller tip and fuselage, and with adequate 
clearance between tunnel wall (or shear layer) on one hand, and 



propeller tip and fuselage wall on the other. For example, 
using the above example with D/<|> = 1.13, r/D = 1.8 and 
D = 0.622m, and assuming a shear-layer thickness of 150mm, it 
is clear that only half the model fuselage can be immersed 
within the open jet flow of the anechoic flow facility at 
NSRDC , Carderock [44], as shown in Figure 75 . (The maximum 
flow Mach Number of the tunnel is M 0 ~ 0.2). Larger fuselage 
models compound the problem. A similar situation exists with 
the United Technologies facility [45], where M 0 = 0.35 can be 
achieved. 

9.3.2 General Aviation Propeller 

Use of a general aviation propeller in conjunction with a 
model fuselage has the advantage that the noise mechanisms 
show certain similarities with those of the propfan. However, 
it is anticipated that there will be several problems associated 
with a test configuration of this type, one problem being the 
general inflexibility in matching the model characteristics 
with those for the baseline airplanes. 

If a fullscale propeller is used with a model fuselage, the 
combination of propeller and fuselage would be too large for 
any wind tunnel test section except the large, low speed 
tunnels which could not reproduce a forward Mach number of 0.8. 
Consequently, the tests would have to be performed under static 
conditions . 

The longitudinal spatial variations of external sound pressure 
level predicted for the three study aircraft can be compared 
with measured variations for general aviation aircraft, as 
shown in Figure 76. It is seen that for a given tip clearance, 
r/D, the rate of change of sound level predicted for the propfan 
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0.2 0.4 

Propfan Calculations 


General Aviation [9] 
Model Tests [46] 


Harmonic 

1 2 3 

Static 
30 Knots 

Hi • A 

O O o 


Propfan [ 6 3 
r/D = 1.2 



r/D = 0.4 


= 0.6 r/D = 0.051 


•0 M f «0.75 r/D = 0.083 
-VM = 0.60 r/D = 0.333 


FIGURE 76. 


SPATIAL VARIATION OF PROPELLER OVERALL 
SOUND LEVEL IN LONGITUDINAL DIRECTION: 
COMPARISON OF PROPFAN PREDICTIONS AND 
GENERAL AVIATION MEASUREMENTS 


'' ) 


- 232 - 



under cruise conditions is much greater than that for the (model 
test) general aviation propeller under static conditions. 

A comparison of pressure field trace velocities also shows 
significant differences between the propfan predictions and 
general aviation measurements (Figure 77), although in this 
case only one set of data is available for general aviation 
propellers. At large separation distances the differences 
between trace velocities for propfan and general aviation pro- 
pellers are associated mainly with the contribution from 
forward motion, whereas at small separation distances the 
controlling parameter is tip clearance. The circumferential 
trace velocity is also a function of tip clearance. 

From this brief comparison it is seen that a general aviation 
propeller under, static conditions cannot reproduce simultaneously 
both sound level spatial decay and pressure field trace velocity - 
two parameters which are considered to be important in the 
analytical model. 

Another problem is associated with the values of the blade 
passage frequencies. If the rotational tip Mach number is 
maintained at 0.8 for the general aviation propeller then the 
blade passage frequency will be too low (because of the small 
number of blades). Consequently the required scaled fre- 
quencies for the propfan blade harmonics will have to be 
modeled by higher-order harmonics of the general aviation 
propeller. Since harmonic level decreases as harmonic order 
increases, the higher order harmonics will be lower in level 
and could create a signal-to-noise ratio problem for the high 
transmission loss tests. The higher order harmonics may also 
have different directivity patterns. 
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This brief discussion has identified several important differ- 
ences between the pressure field of a fullscale general aviation 
propeller and the desired field for the model propfan. These 
differences are sufficiently important that the general aviation 
propeller is not recommended as a model source for direct model- 
ing of the baseline aircraft conditions. It could be used, 
however, as a noise source in a test to validate the analytical 
model for general aviation conditions. 

9.3.3 Electro-acoustic source 

An electro-acoustic source, such as an electro-magnetic driver 
coupled to an acoustic horn offers great flexibility in fre- 
quency content of the excitation. Also, by suitable choice of 
horn characteristics and orientation, the source can be used 
to represent pressure field trace velocities in longitudinal 
and circumferential directions, and sound level spatial varia- 
tion over the model surface. 

Figure 78 presents the directivity required of an electro- 
acoustic source used to simulate the propfan directivity, for 
each of the study airplanes. The Hamilton Standard method 
[6] was used to construct the curve. Comparison with the 
directivity of a point source, located to develop the required 
pressure field trace velocities of the model fuselage surface 
shows that an omni-directional source is not adequate. Strong 
longitudinal directivity is required. 

Ray acoustics has been used to determine the propfan trace 
wavenumbers (see Section 4.6) and is also used to determine 
the position of the electro-acoustic source. Figure 79 
presents a polar plot of the required directivity for the 
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electro-acoustic source for inboard-propeller simulation for 
each study airplane. In order to provide representative 
dimensions, the figure also shows the location of a narrow-body 
panel. It is clear that the electro-acoustic source must have 
pronounced longitudinal directivity. 

Figure 80 shows the predicted variation in circumferential levels 
for an electro-acoustic source, assumed to have point-source 
radiation characteristics in the circumferential direction, 

(X is the center of. the electro-acoustic source). The propfan 
acoustic decay rate is much smaller in the circumferential 
direction than in the longitudinal direction, and is well 
approximated by a point source. Thus the electro-acoustic 
source should have little circumferential directivity. 

A long, flat horn, aligned so that the horn axis is perpendi- 
cular to the model surface, will produce such a pressure field 
spatial variation. The horn width will be approximately two 
acoustic wavelengths long at the lowest model test frequency, 
and the depth will be quite narrow. Possibly two horns will be 
required depending on the frequency range to be tested. An 
alternative configuration would involve a line array of several 
closely-spaced point sources, aligned parallel to the model 
fuselage axis. An electro-acoustic source (or, for that matter 
a general aviation propeller) under conditions of zero forward 
motion could not model the subsonic trace velocity predicted 
for the fuselage region well forward of the propeller plane of 
rotation. However, calculations using the analytical model 
indicate that this forward region is of low importance in terms 
of noise transmission and can be omitted from the validation 
process for the analytical model. 
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9.4 


Validation Experiments with a Model Fuselage 


9.4.1 General 

Sections 9.2 and 9.3 presented an investigation of the excita- 
tion, structure and cabin parameters affecting the design of 
experiments for validation of the analytical model and associ- 
ated noise control predictions. It was assumed in the investi- 
gation that the performer of the validation tests would not 
have ready access to the computer program developed for the 
analytical model. Thus, the design of the model and experiment 
would be carried out in a manner such that the prop— fan air- 
plane conditions could be adequately simulated. Analytical 
results available from the present study could then be used 
directly for comparison with the experimental data. 

It is recognized that the separation of experimenter and 
analyst imposes severe constraints on the validation process. 

In contrast, verification of a similar model for the payload 
bay acoustic environment of the Space Shuttle orbiter vehicle 
[4,10] was achieved by close collaboration between the ex- 
perimenter and analyst. With the Space Shuttle experience 
in mind, care has been taken in the planning of the proposed 
test program to minimize problems associated with a possible 
lack of access to the computer program for the analytical 
model . 

In this section a detailed description is presented for a 
model test program that would provide experimental validation 
of the analytical model and its use in the prediction of prop- 
fan interior noise control. Problems discussed include the 
test configuration (size and number of test models, test en- 
vironment), excitation details (frequency content-, orientation, 
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spectrum shape) , error assessment as it relates to comparisons 
between experiment and predictions, required measurements, and 
data acquisition methods. 

9.4.2 Model Test Configurations 

The most straight forward test program would be the validation 
of the analytical model only for the case of add-on noise 
control treatments. This could probably be accomplished with 
the use of only one structural model. ^ On the other hand, if 
the predictions are to be validated for both add-on and 
advanced noise control methods , several structural models will 
be required, with at least two incorporating the advanced 
noise control concepts. 

The requirements for accurate simulation of the noise trans- 
mission into a prop-fan airplane under cruise conditions have 
been identified in section 9.2.4 and in section 9-3. While 
practical experimental design will force various compromises 
to be made, it is argued that good simulation of full-scale 
fuselages representative of the baseline airplanes is feasible 
and that the use of an electro-acoustic loudspeaker source 
would enable adequate simulation of the major characteristics 
of prop-fan acoustics in the region of maximum transmission 
of acoustic power through the airplane fuselage. 

Figure 81 shows a sketch of the end view of the recommended 
test configuration for the evaluation of the analytical model 
for the baseline structure and the add-on noise control 
treatments. The model consists of a complete cylinder which 
is lined with appropriate sidewall treatment. The cylinder 
is elevated on supports and located in an essentially anechoic 
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environment. The acoustic source, a loudspeaker, is located 
approximately 1.5 model fuselage diameters from the skin 
surface and part way between the cylinder ends. The structure 
of the cylinder would be of skin-stringer-frame construction 
using aluminum material. The dimensions of the structural 
items would be scaled with respect to cylinder diameter as 
indicated in sections 9*2 and 9.3« 

The diameter of the model fuselage is approximately half that 
for the narrow-body study airplane, i.e. 2m. Such a scaling 
allows for convenient measurement of the various important 
physical quantitites and does not introduce significant con- 
struction difficulties which, if the test model is built to 
too small a scale, might change the physics of the noise 
transmission. The length of the model fuselage is selected 
to be between 2 and 3 model diameters, to ensure adequate 
decay of the pressure amplitude over the model length and to 
allow the peak excitation amplitude to be located away from 
the model boundaries. The model length, while not scaled in 
proportion to the fuselage radius, is also required to give 
sufficient modal overlap for the acoustic modes. Care is 
required to position the interior absorptive material to simu- 
late a longer cabin volume: additional absorption will be 
required on the test model forward and aft bulkheads. End 
caps will be required on the test model and the caps will have 
to have high transmission loss characteristics so that there 
are no flanking paths when the high transmission loss treat- 
merits are tested on the sidewall. The tests proposed will not 
require pressurization of the model. Thus the end caps will 
not have to withstand pressurization loads. The caps should 
be removable to allow for introduction of the various noise 
control treatments to the model cabin space. No simulation 
of the cabin furnishings is considered necessary. 
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In order to achieve the correct spatial variation for the 
exterior pressure field, the noise source will be located at 
some position between the two planes defined by the end caps, 
and the directivity of the source will be such that there is 
a propagating pressure field over the model. It should not 
be necessary to provide extension sections to, the test cylin- 
der to correct for diffraction at the ends of the structure. 
However, the test chamber should be treated with absorptive 
material so that the propagating pressure field over the 
model surface is not contaminated unduly-- by any reverberant 
field in the chamber. 

Add-on treatments tested in the model cylinder would concen- 
trate on double-wall systems with porous material in the 
space between the walls. All the parameters of Interest would 
be scaled as required from the analysis in Sections 9*2 and 
9.3. Special consideration should be given to the identifi- 
cation and elimination of flanking paths. Other noise control 
treatments could be tested although only the double-wall was 
found, analytically, to achieve the noise reductions required 
to meet the 80 dB goal for A-weighted Interior sound level. 

It is possible, in principle, to construct scale models for 
any of the advanced noise control methods although, in prac- 
tice, there may be difficulties in achieving the correct 
scaling for the more exotic designs. Structures which involve 
modifications to conventional aluminum skin-stringer-frame 
construction should follow the same scaling laws as for the 
baseline structure. Composite structures will probably need 
some additional analysis to determine the most appropriate 
scaling procedure. 
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The specific choice of test hardware for the advanced concepts 
is less well defined than in the case of add-on treatments. How 
ever, potential candidates include those discussed in Section 8. 

a) Additional ring frames on the baseline model 

b) Frames with increased stiffness achieved by 
redesign of the frame cross-section 

c) Increased skin thickness 

d) Extensive damping treatment 

e) Composite (honeycomb) skin with ring frames, (with 
and without additional damping) . 

It should be noted that the noise control concepts listed above 
have been investigated using the analytical model and in no 
case did the predicted noise reduction approach the required 
reduction of 25 dB for the A-weighted level. Thus validation 
of the analytical model in this case has the somewhat negative 
result of confirming that the noise control approach is not 
very effective. 

The model test structure recommended for validation of the 
analytical model predictions for the baseline structures and 
the add-on noise control treatments is in the form of a 
cylinder with end plates. However, this type of model could 
be extremely expensive for tests of advanced type structures, 
unless these structures could be formed as modifications to 
the basic skin-stringer-frame cylinder. An alternative 
approach would be the use of large curved panels rather than 
complete cylinders. These panels would be placed in one wall 
of an enclosure, which could be constructed as a wooden double- 
wall with high transmission loss material in the space between 
the walls. The cost of construction for the advanced struc- 
tures would then be reduced. 
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The use of test items should not have a large impact on the 
validation of the analytical model since it is assumed when 
making the noise transmission predictions in Sections 6 through 
8 that the fuselage structures can be represented as curved 
panels rather than a complete cylinder. The use of curved 
panels in the test program would mean that the above assumption 
could not be verified, but all other assumptions in the analy- 
tical model could be validated. 

9.4.3 Basic Measurements 

Two sets of measurements should be made for the validation 
tests. One set, basic measurements, represents the minimum 
requirements for validation. The second set, diagnostic 
measurements, will provide the additional information required 
to fully understand the reliability and accuracy of the analy- 
tical model. 

The basic measurements involve the exterior and interior 
sound levels. The exterior measurements ensure that the 
correct excitation is being used and the interior measurements 
indicate whether or not the 80 dB(A) goal is achieved. Exterior 
data required to describe the excitation field will include 
the variation of sound level in the longitudinal and circum- 
ferential directions, and the coherence and phase of the 
pressure field in the two directions. 

In the case of the interior sound levels, the following 
measurement procedure is recommended: 

a) Determine space average sound levels from a series 
of fixed microphones or from microphone traverses. Micro- 
phone locations can be associated with equal volumes for 
ease in space averaging. 
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b) Measure sound level gradients in radial and axial 
directions to obtain estimates of spatial variation. 

c) For given microphone locations, measure sound levels 
as excitation frequency changes by small amounts. Calcu- 
late deviation for sound levels. 

The recommendation to measure both space-average and spatial 
gradients for the interior sound levels is based on experience 
from other validation programs. Factors to be considered 
include : 

a) the model test cavity may not be fully representative 
of the airplane cabin shape and contents; 

b) a detailed analytical representation of the cavity is 
not possible within the accuracy of the present analy- 
tical model; and 

c) the analytical model predicts space-averaged levels 
and then estimates the distribution along the cabin using 
imaging methods. 

Space-averaging has the advantage that inaccuracies in predic- 
tions of the detailed modal response of the cavity are 
smeared-out. However, even space-averaging has to be performed 
keeping in mind that the measured value has, for practical 
reasons, to be based on a finite number of measurement loca- 
tions. Statistical analysis has to be performed to relate 
the measured value to the expected value. This can be done, 
for example, by means of standard deviation and confidence 
limits . 
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9.4.4 Diagnostic Measurements 

It frequently happens that, when experimental data are com- 
pared with an analytical model, deficiencies are found in 
the model. Improvements will then be required to improve the 
accuracy of the model, but usually these improvements can be 
made only if there are adequate diagnostic measurements to 
supplement the basic data. It is recommended that diagnostic 
measurements be made during the validation tests. These addi- 
tional measurements will relate directly to various outputs 
of the analytical model so that appropriate corrections can 
be made to the model. 

Measurements recommended for this diagnostic phase of the 
validation tests include the following items, which describe 
the structural response, and structural and acoustic damping. 

•Structural response: 

Spatial variation of model vibration 

Flexural wavelengths (by means of coherence measurements) 
Modal density (by means of point force excitation) 
Standard deviation of response as excitation fre- 
quency is varied slightly 

• Structural Damping: 

Total loss factors for the structure (by means of 
vibration decay measurements) 

Radiation loss factors (by relating vibration 
response to radiated acoustic power 

• Acoustic Damping: 

Acoustic loss factors (from acoustic decay 
measurements) 
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9.4.5 Performance of the Validation Tests 

In terms of the basic validation tests, the experimenter 
could simulate an exterior spectrum shape and then measure 
the A-weighted sound level in the model cabin. However, 
the resulting overall noise reduction would depend on the 
selected spectrum shape, and a coarse test of this type is 
not recommended for the validation of the analytical model. 
Instead, it is recommended that the tests be conducted on an 
individual frequency (harmonic) basis.. The noise levels or 
noise reductions would be measured, and compared with the 
analytical predictions on a frequency by frequency basis. 

Only then can the analytical model be subjected to an appro- 
priate validation process. 

This individual harmonic approach has several other advantages 
Firstly, the tests can be performed and the model validated 
without concern as to the spectrum shape which will be gener- 
ated by the final propfan design in actual cruise conditions. 
Secondly, if a noise source is used which enables the sound 
level to be adjusted at a given frequency, then signal-to- 
noise ratio problems at higher frequencies will be minimized. 
The proposed electro-acoustic source provides this control of 
the excitation characteristics. 

The proposed tests include the following excitation conditions 

(i) Perform the tests at the model frequencies of 
interest for the three baseline aircraft used in the 
analysis . 

(ii) Repeat at several frequencies either side of each 
frequency used in (i) in order to obtain a measure of 
the statistical variation about the mean. The number of 
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frequencies used will be determined by the measured stan- 
dard deviation and the desired confidence limits. (.Various 
tests on airplane fuselages have shown that the noise re- 
duction characteristics can change rapidly as frequency 
changes.) The statistical analysis is discussed briefly 

in Section 9.^.7. 

(iii) Narrow-band noise, centered at given harmonic fre- 
quencies to be used for comparison with (i) and (ii). 

The first step in the validation tests will be to ensure that 
the excitation pressure field satisfied the model requirements 
in terms of spatial variation in sound level, coherence and 
phase (or trace) velocity. Subsequent steps in the test 
program would then be concerned with measurement of interior 
sound levels or structural response — with adequate care 
being taken to ensure that the excitation pressure field 
characteristics remain unchanged. 

9.4.6 Instrumentation Requirements 

The instrumentation requirements for signal generation, data 
acquisition and data reduction are associated with the use 
of discrete frequency and narrowband noise. Typical schematic 
diagrams are shown in Figures 82 and 83, where both on-line 
and recorded data options are included. The use of on-line 
data acquisition and reduction will probably be the selected 
approach for many of the tests performed for the validation. 

Figure 82 shows potential signal generation and data acquisi- 
tion systems. The excitation signal is generated as a dis- 
crete frequency or narrowband noise and used as input to a 
loudspeaker. The signal can also be recorded for future 
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reference. In the initial set-up tests, when the loudspeaker 
is being located so that it generates the required noise 
field, flush-mounted microphones would be placed on the test 
structure to measure pressure amplitude and cross-spectral 
density. The microphones would be arranged in a cruciform 
array, with 12 to 20 locations selected to adequately define 
the pressure field characteristics. A reference microphone 
would be used to monitor sound levels from test to test. 

These microphones would be, typically, 6mm (0.25 inch) in 

diameter. 

The microphone array inside the model fuselage would consist 
of 8 to 12 microphones positioned, for example, at fixed 
positions in the model cross-section, but capable of being 
moved in the longitudinal direction. The number of measured 
axes would be determined from an initial examination of the 
measured data and from the required measurement accuracy. 
Locations in the cross-sectional plane would be selected so 
that the areas associated with each microphone were equal. 

The structural response would be measured by an array of 
accelerometers attached to the interior surface of the model 
fuselage shell. These accelerometers would measure the 
spatial characteristics of the fuselage vibration. 

A typical schematic for data reduction is shown m Figure 83. 
Spectral analysis is straightforward since the basic filter- 
ing is performed by the signal generation system. The analysis 
consists of measuring the level of the discrete frequency or 
narrowband acoustic signal for the excitation of the interior 
pressure fields. Cross-spectral data reduction requires, for 
example, a- multi-function digital signal processor. Signal 
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FIGURE 83. SCHEMATIC DIAGRAM OF DATA REDUCTION SYSTEM 
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decay measurements, to be used in the estimation of structural 
and acoustic loss factors, requires the use of a decay meter 
(or a level recorder if the decay is not too rapid). 

9.4.7 Data Reliability 

Two factors have to be considered under the heading of data 
reliability. One factor, data repeatability, is important 
because of the use of discrete frequency excitation. The 
second factor, data accuracy, is concerned with the comparison 
between measured and predicted results. 

Because of the use of discrete frequency excitation, micro- 
phone location may be critical in achieving data repeatability 
when there is a small number of modes. This problem is allevi- 
ated by performing space-averaging of the measured sound levels 
in the cabin and by the fact that, in the frequency range of 
interest, there is high modal overlap for the acoustic modes. 
The situation will become somewhat more severe when spatial 
gradients are measured for the interior sound field at a given 
frequency of excitation. Thus, repeat measurements should be 
made for a given test condition, and allowances made when 
comparing spatial gradients for different sidewall configura- 
tions . 

The proposed tests include discrete frequency excitation at 
several frequencies on either side of a given frequency of 
interest. One purpose of this variation in frequency is to 
allow for uncertainties in the prediction of structural and 
cavity resonance frequencies. By the same token, the data, 
when averaged over frequency, will show improved repeatability 

from test to test. 
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The comparison between measured and predicted sound levels or 
noise reductions will have to be made on a statistical basis. 

The importance of such an approach was evident in the valida- 
tion of the Space Shuttle payload bay acoustic model [11], a 
situation for which the excitation was broadband in nature and 
fairly uniform in level. Instead of comparing only the average 
measured and calculated noise levels, consideration was also 
given to the variability of the measurement used to make the 
comparison. Since both the exterior and interior average 
sound levels were based on a finite number of sample measure- 
ments, it followed that the computed estimates would be 
subject to random errors. These errors were taken into con- 
sideration by computing the variances for the measured levels 
and then constructing confidence limits such as those shown 
in Figure 84. A similar approach should be followed for the 
propfan validation tests. 

9 . 5 Aircraft Tests; 

9.5.1 Applicability of Aircraft Tests 

The discussion in Sections 9.2 and 9.3 has concentrated on the 
design of scale model tests which could be used to validate 
results from the analytical model. The tests were envisaged 
as providing answers to the questions posed in Section 9.1. 

An alternative approach is that of using an existing airplane 
and performing either ground or flight tests to measure the 
noise transmission into the fuselage. The question now arises - 
- can the use of an existing propeller-driven airplane provide 
data which will validate the analytical predictions for the 
propfan-powered study airplanes? If the answer to the question 
is negative, then airplane tests can be discarded unless it 
can be shown that they would serve an alternative, but still 
useful, evaluation of the analytical model. 
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DIFFERENCES BETWEEN PREDICTED AND 
MEASURED PAYLOAD BAY LEVELS FOR 
SPACE SHUTTLE VALIDATION TESTS [| |] 






In order to answer the aboye question regarding the direct 
applicability of airplane tests, consider first the relevant 
characteristics of the three aircraft considered in the analyti— 
cal model development. The aircraft are low-winged and have 
two or four propellers, with the plane of rotation for all the 
propellers intersecting the passenger cabin sidewall. The 
propeller diameter is equal to, or greater than the corres- 
ponding fuselage diameter, the propeller has 10 blades and a 
tip rotational Mach number of 0.8, the separation distance 
between propeller tip and fuselage structure is greater than 
0.3D where D is the propeller diameter', the fuselage diameter 
is greater than 2.44m (8 feet), and the airplane Mach number 
is 0.8. 

Obviously there is no airplane at the present time which meets 
all these requirements, although there are some aircraft which 
meet some of the requirements. Two such airplanes are the 
Lockheed Electra/P3 and the Swearingen Metroliner. They are 
both low-wing, turboprop aircraft; other typical characteris- 
tics are listed in Table 16. 


Table 16. Comparison of Study and Existing Aircraft 


Aircraft 

Propfan 

Baselines 

Lockheed 

Electra/P3 

Swearingen 

Metroliner 

Number of Propellers 

2-4 

4 

2 

Number of Blades 

10 

4 

3 

Propeller diameter D (m) 

3.37-5.83 

4.11 

2.49 

Fuselage diameter cf> (m) 

2.44-6.02 

3.45 

1.73 

Tip Clearance (Inbd.Prop.) 

0.3D-1.3D 

0. 2D 

0 .12D 

D/ <(> 

0.97-1.38 

1.20 

1.44 

Tip Rotational Mach No. 

0.8 

0.7 

0.8 

Airplane Mach No. 

0.8 

0 ..6 

0.4 

Typical Cruise Alt. (ni) 

9140 

7620 

6100 

Blade Passage Freq. (Hz) 

133-231 

68 

96 
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The main differences between the baseline aircraft and the 
Electra and Metroliner occur in the number of blades, tip 
clearance and airplane Mach number. The number of blades (and 
propeller rpm) Influences the blade passage frequency, tip 
clearance influences the spatial decay and trace velocities 
of the external pressure field, and airplane Mach numbers in- 
fluence the pressure field longitudinal velocity. All these 
parameters have been shown to be important in the analytical 
model in determining the noise reduction from exterior to in- 
terior of the fuselage. 

Differences in blade passage frequency can be off-set to some 
extent by considering the higher order harmonics for the 
Electra or Metroliner. For example, if one selects as a 
criterion, a constant value of k<j> where k is the acoustic wave 
number of the pressure field and $ is the fuselage diameter, 
then the blade passage frequency for a propfan airplane is 
represented approximately by the third order harmonic of the 
Electra or Metroliner propeller. However, the use of higher 
order harmonics runs into the problem of signal-to-noise ratio 
for the tests on high-transmission loss sidewalls, since the 
external levels of the higher order harmonics for the subsonic 
propellers in the Electra and Metroliner will be relatively 

low. 

Trace velocities in the longitudinal direction are influenced 
by tip clearance and airplane velocity. The effects of de- 
creased tip clearance and decreased airplane velocity both act 
to reduce the trace velocity in the region aft of the plane of 
rotation of the propeller. In addition, using Aero Commander 
data [9] as a basis, the decreased tip clearance also reduces 
the circumferential trace velocity. This velocity may be 
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subsonic for the Electra and Metroliner whereas the analytical 
model predicts a supersonic circumferential trace velocity. 

Finally, as tip clearance decreases, the longitudinal decay 
of the pressure field becomes more rapid. Thus, if one takes 
the directivity patterns to be those predicted by Hamilton 
Standard the longitudinal decay factors associated with the 
tip clearance of the Electra and Metroliner aircraft will be 
at least twice as large as those for the corresponding propfan 
designs. It is also reasonable to exp-ect that the small tip 
clearances will increase the spatial decay in the circumferen- 
tial direction. Typical calculations for propfan aircraft 
indicate that a 50 % increase in the spatial decay factors can 
produce about a 5dB change in joint acceptance. Thus the 
decay parameters are important in determining noise trans- 
mission and, hence, the spatial extent of noise control treat- 
ment . 

From the above discussion it is apparent that significant 
acoustic differences exist between the baseline propfan air- 
craft on one hand and present day propeller aircraft such as 
the Electra and Metroliner on the other hand. Consequently, 
flight tests in current aircraft would have to be viewed as 
validating the model for the specific test airplane, rather 
than verifying the predictions for the study aircraft. Cor- 
responding calculations would have to be performed to apply 
the analytical model to the particular test airplane. 

Some comment is also appropriate with regard to the structure 
of the airplane fuselage, particularly that of a small-diameter 
airplane such as the Metroliner. Analysis of the three study 
airplanes indicates that the small— diameter fuselage skin is 
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thicker than it would be for a similar diameter scale model of 
the wide or narrowbody fuselage. Consequently a small-diameter 
fuselage does not form a completely representative scale model 
of a larger diameter fuselage. Differences in calculated 
radiation efficiency, modal density and band-averaged joint 
acceptance are shown in Figures 85 and 86. The situation 
becomes worse if the smeared-out mass of the small diameter 
fuselage is reduced to that of the scale model by removal of 
the stringers. Thus, if a small-diameter fuselage is used as 
the test hardware, the validation cannot be considered as a 
direct verification of the analytical model predictions for 
the narrow or wide body fuselages. Some adjustments will be; 
necessary, for example, by use of figures such as Figure 86. 

There are other disadvantages in using a current airplane as 
a means of validating the analytical model. The airplane 
provides little or no flexibility in varying excitation fre- 
quency, pressure decay rate or trace velocity. It would be 
highly desirable to have such flexibility in order to check 
the analytical model over a range of conditions. Also, while 
a test program in an actual airplane is highly desirable as 
a means of reducing interior noise levels which are trans- 
mitted by a variety of paths from a variety of sources, it is 
not desirable as a means of validating an analytical model 
directed to a particular noise source and a particular noise 
transmission path. It is recognized that the propf.an air- 
plane may well have numerous transmission paths, to which 
noise control methods must be applied, but such noise control 
approaches should be the subject of separate programs, or 
await testing on a prototype propfan aircraft. 
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9.5.2 Airplane Tests 


Although It Is concluded above that airplane tests would not 
validate the analytical predictions for the study aircraft, some 
useful Information could be obtained from such tests. Thus It 
Is appropriate to review briefly the type of tests which could 
be performed and the measurement procedures which should be 
followed. 

The simplest test would be one In which several noise control 
treatments were installed in a propeller-driven airplane and 
the noise reduction measured for each treatment. Such a test 
gives very little information on the validity of the present 
analytical model. A more useful experimental program would 
provide information on the exterior pressure field, structural 
response, and cabin noise levels for a test airplane, and the 
noise reductions due to add-on treatments. Data appropriate 
to these test conditions would be used in the analytical model 
to predict the structural vibration and cabin noise levels. 

Then the experimental and analytical results would be compared 
and necessary modifications, if any, made to the analytical 
model. 

As the tests are concerned with propeller noise only, the test 
airplane should be powered by turboprop engines, and the plane 
of rotation of the propellers should intersect the passenger 
cabin. It might at first seem appropriate to install bulkheads 
in the cabin so that the noise transmission region can be con- 
trolled to a specific region of the fuselage structure. There 
are, however, restrictions on this approach because the analy- 
tical model assumes high modal overlap for the acoustic modes 
for all frequency bands of interest. Consequently the size of 
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the receiving cavity cannot be made too small. This assump- 
tion of high modal overlap also influences the choice of 
propeller harmonics for the noise tests. It is probable that 
the lowest order harmonics of the propeller will occur at 
frequencies which do not satisfy the condition of high acoustic 
modal overlap. The measurements will then have to be per- 
formed for the higher order harmonics. High transmission loss 
treatments should be applied to the floor. 

It is known that the effect of forward motion on the noise of 
subsonic propellers is to reduce the levels of the higher 
order harmonics, a trend which may have an adverse effect on 
the signal-to-noise ratio when treatments with a high trans- 
mission loss are tested. One possible way of alleviating the 
problem is to perform the tests under static conditions. This 
has the added advantage that the tests are more convenient to 
perform but there is the disadvantage that forward motion is 
not simulated in the trace velocity of the exterior pressure 
field. Resolution of these arguments would probably need an 
evaluation of the specific airplane proposed for the tests. 

Assuming that the objective of the test program is a validation 
of the analytical model as applied to the specific test air- 
plane, the following measurements should be made. 

(i) Exterior pressure field: 

Use a cruciform array of flush-mounted pressure 
transducers to measure the spatial decay, coherence, 
and phase angle parameters for the pressure field 
in the longitudinal and circumferential directions. 

The array would be similar to that used in the Aero 
Commander tests [8], but more extensive in terms of 



number of transducers and spatial extent. The 
results of these measurements would be used as a 
basis for the representation of the exterior 
pressure field in the analytical model; 

(ii) Structural vibration: 

Analysis of the study airplanes indicates that the 
acoustic power flow is dominated by resonant re- 
sponse. Thus vibration measurements — which give 
the response — can be used for diagnostic purposes. 

The vibration measurements would be made using 
light-weight accelerometers mounted mainly on fuse- 
lage skin panels in the neighborhood of high ex- 
terior noise levels. 

(iii) Cabin Noise Levels: 

Noise levels would be measured in the cabin in the 
region of high exterior sound levels. The measure- 
ments would be repeated for a) bare interior (with 
some absorption), b) baseline sidewall treatment and, 
c) one or more high-transmission-loss treatments. 

The measurement locations would be chosen such that 
space-averaged levels could be calculated for an 
appropriate cabin volume, and such that longitudinal 
and radial gradients could be determined. 

The test conditions would be limited by the performance of the 
selected airplane. Assuming that flight tests were performed, 
the airplane would be operated at high speed with high pro- 
peller rpm. Tests would be repeated at different engine rpm, 
provided that rpm variations can be accomplished. To avoid 
beat effects, the right and left hand propellers should be 
operated at different rpm. Measurements could be performed 
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with and without cabin pressure differential although the 
analysis indicates that pressure differential is not very 
important for the frequency range which dominates the noise 
transmission through the sidewall. Consideration could be 
given to performing certain ground, static tests for compari- 
son with the flight results. 



10.0 CONCLUSIONS 


An analytical model has been developed to predict the sound 
levels in the interior of high-speed propeller-driven aircraft. 
Although the assumptions contained in the model have not yet been 
validated experimentally, the analytical model has been used to 
predict the sound levels in three study aircraft and to explore 
potential noise control concepts. 

The analytical study indicates that, if the high-speed propeller- 
driven aircraft have fuselage structures and cabin sidewall 
treatments similar to those in current-day turbofan-powered air- 
craft, the cabin noise levels will reach a maximum of about 105 
dB(A). These levels are about 25 dB higher than the goal of 80 
dB(A) . Consequently, a significant increase in cabin sidewall 
transmission loss is required, and this increase has to be 
achieved at relatively low frequencies (130-230 Hz) since the 
cabin A-weighted sound levels are dominated by contributions at 
the propeller blade passage frequency. 

The noise control methods investigated in the analytical study 
include add-on and advanced concepts. The add-on treatments do 
not involve changes to the fuselage primary structure whereas 
the advanced concepts usually do require, structural changes. Of 
the concepts studied, only one shows the capability of achieving 
the required reductions in transmission loss. This treatment 
involves the use of a double— wall system with a heavy , limp 
trim panel which is isolated from fuselage vibrations. All the 
other concepts studied provided relatively small noise reduc- 
tions. The weight penalties associated with the double— wall 
system are significant, ranging from 0.75$ to 1.5 % of take-off 
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gross weight, depending on fuselage diameters. In terms of 
passenger payload, the penalties lie in the range of 4 to 19% 
of payload weight. In both cases the highest penalties are 
associated with the smallest airplane. 

It should be noted that an extensive weight-optimization study 
was not performed and that, by suitable combination of different 
noise control treatments, the weight penalties could be reduced. 
On the other hand, it should be noted that the achievement of 
the required large transmission losses in a practical situation 
has yet to be proven, and the effects of windows and other po- 
tential flanking paths have been excluded from the present 
analytical model. 

Experimental validation of the analytical model and the predicted 
noise reductions is an essential part of the overall development 
of the model. Such validation should be possible with laboratory 
tests, and experimental procedures for performing these tests are 
presented in this report. 
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APPENDIX A 

JOINT ACCEPTANCE EXPRESSIONS FOR NON-HOMOGENEOUS PRESSURE FIELDS 

The joint acceptance function describes the coupling between the 
excitation field and the structure and can be defined for an exci- 
tation tonal frequency as 




A2 


< Pxo ( “ )> 


II 


X X' 


C ( x , x ' 

p 


;u))i|; r ( x)ip r ( x’ )dx dx' 


( A . 1 ) 


where A is the panel area, <p* Q (w)> is the mean square pressure 
at Xo attributable to the b th harmonic, C (x,x f ;w) is the co- 
spectrum of the non-homogeneous pressure field (defined relative 
to Xo) , and ^ r (x) is the shape of the r th mode. To evaluate this 
joint acceptance function, expressions are required for the mode 
shapes of the structure and for the excitation pressure field. 


The excitation is represented as a convecting pressure field with 
pressure amplitude decaying with distance from the location of 
peak levels Xo and with spatially-decaying coherence. Exponential 
forms for both amplitude and coherence decays have been modeled 
to enable convenient mathematical manipulation. Then 


■a I x-Xo I -a I x* - Xo! 


C ( x,x’ : w) = e 

p 5 ’ 


-C ! X* -X 


1 cos [k(w)( x* - x)]<pj 0 (w)> 

(A. 2) 

where a is the amplitude decay rate, c is the coherence decay rate, 
and k is the excitation wavenumber at frequency w. 


Consider the example of a one-dimensional system in which a beam 
of length L is driven to vibrate in itsm th mode shape by a non- 
homogeneous pressure field. In general, Xo , the location of peak 
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pressure levels, will lie at some arbitrary point along the beam. 
The expression for the joint acceptance function has been evaluated 
for such a case; but for present demonstration purposes a simpler 
case will be considered, where the maximum Intensity Is located at 
Xo = 0. Figure A.l describes this situation. 



FIGURE A.l EXCITATION OF BEAM MODE BY 
NON-HOMOGENEOUS PRESSURE 
FIELD 

With Xo = 0, the joint acceptance (one-dimensional) becomes 


j * ( w ) 
mm 



cos[k(w)£] ^ m (x).^ m (x' )dxdx* 

(A. 3) 


where £ = x ' - x. 


Using sinusoidal mode shapes, j z (w) becomes 

’ "mm 


J" (w) = 


mm 


■ 


-ax . , , 

e sin k x dx 
m 




e ax e c I 5 1 C os(kC) sin k x* dx’ 

m 


(A. 4) 


where k Is the wavenumber of the structural mode, 
m 
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Evaluation of j 2 (w) leads to the result: 
mm 


A’ / TNin-bL ' 

■ TE2 k- 1 ’ e L 


•b sin kL - (k + k)cos kL b sin kL - (k - k)cos kL 


+ (k + k)/b + + (k -k)/b“ 
v m m 
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where: 


TtfTlV 2a 
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x / _ -2aL x k (l - e" 2aL ) 
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k + k 
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k 





(k m + k ) . (k m ~ k) 
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(_D m + 1 e -bL[- b sln kL + 

, 1 >.m+l -bL r . , T 

(-1) e [b cos kL - 

(_l) m+ -*- e ”kk[b sin - 

(-1) e [b cos kL + 


(k + k)cos kL] 

(k + k)sin kL] 
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m 

(k - k)sin kL] 



Special cases: 

a) Zero amplitude and correlation decay rates: a = c = 0 

2 [1 - (-l) m cos kL] 

J 'mm (w) = [1 - (k TiT^T— 1 ' 

which is the well-known result, corresponding to plane acoustic 
waves propagating with zero spatial decay along a beam of length 
L [A.l] . 

At coincidence where k = k , 

i 2 (co) = 1/4 for all modes. 
u mm 

b) Zero correlation decay ,, but a is finite: c = 0, a / 0. 


n 


~) 




1 


j 2 (u) = [A" 2 + B" 2 ]/4L 2 
J mm 


(A. 6) 


where: 
A" J 
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The action of strong amplitude decay is to dampen the oscillation 
° f J *mm^ as k//k m varles through coincident conditions (see Pig. 2). 


Asymptotic limits can be developed for the case where strong 
amplitude decay occurs , i.e. when e a ^<< 1. Three examples are 
given below: 

For k m >> k and k L > aL, 


j 




(A. 7) 


For k << k and k L > aL, 
m m i 



(A. 8) 


For k m - k (i.e. close to coincidence) . and kL >> aL, 


l 2 

WiJ (A. 9) 

where l is that equivalent length over which the excitation may be 
considered to be of equivalent constant amplitude. % is analogous 
to a correlation length for homogeneous excitations in which the 
coherence decay is rapid relative to the beam length. 


j 2 ( w ) 
mm ' 


1 

4a 2 L 2 
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EXPECTED VALUE OF THE MODAL ADMITTANCE FUNCTION g(w b »w n .) 


It is required to solve for the expected value of s(w b >“ n )> which 
is the expected value that g(w b ,w n ) takes in the frequency band 
Aw as, for example, w fe varies within Aw. 


By definition [B.l]., 


E[g(w b ,w n ) ] 



■ (w, ,w ) 
1 b ’ n 


p ^“b ,a) n^ du b 


dw 

n 


(B.l) 


where w fc and w n are statistically independent processes, in which 
case Eq. (B.l.) reduces to 


E[g(w b ,w n ) ] 



g (<V%) p ^ w b^ 


p(w n )dw b dw n . 


(B.2) 


The probability density function p(w b ) describing w fe in Aw is 
(Aw)" 1 [B.2]. The number of acoustic modes with frequencies w n 
below the band Aw for a rectangular volume is approximately 


N n<Aw 


Kw 3 

c 


where w is the band center frequency and K is a constant of 
c 

proportionality. 


The probability distribution function P( w n ) ? or acoustic modes 
located in the frequency band Aw may be calculated as [B.2] 
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where Aw = 2ew c and e<<l. Subscripts u and l refer to the upper 
and lower limits of the frequency band Aw. 


By definition, the probability density of w n [B.2] is 
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Consider integration with respect to w : 
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The function contained in { } brackets varies between tt/2 and tt , 

as a) varies between w (1 - e) and w (1 + e) , with values generally 
n c o 

close to tt. Consequently, 
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APPENDIX C 


rev 

RELATIONSHIP BETWEEN Im[I rr (u)) ] AND f 2 (n,r), j 2 (w) 

The expected value of the joint intermodal coefficient, which 
expresses the coupling between structural modes and modes of 
the cavity, is defined by Eq. (20) as 

E[I rr (u>) ] = e[£ ^ ( ^ E - k 2 ) - f 2 (n,r)j (C.l 

where f(n,r) = j J $ n (x)ip r ( x)dx is the coupling factor between 

the structure and the cavity. 

The imaginary part of E[I r:r (w)] is required for calculation of 
the real power flowing inwardly [see Eq. (17)]; then 


Im{E[I rr (o)) ]) = E {lm[I rr ( w) ] } = e\~ 


e k k 

n n n ^ f2(n 
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n i n 
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r A 2 c ? 
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n n * 
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- ♦ -i 


'( c . 2 ) 


using Eqs. (22) and (23). 


a 2 c ? 


Vo^ <e n f 2 (n.r)>n n E[G(«., t o rl )] 


(C.3) 
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where E[G(oo,w n )] is defined and derived in Section (3.5.1) 

[Eqs. (39) and (40)]: then 

ttA 2 c? 

E{Im[I rr (w) ] > = -gv— <e n f 2 (n,r)> n n (w). (C.4) 

■jo •jo 

Alternatively, an approximation for Im[I (to)] can be obtained in 
rev 

terms of j (w), the joint acceptance function for a reverberant 
pressure field. This approximation is obtained by analogy with 
results from [C.l] and [C.2], as discussed in [C . 3 ] to give 


, rev 

E(Im[I rr (m)]} = 77 77 A 2 j 2 (w). 


(C. 5) 


Then, comparing Eqs. (C.4) and (C.5), one finds that 
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rev 
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APPENDIX D 


FINAL SIDEWALL DESIGN UTILIZING ADD-ON TREATMENTS 

It is demonstrated in Section 7 that the only add-on treatment 
capable of providing the noise reductions required to achieve 
the goal of an A-weighted sound level less than 80 dB is a 
double-wall- system with a limp-mass inner wall. In addition, 
the calculations in Section 8 do not identify any advanced 
concept which can achieve the noise reduction predicted for 
the double wall. Thus it is appropriate to summarize the final 
double-wall designs that are predicted to provide cabin noise 
levels which do not exceed 80 dB(A) in the three study airplanes. 

For the purposes of the investigation, three baseline air- 
craft were chosen on the basis of existing technology and 
conventional fuselage construction. To some extent the choice 
of airplane influences the amount of treatment required on the 
sidewall, and this influence will be reflected in the details 
of the three double-wall systems described in this Appendix. 
Furthermore, certain simplifying assumptions were made, with 
the result that noise transmission through the floor and 
windows was neglected, as was structureborne transmission 
through the sidewall itself. It is probable that, in practice, 
these paths will require some noise control treatment, or the 
performance of the double-wall sidewall will be compromised. 

The construction of the three sidewall systems is described 
below in terms of the different components. The sidewall is 
shown schematically in Figure D.l. 
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Fuselage Structure : - 


In all the study aircraft the fuselage structure Is conventional 
aluminum skln-stringer-frame construction. For two of the 
study aircraft (wlde-body and small-diameter) the skin thick- 
ness is constant along the cabin length, while for the narrow- 
body design the skin thickness Is allowed to vary. Details of 
the fuselage structure are given In Sections 4 . 5 and 5*7 9 
Figures 17 to 19, and Tables 3 and 7. The structures have no 
additional damping treatment although it might be advisable to 
treat the frames and stringers in order to minimize structure- 
borne noise transmission through flanking paths such as the 
floor and the trim panel mounts. 

Add-On Insulation: - 

The add-on insulation consists of a limp-mass trim panel, 
which forms the second wall of the double wall system, and 
acoustic insulating material which fills the cavity between 
the fuselage skin and the trim panel. 

The trim panel is installed at a distance of at least 12.7min 
from the fuselage skin panel, has limp-mass characteristics 
such as those of lead-impregnated vinyl, and has high internal 
damping, with a loss factor of about 0,1. The damping could 
be provided by the use of visco— elastic material. A framework 
may be necessary to support the trim panels, but the framework 
should be isolated from the fuselage structure so that there 
is negligible structureborne noise transmission. 

The acoustic insulation in the cavity is a light-weight, porous 
material such as glass fiber wool. The flow resistivity should 
be approximately 50,000 mks rayls/m and the bulk density about 
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9.6 kg/m 3 . For a thickness of 127 mm, the surface density 
is 1.2 kg/m 2 . 

The variation in surface density of the trim panel plus porous 
material along the length of the cabin is shown in Figure 56. 
Since there is a uniform distribution (1.2 kg/m 2 ) of the porous 
material, the longitudinal variation is due solely to changes 
in trim panel surface density. There is a uniform circumferen- 
tial distribution of porous material and trim panel surface 
densities. The maximum values of the trim panel surface density 
are approximately 36 kg/m 2 for the wide-body airplane, 14 kg/m 2 
for the narrow-body airplane, and 26 kg/m 2 for the small- 
diameter airplane. The maximum trim panel weight for the 
narrow-body airplane is lower than those for the other two 
aircraft because of the heavier skin panel weight in the region 
of maximum exterior sound pressure level. 

Good isolation of the trim panel from the fuselage structure 
(skin, stringers, frames and floor) is essential if the acous- 
tic performance of the add-on sidewall is not to be degraded. 

The design of vibration isolation mounts was not included in 
the present study and work is required on this topic. It is 
possible that the low transmissibilities required at the 
fundamental blade passage frequencies will require the use of 
constrained damping material on the frames and stringers, 
to increase the total damping loss factor to lie in the range 
0.05 to 0.1. 

Interior Furnishing s : - 

The acoustic absorption coefficients for the interior surfaces 
and furnishings should be equal to, or greater than, the values 
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shown in Figure 22. To achieve these values, the carpets, 
seats and the surface of the trim panels will have to be 
at least as acoustically absorptive as are the furnishings 
in current turbofan aircraft. 
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APPENDIX E 
LIST OF SYMBOLS 
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Structural area 

Cross-sectional area of ring frames and 
stringers, respectively 

Areas of absorbing walls 

Pressure amplitude decay rate parameter 

Cavity dimensions in x,y,z directions (Figure 4(a)) 

Propagation constant in porous material b - ot+j 3 

Co— spectral density of exterior blocked pressure 

Correlation decay parameters for exterior 
pressure field 

Speed of sound in air in the interior and 
exterior media 

Compressional wavespeed in material 

Propeller diameter; or bending rigidity of 
structure 

Bending rigidities of structure and skin 
Cross rigidity of structure 
= Gh 3 /12 

Elastic modulus for ring frame, stringer and 
skin, respectively 

Expected value of function g( u b , u n ) 

Expected value of the acoustic power flow into 
cavity at frequency centered in band Aw 


f (n,r) 


Internal coupling factor for structure and cavity 


G 


Modulus of rigidity for shear stresses of 
structure 


Gp(x|x';<Jj) Green’s function for the airplane cabin interior 

G°(x|x';w) Green’s function for the exterior field 

Jr 


V a s 


s'V’rP 

H 


Torsion modulus of ring frames and stringer, 
respectively 

Modal admittance function at frequency 
Cabin height 
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Skin thickness and effective thickness respectively 



Acoustic intensity at 
= <P?>/ P j c. 


1 1 
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Moment of inertia of ring frame and stringer, 
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Internal acoustic coupling factors for structural 
modes (Section 3.2.1) 


.rs .rr 
j , J 


External acoustic coupling factors for structural 
modes (Section 3.2.1) 


J ,J 
r * s 


Polar moment of inertia of ring frames and 
stringers respectively 
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tudinal §nd lateral directions for progressive 
wave excitation 

Joint acceptance for reverberant field excitation 
Complex compressibility of the porous material 
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x J y 


Wavenumber = w/c (k=k ,k ) 

o x* y 

Structural wavenumber, where k = mTr/L 

in x 

Complex eigenvalue for cavity 
Real part of k n 

Dimensions of rectangular curved structural 
element 
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x* y 


M r 

m,n 


Sidewall cavity depth 

Subpanel element dimensions 

Frame Pitch 

Stringer Pitch 

Modal mass of structure 

Number of axial halfwaves and number of 
circumferential full waves in structure 


N 


n 


NR 


Number of modes in band 

Noise reduction = 10 log (<?*>/<??>) 


n r’ 


n 


n 


Structural and acoustic modal densities, 
respectively 


P°,P^,P Fourier transform of external, interior and 

1 1 bll blocked pressures of truncated signal of 

length T, respectively 

<p?> Space-averaged, band-limited, mean square pressure 

1 in cabin volume 


<p 2 > Space-averaged band limited mean square pressure 

n in volume V, for nth mode 

p 3 q Mode orders; also indices in image array formula- 
tion 


qpCxjx' ) 


Spatial distribution function for power spectral 
density of the exterior blocked pressure field 


R 


Radius of curvature; or distance between source 
and observer in cabin 


R,3 


Ring frame and stringer bending parameters, 
respectively 
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R 

R 


rad 
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rad 
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Radiation resistance 

Internally-looking radiation resistance 
Externally-looking radiation resistance 
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Plow resistivity of the porous material 


R p (x,x' jt) 
r 


Average cross-correlation of the blocked 
pressure over the exterior 

Distance between propeller tip and fuselage 
sidewall 


(Sa) Absorption, in cabin interior, of absorbing 

area S 


V“> 


Power spectral density of exterior blocked 
pressure 


Sp( x,x ;w) 


Cross power spectral density of exterior blocked 
pressure 



One-sided power spectral density of exterior 
blocked pressure at reference location Xo 


S w (x,x' ;w; Cross-power spectral density of panel dis- 
placement 


Sp 1 v(x,x' ;m) Cross-power spectral density of pressure at x 
and velocity at x' 


V 
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w T 
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rad^ ' 


[W in (w)] Am 
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w 


int 

abs 


Average trace wavespeed of exterior pressure field 
Volume of airplane cabin 

Net power flowing through structure; or cabin 
width 

Fourier transform of displacement truncated 
signal of length T 

Spectral density of power radiated into cabin 

Power inflow to acoustic structural modes 
resonant in band Am 

Power radiated from jth structural element 
Power absorbed on cabin surfaces 
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z , z 
V* s 


= Xo, Yo, is the location of peak pressure on 
the structural element 

Porosity of porous material 

Structural receptance for mode r 

Characteristic impedance of porous material, 

Z q = -jKb/wy , see Eq. (65), Section 3.9.2 

Surface specific acoustic impedance ( = p/v) 

Infinite panel impedance, Eq. (68) 

Distance from skin middle surface to centroids 
of ring frame and stringer 


VVV a z>VV a f 


V e q’ £ r 


Absorption coefficients of surfaces 

Surface admittances, 8=5 - i a • Also complex 
wavenumber in porous material 

Generalized blocked force 

Defined in Eq. ( 26) 

Frequency bandwidth 

Pressure differential across cabin fuselage 
= nL /imtR 

X 

Mode participation factors, (= 1 jP,q, r=0 : = 2; p,q,r>0) 
= e e e 

P q r 

Also relative error - Eq. (46) 

Separation distance in y direction U = y' - y) 
Structural mode loss factor 
Average value of over band Aw 
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a rad 
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? (x) 
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Acoustic mode loss factor 

Average value of r) over band 

Radiation loss factor, = pca rad /ioy 

Loss factor associated with structural dissipation 
mechanisms 

Radius of gyration of panel, structure 
Imaginary part of 

Parameters used in frequency equation for 
simply-supported cylinder 

Skin critical wavelength 

Average mass surface density of structure 

Average mass surface densities for two walls of 
a double wall structure 

Structural resonance frequency, non-dimensionalized 
by the ring frequency 

Poisson's ratio of structure, for x,y axes 

Separation distance in x direction, £ = x' - x 

Normalized surface conductances for surfaces 
of cavity 

Ambient densities in the interior and exterior 
media 

Variance of function G 

Normalized susceptances of surfaces of cavity 

in t* 

Radiation efficiency, a rad = R ra d/ p i c i A 

Bare fuselage transmission coefficient 

Transmission coefficient of add-on sidewall 
treatment 

Complex mode shape for cavity 


tl (3) 

ip r (x) 


CO 

w b 


CO 


CO 


O 


W. 


V 


0 ) 


n 


< > 

r 


Real part of complex mode shape <j> n (x) of cavity 

Mode shape, or eigenfunction, of the rth mode 
of the structure 

Angular frequency 

= bw 0 , bth harmonic of the blade passage 
frequency co q 

Center frequency in band Am 

Angular blade passage frequency, 2irf o 

Structural mode resonance frequency; also 
fuselage ring frequency 

Acoustic mode resonance frequency 

indicate averages over modes in band 
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